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ABSTRACT 


This  report  explains  the  fundamental  theory  and  equations 
used  in  writing  a  set  of  software  applications  which  predict 
antenna  radiation  parameters.  Each  application  predicts  the 
radiation  characteristics  of  a  particular  type  of  antenna  over 
a  planar  surface  which  serves  as  a  model  of  either  earth  or 
seawater.  The  radiation  parameter  predictions  are  based 
solely  rn  an  antenna's  physical  dimensions,  the  properties  of 
the  uv.uerlying  surface,  and  electromagnetic  theory.  Existing 
electric  field  equations  provide  the  basis  for  radiation 
parameter  predictions,  and  the  accuracy  of  the  predicted 
radiation  parameters  is  totally  dependent  upon  the  extent  to 
which  the  electric  field  equations  used  realistically  model 
the  actual  radiated  electromagnetic  fields  of  the  antennas. 

In  addition  to  a  review  of  applicable  electromagnetic 
field  theory,  this  report  is  also  intended  to  be  a  user's 
guide  for  the  corresponding  computer  applications.  The 
appendices  contain  computer  hardcopies  of  sample  calculations 
for  several  antenna  types  and  remarks  regarding  the  conformity 
of  predicted  radiation  parameters  to  expectations.  Radiation 
parameters  computed  thus  far  are  consistent  with  expectations 
based  on  other  computational  programs  and  empirical 
measurements . 
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I.  BACKGROUND  AND  PURPOSE 


This  report  and  associated  computer  software  applications 
are  submitted  in  fulfillment  of  the  thesis  requirements  for 
the  degree  of  Master  of  Science,  Aeronautical  Engineering  from 
the  Naval  Postgraduate  School  in  Monterey,  CA.  The  thesis 
text  explains  the  fundamental  theory  and  equations  used  in 
writing  a  set  of  software  applications  which  predict  an 
antenna's  radiation  parameters.  Each  application  predicts  the 
radiation  characteristics  of  a  given  type  of  antenna  over  a 
planar  surface  which  serves  as  a  model  of  either  earth  or 
seawater.  The  thesis  requirement  was  generated  by  a  request 
from  the  Naval  Maritime  Intelligence  Center  (NAVMARINTCEN) . 
At  the  request  of  NAVMARINTCEN,  Mathcad  mathematical  software, 
DOS  version  3.1,  was  used  to  write  the  applications.  The 
NAVMARINTCEN  request  was  extensive  enough  to  be  tasked  as  two 
separate  thesis  topics.  This  thesis  fulfills  the  first  half 
of  the  NAVMARINTCEN  requirements.  Mathcad  applications  are 
complete  for  the  following  types  of  antennas: 

1.  Vertical  Monopole  Antenna 

2.  Elevated  Vertical  Dipole  Antenna 

3.  Elevated  Horizontal  Dipole  Antenna 

4.  Arbitrarily  Oriented  Dipole  Antenna 

5.  Inverted  L  Antenna 

6.  Long-Wire  Antenna 

7.  Terminated  Sloping  V  Antenna 

8.  Side-Loaded  Vertical  Half  Rhombic  Antenna 

9.  Terminated  Sloping  or  Horizontal  Rhombic  Antenna 

10.  Sloping  Double  Rhomboid  Antenna 

11.  Vertically  Polarized  Log -Periodic  Dipole  Array 

12.  Horizontally  Polarized  Log-Periodic  Dipole  Array 

13.  Horizontal  Yagi-Uda  Array 
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II.  INTRODUCTION 


New  American  Heritage  Dictionary  defines  antenna  as  'an 
apparatus  for  sending  and  receiving  electromagnetic  waves'. 
An  electromagnetic  wave  is  a  time-varying,  self -propagating, 
interrelated  combination  of  electric  and  magnetic  fields. 
Antennas  radiate  electromagnetic  fields  as  a  function  of  their 
time-varying  surface  charge  density  and  surface  current 
density  distributions.  There  are  many  different  types  of 
antennas,  each  with  its  own  set  of  radiation  characteristics 
or  radiation  parameters.  These  radiation  characteristics  are 
related  to  an  antenna's  radiated  electromagnetic  field 
distribution,  and  they  determine  the  useful  applications  for 
a  particular  antenna  design. 

Current  personal  computers  with  high-speed  processors  can 
quickly  and  accurately  predict  antenna  radiation  parameters 
from  electromagnetic  theory.  The  Mathcad  routines  described 
in  this  report  predict  an  antenna's  radiation  parameters  based 
solely  on  its  physical  dimensions  and  electromagnetic  theory. 
Accurately  predicted  radiation  characteristics  provide 
intelligence  analysts  with  a  reliable  estimate  of  an  antenna's 
capabilities,  limitations,  and  vulnerabilities.  Chapters  III 
and  IV  are  a  review  of  applicable  electromagnetic  theory,  and 
the'  remaining  chapters  describe  the  calculations  of  the 
associated  Mathcad  computer  code,  and  are  intended  to  be  a 
user's  guide  for  the  associated  Mathcad  applications. 
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III.  FUNDAMENTAL  ELECT~ “MAGNETIC  THEORY 
A.  MAXWELL'S  EQUATIONS  AND  HELMHOLTZ'S  EQUATIONS 

Maxwell's  equations  are  the  basis  for  electromagnetic 
field  calculations.  Their  differential  form  is  given  by  [Ref 
1:  pp.  321-325] 

V  X  E  =  -\ (3.1a)  V  •  E  =  -£  (3.1c) 

dt  e 

V  X  If  =  J- +  e-f^  (3.1b)  V  •  H  =  0  (3.  id) 

ot 

Equation  3.1a  is  Faraday's  Law,  3.1b  is  Ampere's  Circuital 
Law,  3.1c  is  Gauss's  Law,  and  3. Id  postulates  the  inexistence 
of  magnetic  charge. 

Maxwell's  equations  dictate  the  relationship  between  the 
electric  and  magnetic  field  intensities  (E  and  H)  and  the 
charge/current  density  source  distributions  (p  and  J)  (i.e., 
the  antenna).  Maxwell's  equations  are  often  applied  in  a 
source-free,  current-free  region,  and  an  e*"'  time  dependence 
is  assumed.  The  time-harmonic,  free  space  Maxwell's  equations 
are  [Ref  1:  p.  340] 

V  X  E  =  -jo>y0ff  (3.2a)  V  •  £  =  0  (3.2c) 

V  X  H  =  jwe0E  (3.2  b)  V  •  H  =  0  (3. 2d) 

where 

eQ=free  space  permittivity =— i-  •  io-9  /  ^radsv 

0  3 6 tv  \  meter  ) 

I xQ=free  space  permeability=  4n  •  10*7 
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Rather  than  solve  Maxwell's  equations  directly,  the 
method  of  vector  potentials  is  often  used  where  the  scalar 
potential  (V)  and  magnetic  vector  potential  (A)  can  be 
obtained  from  the  nonhomogeneous  Helmholtz  equations 

V*V  +  $2V  =  (3.3)  V2*  ♦  P2A  =  (3.4) 
eo 

where  H=(1/m0)  (V  X  A)  [Ref  1:  pp.  338-340]. 

The  solutions  to  Helmholtz's  equations  are  the  retarded 
scalar  potential  and  the  retarded  vector  potential .  The  value 
of  V  and  A  at  some  distance  (R)  from  the  source  depends  on  the 
source's  charge/current  density  at  an  earlier  time  (t-R/c) , 
where  c=*l/VV0e  0  is  equal  to  the  speed  of  light  in  free  space. 
The  delay  (R/c)  is  the  time  required  for  electromagnetic. waves 
to  propagate  through  a  distance,  R,  in  free  space  from  source 
to  observation  point.  The  solutions  to  the  nonhomogeneous 
Helmholtz  equations  are  [Ref  1:  pp. 338 -340] 


1  r  r\a~j$\8~8'  I  . 

V{R)=—±—l  Pe _  - - dv' 

4 ite0Jv*  I  R-R'  I 


A(R)=i^ 

4  it  Jv 


L 


r-r' 


■  dv' 


(3.5) 

(3.6) 


where  the  wavenumber  (0)  is  c^t0e  0  •  in  free  space,  R  is  th<S 
observation  vector,  and  R'  is  the  source  vector.  The  *?ave-\ 
number  is  most  often  expressed  as  w/c  =  2i r/X,  where  \  is  the 
wavelength.  The  Mathcad  applications  use  the  expression  2n/\ 
to  calculate  0. 
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The  Mathcad  applications  described  herein  assume  that  the 
antenna  under  analysis  is  radiating  in  free  space  above  a 
defined  ground  plane.  The  free  space  Maxwell's  equations  and 
nonhomogeneous  Helmholtz  equations  (using  the  free  space 
wavenumber)  provide  the  theoretical  basis  for  deriving  the 
equations  for  an  antenna's  radiated  electromagnetic  fields. 
The  radiated  E  and  H  fields  can  be  calculated  directly  from 
their  source's  charge  and  current  densities  using  Maxwell's 
equations.  However,  as  previously  mentioned,  it  is  often 
simpler  to  find  the  magnetic  vector  potential  using 
Helmholtz's  equations.  Once  the  vector  potential  is 
determined,  the  electric  and  magnetic  fields  are  calculated 
from  k=(l/nQ)(V  X  A)  and  E=(l/jwe0)(V  x  1)  [Ref  1:  pp.  338- 
341] .  Existing  theoretical  equations  for  the  radiated 
electric  fields  of  many  antenna  types  have  been  derived  in 
this  manner.  These  equations  are  used  in  the  Mathcad 
applications  to  predict  an  antenna's  radiated  electric  field 
distribution.  Other  radiation  par  meters  are  then  calculated 
on  the  basis  of  the  predicted  electric  field  distribution. 

B.  ELECTROMAGNETIC  FIELD  REGIONS 

The  space  surrounding  an  antenna  is  divided  into  three 
regions:  the  reactive  near- field,  radiating  near- field,  and 
far-field  [Ref  2:  pp.  22-24].  The  reactive  near-field 
occupies  the  space  immediately  surrounding  the  antenna  out  to 
a  radius  of  about  0 . 62  •  (D3/X) 'S  where  X  is  the  wavelength  and 
D  is  the  maximum  dimension  of  the  antenna.  Reactive 
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electromagnetic  fields  dominate  in  this  region.  The  radiating 
near-fiel  ccupies  the  space  from  the  boundary  of  the 
reactive  field  out  to  a  radius  of  about  2-(D2/X).  Radiation 
fields  dominate  in  this  region,  and  radiation  patterns  are 
often  a  function  of  both  radial  and  angular  coordinates.  The 
far- field  region  occupies  all  space  outside  the  radiating 
near- field  which  meets  two  additional  far- field  requirements: 
27rR/X  >>  1  and  R  >  5D.  Radiation  fields  also  dominate  in  the 
far- field,  but  field  components  are  primarily  transverse,  and 
radiation  patterns  are  normally  a  function  of  directional 
variables  only.  The  Mathcad  antenna  applications  are  valid 
only  in  the  far -field  region  of  a  given  antenna. 

C.  RADIATION  PATTERNS 

A  radiation  pattern  is  a  three-dimensional  representation 
of  an  antenna's  radiated  electromagnetic  field  distribution  or 
power  distribution  at  a  fixed  distance  from  the  antenna. 
Because  it  is  difficult  to  depict  three-dimensional  images, 
the  patterns  are  most  often  plotted  in  a  defined  vertical  or 
horizontal  plane.  The  Mathcad  antenna  applications  plot  two- 
dimensional  radiation  patterns  in  polar  coordinates  depicting 
the  far-field  electric  field  distribution.  Patterns  are 
plotted  in  the  <f>=0  and  <f>= tt/2  vertical  planes,  and  also  in  a 
horizontal  plane  at  an  elevation  angle  set  by  an  index  from 
the  applicable  'Elevation  Angle  Index  Table'  (Table  31  or 
3.2).  Field  magnitudes  are  normalized  with  respect  to  the 
maximum  radiated  electric  field  intensity.  The  magnitude  to 
which  each  pattern  is  normalized  is  displayed  below  its  plot. 
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ELEVATION  ANGLE  INDEX  TABLE 
FOR  LOG- PERIODIC  AND  YAGI-UDA  ARRAYS 


ELEVATION 

INDEX 

ELEVATION 

INDEX 

314 

45.71° 

235 

2.28° 

311 

48.00° 

231 

4.57° 

307 

50.28° 

227 

6.86° 

303 

52.57° 

223 

9.14° 

299 

54.86° 

219 

11.43° 

295 

57.14° 

215 

13.71° 

291 

59.43° 

211 

16.00° 

287 

61.71° 

207 

18.28° 

283 

64.00° 

203 

20.57° 

279 

66.29° 

199 

22.86° 

275 

68.57° 

195 

25.14° 

271 

70.86° 

191 

27.43° 

267 

73.14° 

187 

mmm 

263 

75.43° 

183 

32.00° 

259 

77.71° 

179 

34.29° 

255 

80.00° 

175  “ 

36.57° 

251 

82.29° 

171 

38.86° 

247 

84.57° 

167 

41.14° 

243 

86.86° 

163 

43.43° 

239 

89.14° 

159 

8 


\ 


The  Mathcad  applications  plot  the  far-field  radiation 
patterns  for  an  antenna's  space  wave,  surface  wave,  and  total 
radiated  electric  field.  The  space  wave  is  comprised  of  the 
electromagnetic  waves  which  propagate  directly  from  antenna  to 
the  field  point  and  electromagnetic  waves  reflected  from  the 
surface  below  the  antenna.  The  surface  wave  is  composed  of 
electromagnetic  waves  ducted  along  the  surface,  much  like 
waveguide  propagation.  The  total  field  is  composed  of  the 
space  wave  and  surface  wave  combined. 

Space  wave  patterns  typically  vary  with  directional 
spherical  coordinates,  6  and  <p,  but  they  are  predominately 
invariant  with  distance  from  the  antenna.  The  surface  wave 
attenuates  exponentially  with  distance,  and  far-field  surface 
wave  patterns  are  a  function  of  distance  from  the  antenna. 
Total  field  patterns  are  still  predominately  invariant  with 
distance,  because  the  maximum  surface  wave  intensity  is 
usually  many  orders  of  magnitude  less  than  the  maximum  space 
wave  intensity.  Even  at  distances  which  just  meet  far-field 

requirements,  the  surface  wave  contribution  to  the  total  field 

I 

is  barely  significant  for  parameter  values  of  interest. 

D.  POWER  DENSITY,  POYNTING  VICTOR  AND  RADIATION  RESISTANCE 
The  instantaneous  real  power  flux  density  of  an  electro¬ 
magnetic  wave  is  called  the  Poynting  vector  (<?) , 

P  (x.y.z)  =  EX  B  =  3te{£U,y,z)  ^u£}  X  &e{HU,y,  z)  ejat)  (3.7) 

=  -i-8ie{ff  X  H*}  +  -|3t e{E  X  He2^wt} 
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The  time  invariant  term  on  the  right  hand  side  of  equation  3.7 
is  the  time  average  power  flux  density  or  average  Poynting 
vector.  Far- field  E  and  H  fields  are  related  by  the  simple 
expression  |E|  =  (1/ij0)  *  |H|  ,  where  »?o=vVo/e  o  =  i207r  is  the 
intrinsic  impedance  of  free  space.  The  average  Poynting 
vector  solely  in  terms  of  the  radiated  electric  field 
intensity  is  given  by  [Ref  1:  pp.  382-385]  . 

Vavfr(x.y,  z)  =±Sle{E  X  H'}  =  ±8le^E  X  =  (3.8) 

The  total  power  flax  density  is  composed  of  both  the 
Poynting  vector  and  the  imaginary  part  of  (E  X  H*)  which 
represents  the  reactive  power  of  the  radiated  electromagnetic 
fields.  Far- field  power  flux  density  is  predominately  real 
and  is  approximately  equal  to  the  average  Poynting  vector.  An 
antenna's  total  average  radiated  power  is  calculated  by 
integrating  the  poynting  vector  over  a  Gaussian  surface  which 
surrounds  the  antenna  in  the  far- field  [Ref  1:  pp.  382-386] . 

Pr ad  =  Pavg  «  fs  «W  dB  =  jfs*e{E  X  H')-  ds  (3.9) 

The  Mathcad  applications  use  equations  3.8  and  3.9  for  power 
calculations  by  integrating  over  the  hemisphere  which  encloses 
the  antenna  with  the  ground  plane  as  the  lower  boundary. 

Radiation  resistance  relates  an  antenna's  total  average 
radiated  power  to  its  peak  input  current.  If  a  DC  current 
equal  to  the  rms  input  current  were  applied  to  a  resistance 
equal  to  the  radiation  resistance,  the  dissipated  power  would 


equal  the  antenna's  total  average  radiated  power  [Ref  3:  pp. 
47-48] .  For  the  Mathcad  antenna  applications,  the  radiation 
resistance  is  defined  in  terms  of  the  total  average  radiated 
power  and  input  current  by 


^zad~ 


2i>,ad 


or 


^ zad~ 


Kzad 


A  sinusoidal  input  current  {or  voltage  response  across  the 
input  terminals  for  Yagi-Uda  Arrays/  with  a  maximum  of  unity 
is  assumed  in  the  Mathcad  applications. 

E.  DIRECTIONAL  GAIN,  DIRECTIVITY  AND  POWER  GAIN 


An  isotropic  radiator  is  a  hypothetical  antenna  which 
radiates  electromagnetic  energy  equally  in  all  directions. 
Its  radiation  pattern  in  any  plane  passing  through  the  source 
is  a  circle.  All  real  antennas  are  directional  radiators 
which  radiate  or  receive  electromagnetic  energy  more 
effectively  in  some  directions  compared  to  others.  Some 
antennas  radiate  an  omnidirectional  pattern  which  is 
directional  in  elevation,  but  non- directional  in  azimuth. 

Directive  gain  (Df)  is  the  ratio  of  an  antenna's  Poynting 
vector  in  a  given  direction,'  divided  by  the  power  flux  density 
of  an  isotropic  source  with  equal  total  average  radiated 
power.  Directivity  (D0)  is  the  maximum  value  of  the  directive 
gain,  and  it  is  the  primary  measure  of  an  antenna's 
directional  properties.  The  power  flux  density  of  an 
isotropic  radiator  is  equal  to  the  total  average  radiated 
power  divided  by  the  surface  area  of  a  sphere  at  a  given 
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distance  from  the  antenna  (47rR:)  .  The  expression  for 
directive  gain  is  [Ref  2:  pp.  29-31] 

zye,<t>)  =  9  =  ±*?221*.>  4> )  o.id 

Viso  °rad 

The  Mathcad  applications  calculate  only  directivity  using 

max  )  _  47t^2p(emax<4>max)  (312) 

^iso  ^rad 

Power  gain  accounts  for  an  antenna's  radiation  efficiency 
in  addition  to  its  directional  properties.  Power  gain  (Gd)  is 
defined  as  47rR2  times  the  ratio  of  the  Poynting  vector 
magnitude  in  a  given  direction  to  the  total  input  power  to  the 
antenna.  Total  input  power  includes  losses  due  to  ohmic  power 
dissipation  and  impedance  mismatches  between  transmission  line 
and  antenna.  The  ratio  of  total  average  radiated  power  to  the 
total  input  power  is  the  radiation  efficiency,  peff.  The 
relationship  between  power  gain,  directive  gain,  and  radiation 
efficiency  is  given  by  [Ref  2:  pp.  42-45] 

P8«=J^  (3.13) 

r in 

Gd(e,<|>)=pef^(0rf4))  (3.14) 

Go=G^x(e,4»)=P8^?MJ<(0,4>)=P8f^o  (3.15) 

Since  radiation  efficiency  and  power  gain  account  for 
ohmic  losses  and  impedance  mismatches,  these  parameters  are 
difficult  to  predict  based  on  antenna  dimensions  alone.  For 
these  reasons,  the  Mathcad  applications  do  not  predict  power 
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gain  and  radiation  efficiency.  Gain  calculations  are  limited 
to  directivity  since  it  depends  only  on  total  average  radiated 
power  and  the  radiated  electromagnetic  field  distribution. 
The  maximum  Effective  Isotropic  Radiated  Power  (EIRP)  is  also 
calculated  from  EIRP=Pra<i-D0  [Ref  3:  p.  62]. 

F.  HALF -POWER  BEAMWIDTH 

Radiation  lobes  are  portions  of  the  radiation  pattern 
where  a  local  maximum  field  intensity  is  bounded  on  either 
side  by  relative  minimum  intensities.  The  major  lobe,  or 
main  beam,  is  the  lobe  in  the  direction  of  maximum  radiated 
power  flux  density.  Split  beam  antennas  have  more  than  one 
main  beam.  A  minor  lobe  is  any  lobe  except  the  main  beam. 
Beamwidth  refers  to  the  half -power  beamwidth  of  an  antenna's 
main  beam  (or  beams) ,  and  it  can  be  measured  in  the  vertical 
or  horizontal.  For  antennas  with  a  defined  beamwidth,  the 
power  intensity  drops  off  on  either  side  of  the  point  of 
maximum  radiated  average  power  density  until  it  is  one  half  of 
the  maximum  value.  The  angle  between  the  radials  which 
intersect  the  half  power  points  is  defined  as  the  half -power 
beamwidth.  Main  beams  are  not  necessarily  symmetrical  since 
the  radial  which  intersects  the  point  of  maximum  power  density 
dees  not  necessarily  bisect  the  half -power  radials. 

The  Mathcad  applications  do  not  calculate  half -power 
beamwidth  directly,  but  radiation  patterns  provide  a  good 
estimate.  Since  the  radiation  patterns  depict  the  electric 
field  distribution  normalized  to  unity,  a  magnitude  of 
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approximately  0.7  represents  the  half  power-points  of  the 
pattern.  If  two  radials  are  drawn  from  the  origin  through  the 
half-power  points,  the  angle  subtended  by  the  two  radials  is 
the  antenna's  half -power  beamwidth. 

G.  POLARIZATION  AND  BANDWIDTH 

Polarization  describes  the  time  varying  properties  of  an 
antenna's  radiated  electric  field  at  a  set  spatial  coordinate 
as  electromagnetic  waves  propagate  outward  from  the  antenna. 
Signals  are  usually  transmitted  with  linear,  elliptical,  or 
circular  polarization.  Linear  polarization  is  where  the  field 
vector  at  the  spatial  coordinate  is  directed  along  a  fixed 
linear  path.  An  elliptically  polarized  signal  is  one  in  which 
the  tip  of  the  electric  field  vector  traces  an  ellipse  around 
the  fixed  coordinate  as  the  wave  propagates  outward.  Circular 
polarization  is  a  special  case  of  elliptical  polarization 
where  the  magnitude  of  the  electric  field  remains  constant  as 
it  rotates  about  the  spatial  coordinate.  Clockwise  rotation 
of  the  electric  field  vector  is  right-hand  polarization,  and 
counter-clockwise  rotation  is  left-hand  polarization.  The 
Mathcad  applications  do  not  quantitatively  predict  the 
polarization  at  a  selected  spatial  coordinate  because  of  the 
computational  intensity  of  such  a  calculation.  Only  the 
predominant  type  of  polarization  expected  is  addressed  briefly 
in  the  introductory  remarks  of  each  application. 

The  bandwidth  of  an  antenna  is  the  range  of  frequencies 
over  which  its  radiation  characteristics  meet  or  exceed  a 
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specified  acceptable  performance  standard.  Because  radiation 
characteristics  of  dissimilar  antennas  vary  differently  with 
frequency,  there  are  no  set  rules  or  guidelines  which  define 
bandwidth  for  all  antennas  in  general.  Bandwidth  standards 
vary  depending  on  antenna  design  and  intended  application. 
Therefore,  the  Mathcad  applications  do  not  predict  operational 
bandwidth.  Instead,  multiple  frequencies  can  be  input  for  a 
given  configuration  to  determine  the  bandwidth  for  whatever 
standards  are  used  to  define  acceptable  performance. 

H.  EFFECTIVE  LENGTH  AND  EFFECTIVE  AREA 

The  effective  area,  or  effective  aperture,  gives  the 
power  delivered  to  the  antenna  load  when  multiplied  by  the 
incident  average  power  flux  density.  It  is  defined  by  [Ref  2: 
pp  59-63] 

(3.16) 

The  maximum  effective  area  is  related  to  directivity  for  all 
antennas  by  [Ref  4:  pp.  46-47] 

.  <3.17)  . - 

This  is  a  theoretical  maximum  effective  area  since  directivity 
does  not  account  for  radiation  efficiency.  Like  directive 
gain  and  power  gain,  actual  effective  area  and  theoretical 
effective  area  are  related  by  radiation  efficiency 

x2g 

A*rt-=“4^-=Pe«A»ax  (3.18) 
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The  maximum  theoretical  effective  area  is  a  good  approximation 
of  actual  effective  area  for  frequencies  within  the  bandwidth 
of  a  well  designed  antenna  with  a  high  radiation  efficiency. 

Effective  length  is  defined  as  the  ratio  of  the  voltage 
induced  across  the  antenna  terminals  by  the  incident  electric 
field  and  is  given  by  [Ref  4:  pp.  40-42] 


laff= 


(3.19) 


Effective  length  is  related  to  effective  area  by 


W2 


M 


■^radAnax 


(3.20) 


Substituting  equation  3.17  into  equation  3.20,  we  get  the 
maximum  theoretical  effective  length  in  terms  of  wavelength, 
directivity,  radiation  resistance  and  intrinsic  impedance  of 
free  space  as 


-^max=2o 


471  tlo 


(3.21) 


Finally,  the  actual  effective  length  is  related  to  the 
theoretical  effective  length  by 

max  (3.22) 

Effective  area  and  effective  length  are  actually  a 
function  of  directional  coordinates  with  respect  to  the 
iource.  The  equations  given  here  are  for  maximum  values  of 
effective  area  and  effective  length  along  the  radial  of 
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maximum  directive  gain  (directivity)  .  While  effective  area  is 
a  parameter  more  generally  applied  to  all  antennas,  the 
Mathcad  applications  use  equations  3.17  and  3.21  to  calculate 
both  a  theoretical  maximum  effective  area  and  effective 
length. ■  If  an  antenna's  radiation  efficiency  for,  a  given 
frequency  is  known,  the  actual  effective  area  and  length  can 
be  calculated  by  multiplying  the  theoretical  values  by  the 
radiation  efficiency. 

The  next  chapter  on  the  dipole  antenna  explains  in  detail 
the  derivation  of  the  radiated  electric  field  equations  for  a 
vertical  dipole  antenna.  The  elevated  vertical  dipole  Mathcad 
application  uses  this  equation  to  obtain  the  electric  field 
distribution  and  total  average  radiated  power  (P„d)  based  on 
an  assumed  sinusoidal  current  distribution  on  the  antenna  with 
a  maximum  of  unity.  Total  average  radiated  power  is  then  used 
to  predict  the  radiation  parameters  discussed  in  the  preceding 
paragraphs.  In  follow-on  chapters,  only  the  final  electric 
field  equations  are  presented  because  the  derivation  is 
similar  for  all  types  of  antennas.  They  are  primarily  summary 
chapters  which  describe  the  calculations  of  each  individual 
application,  but  any  new  or  important  concepts  not  presented 
in  previous  chapters  will  be  discussed  as  required. 


IV.  THE  VERTICAL  DIPOLE  ANTENNA 
A.  TEE  ELEMENTAL  DIPOLE 

The  current  element,  or  elemental  dipole,  is  the 
building  block  for  calculating  antenna  radiation  parameters. 
The  elemental  dipole  is  a  hypothetical  differential  length 
thin-wire  antenna  with  a  sinusoidally  oscillating  surface 
current.  The  antenna  length  (1)  is  many  orders  of  magnitude 
smaller  than  the  wavelength  (X) ,  so  the  instantaneous  current 
is  considered  uniform  over  the  length  of  the  dipole.  Figure 
4 . 1  depicts  the  elemental  dipole  oriented  along  the  z  axis  in 
a  rectangular  coordinate  system  with  its  center  at  the  origin. 


FIGURE  4.1:  Elemental  dipole  spatial  orientation 


The  surface  current  is  given  by 

i  ( t)  =2  cos  ( wt)  =  9te(  Ie'“c}  (4.1) 

and  the  magnetic  vector  potential,  A,  is  found  from  equation 
3 . 6  to  be 


Ho  r  .  Ho Idl[  e-ilR) 


(4.2) 


Primed  coordinates  refer  to  source  points,  R  is  the  distance 
from  the  origin  to  the  field  point,  0-(j/n 0e 0  is  the  free  Space 
wavenumber,  and  the  magnitude  of  the  current  density  (j)  is 
given  by  I  (Ref  2:  pp.  100-102].  The  magnetic  vector 
potential  in  spherical  coordinates  is 


AR-AZ  cos  (0)  = 


ji0  Jdlf  e-j*R) 

cos  (0) 

(4.3a) 

47t  V  *  J 

HoIdJ(  e-iP^ 

(sin  (0) 

(4.3  b) 

4>i  l  R  ) 

A^-Q 

(4.3c) 

The  magnetic  field  intensity,  H,  is  found  from 


if=— V  X  A=a4— i- 
Ho  *  Ho* 


<9  /  v  dil  p 


(4.4) 


P2  sin(0) 

9  47C 


jP*  (jPi?)2 


,  -m 


(4.5) 
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and  the  electric  field  intensity,  E,  is  found  from  Ampere's 
Law  by 


x  H=— — far  ;x  j.  [fllsin  (0)  ]  -a04^  (JOT*)]  (4.6) 

jue0[  ri?sm(8)  dQ  9  B  p  *  I 


£*=~*loP2cos(e) 


.  +  ±  e-j'P* 

(j$R)2  (Jpi?)3. 


£e=*-~rTloP2sin(0) 


I  (Jpi?)  (j’Pi?)2  (Jpi?) 


(4.7a) 


(4 .72?) 


Hr=H6=E^= 0  (4.7  c) 

At  far- field  distances,  the  only  significant  components  are 
the  first  terms  inside  the  parenthesis  on  the  right  hand  side 
of  equations  4.5  and  4.7b.  The  far- field  E  and  H  fields  are 
functions  of  a  single  angular  coordinate  (0)  for  a  vertically 
orientated  element  [Ref  1:  pp.  602-607]. 

Retaining  only  those  terms  that  vary  as  1/R,  we  get  the 
far- field  electric  and  magnetic  field  intensities  for  the 
current  element  as 


^=J’^(‘£T")psin(0)  (4‘8) 

£'®SJ^(‘1F")T,oPsin(0)  (4,9> 


Note  that  the  far- field  relationship  |h|=t?0|e|  holds  true  for 
the  elemental  dipole's  radiated  fields  in  the  far- field 
region. 

The  elemental  dipole's  radiated  average  Poynting  vector 
(power  density)  is  [Ref  1:  607-612] 
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9*  x  •§  ^8  ~~~r-  ^n.p-sin-  (6)  (4.10) 


The  P.oyr.t ir.g  vector  integrated  over  a  Gaussian  surface  in  the 
far- field  yields  the  total  average  radiated  power 

^avi  ds  =  T1  p  :f~  'j*  [sin2(0) ;  Sin(0)  o6c?4>  (4.  ll) 

p  =2it-^^4-Ti.p:  ■— v-9-  -cos  (0)  =2n-(-— ---H-P!|4) 

■  32n*  3  J0  32n2  '  '  3  / 


AlAlll  n  -  G 2  =  JlLALU  n  .  ft 2 

3  32n2  12n 


The  directive  gain  and  directivity  are  found  from  the 
power  density  and  total  radiated  power  as 


(4.12) 


D: (0,  <J>)  ~DJ^  .4>|=1  •  5  or  10  log13  1.5  =  1.76  dB  (4.13) 


Radiation  resistance  is  found  directly  from  Pftd  and  equation 
2.10  as  [Ref  1:  pp.  607-612]. 


p  -Up  -  (Idi)2_  n 2 
pr»d-  2  «r*d-  12n  ’loP 


(4.14) 


-£^  1207t| 

‘  =~'80n2(— -)2 

6  n  \ 

l  X  / 

2  [  l  1  1 

Rt,d=80n2[ 

di\2 

X  j 

(4.15) 
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B.  THE  ELECTRICALLY  SHORT  DIPOLE 


The  derivation  of  the  electric  field  for  an  electrically 
short  dipole  is  the  next  step  in  deriving  the  field  equations 
for  a  practical  finite  length  dipole.  The  electrically  short 
dipole  is  a  center  fed  thin-wire  antenna  oriented  as  in  Figure 
4.1  with  length  greater  than  0.02X,  but  less  than  0.10X.  The 
current  distribution  is  estimated  as  triangular  where 


a,1*,  0*z'*i 

Kx'.y',  z')  =  I(z')  =  ■  ■  (4.16) 

a.Im(l  +  jz'),  -pz'i  0 

where  Im  is  the  current  maximum  at  the  center  feed. 


FIGURE  4.2i  Current  distribution 
for  an  electrically  short  dipole. 


The  dotted  line  in  Figure 
4.2  is  a  more  exact 
representation  of  the 
actual  current  distri¬ 
bution,  but  the  triangular 
estimate  is  a  reasonable 
approximation  for  short 
dipoles. 

The  electrically 


short  dipole  has  a  magnetic  vector  potential  given  by 
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Because  dipole  length  is  small  in  comparison  to  the  far- field 
distance,  R  is  considered  constant  throughout  the  integration 
over  the  length  of  the  dipole.  Note  that  the  magnitude  of  the 
magnetic  vector  potential  for  the  electrically  short  dipole  is 
half  that  of  the  elemental  dipole  for  the  same  current  input. 

The  equations  for  the  radiated  electromagnetic  fields  are 
derived  in  the  same  way  as  for  the  elemental  dipole,  and  the 
final  expressions  for  the  electric  and  magnetic  field 
distributions  of  the  electrically  short  dipole  are  given  by 

fir  ie~jtR 

Se=in° — — sin(0)  (4,18) 

BJ  le'M* 

H^j  sin  (Q)  (4.19) 

Et=E^Hz=H9= 0 

Directive  gain  and  directivity  for  the  electrically  short 
dipole  are  the  same  as  those  of  the  elemental  dipole. The 
magnitudes  of  the  radiated  power  density,  total  radiated 
power,  and  radiation  resistance  for  the  short  dipole  are  given 


by  [Ref  2:  pp.  109-112] 

(4.20) 

P-= 

(4.21) 

,  /  7  i2 

*,.<.20*’  ( j.) 

(4.22) 
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C.  THE  FINITE  LENGTH  DIPOLE 

Finite  length  dipoles  are  center  fed  thin- wire  antennas 
with  negligible  diameter  (diameter  <<  X  and  radius  <<  length) . 
Dipoles  with  a  length  greater  than  0.10X  can  no  longer  use  the 
assumption  of  triangular  current  distribution  with  acceptable 
accuracy.  The  finite  length  dipole  must  account  for  the 
sinusoidal  current  variation  along  the  length  of  the  antenna. 
The  sinusoidal  surface  current  distribution  for  a  finite 
length  dipole  oriented  as  depicted  in  Figure  4.1  is 

asIa  sin[P(h-zO], 

Kx'.y'.z')  =  Hz')  =  ■  (4.23) 

a,I0  sin[P(h+z/)],  -^zz'zO 

where  l0  is  the  peak  current  input  and  h 

Considering  the  finite  length  dipole  as  an  aggregate  of 
elemental  dipoles,  we  get  the  expression  for  the  electric  far- 
field  from  equation  4.7b  as  [Ref  2:  pp.  118-120] 

£,0=j‘Iid£’0=j-5l|sin(e)r|^J(x/,y/,z/)  e’^dz'  (4.24) 

1  '■ 

For  purposes  of  the  radiated  electric  field's  phase 
angle,  the  radial  field  distance  (R)  cannot  be  considered 
constant  while  integrating  over  the  length  of  the  dipole.  The 
variable  r  (which  stands  for  } R- R ' | )  in  the  exponential 
(phase)  term  is  approximated  as  [R-z* cos (9) ] ,  where  R  is  the 
distance  from  the  origin  to  the  field  point,  and  z'  is  the 
source  coordinate  of  the  vertical  dipole  along  the  z  axis 


[Ref  2:  pp.  112-114).  The 
approximation  is  depicted 
in  Figure  4.3,  and  is 
accurate  in  the  far- field 
for  r  >  21:A-  For 


Far  Field  Approximation  purposes  of  calculating 
r=R-z'cos(0) 

- electric  field  magnitude, 

FIGURE  4.3:  Phase  approximation  the  variation  of  r  over 
for  distance  variable,  r. 

the  length  of  the  dipole 
is  negligible,  and  the  distance  from  the  origin  to  the  field 
point  (spherical  coordinate  R)  is  used  in  the  denominator  of 
equation  4.24.  Substituting  equation  4.23  for  I(x',y',z') 
into  equation  4.24  and  (R-z'cos(0)]  for  r  in  the  exponential 
term  of  equation  4.24  and  integrating,  we  get 


1oPe 


sin (0) °  J0sin(P  (h+zO  ]e*:’Pz,C0S(e) dz'  (4.25) 

+  r*J0sin[P  (h-z')]  e*^*'cos(e)  dz ' 

J  0 


iloPe' 


■sin (0) 


|J^I0sin{P  (h+z')  ][cos  [pz'cos  (0) )  +jsin  [Pz'cos  (0)  ]  ]dz' 
+  f  AJ0sin[P  [h-z')  ][cos  [Pz'cos  (0)  ]  +js in  [Pz'cos  (0)  ]  )dz' 

J  0 


Since  the  current  distribution  is  an  even  function,  the 
odd  term  of  j  ■  sin (/3 •  z '  •  cos  ( 6) )  integrates  to  zero  over  the 
interval  from  -h  to  h.  Since  the  cos  03 -z'- cos  (0))  term  is 
even,  the  integral  can  be  written  as 
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n  Rj  p  *  ft 

E6=j  — - — sin.  (0)  |  sin[P  (h-z1)  lcos  [Pz'cos  (0)  ]  dz'  (4.26) 

2nR  Jo 

After  integrating,  the  final  expressions  for  radiated  electric 
and  magnetic  .  field  distributions  of  a  finite  vertical  dipole 
are  given  by  [Ref  2:  pp.  118-120] 


E*  -j 


Tlo-*oe"j(>*r  cos[p,hcos  (8)  ]  -  cos(Ph) 


2  nR 


sin (0) 


(4.27) 


E ’  j 

'  cosfP^cos  (8)  ]-cosPh 

t]  o  2nR 

sin(6) 

(4.28) 


The  coefficient  outside  tlje  brackets  is  called  the  element 
factor,  and  the  expression  inside  the  brackets  is  the  space 
factor. 

The  electric  field  radiation  pattern  is  a  function  of 
elevation  (0)  and  wavelength  (X)  .  For  a  fixed  wavelength,  the 

i 

radiation  pattern  for  any  vertical  plane  through  the  antenna 
axis  is  obtained  by  plotting  the  magnitude  of  equation  4.27  in 

r 

polar  coordinates  as  0  varies  from  zero  to  27t.  The  pattern  is 
normalized  by  multiplying  the  electric  field  magnitude  by  a 
factor  such  that  the  maximum  result  is  one.  Since  electric 
field  distribution  is  independent  of  4>,  the  pattern  for  -tts0sO 
is  symmetrical  with  that  for  Os0stt.  Figure  4.4  depicts  the 
normalized  radiation  patterns  of  four  finite  dipoles  with 
successive  lengths  of  X/4,  X/2,  3X/4,  and  X.  These  lengths 
simplify  calculation  of  the  electric  field  patterns  since  /?h 
is  7r/4,  7t/2,  37t/4,  and  ir  for  each  successive  dipole.  If 

dipole  length  is  greater  than  one  wave-length,  the  patterns 
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The  Mathcad  applications  evaluate  equation  4.30  to  obtain 
total  radiated  power  for  a  sinusoidal  input  current  with  a 
maximum  of  unity.  An  expression  for  radiation  resistance  is 
also  derived  by  Balanis  [Ref  2:  p.  124]  .  However,  the  Mathcad 
applications  calculate  radiation  resistance  from  total  average 
radiated  power  and  equation  2.10. 

When  the  dipole  is  elevated  above  the  origin  along  the  z 
axis,  an  additional  phase  shift  term  of  exp  [ -  j /UHqCos  ( 6 )  ]  must 


be  included  in  equations  4.27  and  4.28,  where  H0  is  the 
elevation  of  the  feed  at  the  center  of  the  dipole;  hence, 


..  lo-^o®  ffpcos  cosfphcos  (8)  ]  -  cos  (fth) 
2  rcR  sin(0) 
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(4.31) 
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(4.32) 


D.  GROUND  PLANE  EFFECTS  ON  THE  RADIATED  SPACE  WAVE 

The  electric  field  equations  derived  in  Chapters  3A-3C 
are  free  space  equations  which  do  not  consider  the  presence  of 
a  ground  plane.  Since  antenna  orientation  is  described  with 
respect  to  the  ground  plane,  orientation  descriptions  are  not 
really  valid  until  a  ground  plane  is  introduced.  The  free 
space  electric  field  distribution  for  a  dipole  oriented  along 
the  x  axis  or  y  axis  is  identical  to  the  distribution  for  the 
same  dipole  along  the  z  axis.  The  radiation  pattern  is 
displaced  to  correspond  to  the  new  position  of  the  dipole,  but 
the  radiated  electric  field  distribution  relative  to  the 
antenna  is  unaffected. 

The  free  space  electric  field  equations  no  longer  provide 
an  accurate  description  of  an  antenna's  radiated  electric 
field  distribution  above  a  ground  plane.  The  free  space 
equations  only  provide  the  direct  wave  component  of  the  space 
wave  which  propagates  directly  from  the  antenna  to  the  field 
points  above  ground  (-ir/2  <  6  <  n/2)  .  They  are  not  valid  for 
|0|  >  tt/2  (points  below  the  ground  plane). 


A  complete  description  of  the  space  wave  above  a  ground 
plane  must  include  the  reflected  wave  component.  The 
reflected  component  consists  of  electromagnetic  waves  radiated 
at  elevation  angles  corresponding  to  |0|  >  n/2  in  the  free 
space  equations.  These  waves  strike  the  ground  plane  and  are 
absorbed,  partially  reflected,  or  completely  reflected  before 
reaching  a  field  point  above  the  ground  plane.  The  total 
space  wave  is  obtained  by  adding  the  direct  wave  and  reflected 
wave  components  for  field  points  defined  by  -n/2  <  6  <  ir/2. 
The  distribution  of  the  reflected  waves  vary  relative  to  the 
antenna  as  its  spatial  orientation  changes,  so  the  electric 
field  distribution  (radiation  pattern)  is  also  affected  when 
the  antenna  is  repositioned  with  respect  to  the  ground. 

The  reflected  wave  component  is  determined  using 
reflection  coefficients  and  image  theory.  Reflection 
coefficients  determine  how  much  of  an  incident  electric  field 
is  reflected  by  a  given  ground  plane,  and  they  are  a  function 
of  6  and  the  complex  index  of  refraction.  The  complex  index 
of  refraction  (n)  in  terms  of  frequency,  conductivity  (a),  and 
relative  dielectric  constant  (sr)  of  the  ground  plane  is 
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The  vertical  reflection  coefficient  (rv) 


T  _  •n2c°s  (9r>  -yfa2- sin2  (9f) 
n2cos  (8r)  +^n2- sin2  (©TT 


(4.34) 
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is  used  to  find  the  reflected  wave  contribution  to  the  0 
component  of  the  space  wave.  The  horizontal  reflection 
coefficient  (rh) 


r  cos  (0f)  -Jn2- sin2  (0r) 
1  n~. 


(4.35) 


cos  (0r)  +^n2- sin2  (0r) 

is  used  to  find  the  reflected  wave  contribution  to  the  <f> 
component  of  the  space  wave  [Ref  5:  pp.  629-633]. 

Antenna  image  theory  accounts  for  phase  differences 
between  the  direct  and  reflected  wave  components.  Phase 
differences  are  caused  by  the  different  distances  travelled  by 
the  two  components  to  a  given  field  point.  The  image  is  a 
'mirror  image'  located  below  the  ground  plane  at  a  distance 
equal  to  the  antenna  height  above  ground.  The  distance  from 
the  image  to  a  given  field  point  is  the  distance  travelled  by 
the  reflected  wave  to  that  field  point.  This  distance  is 
approximated  by  R+HoCos(0),  similar  in  derivation  to  the 
approximation  illustrated  in  Figure  4.4  for  the  direct  wave 
distance  of  R-Hocos(0).  The  distances  travelled  by  the  two 
components  differ  by  2HoCos(0),  but  the  direct  wave  expression 
already  includes  the  term  exp [+j/?HoCos  (0)  ]  ,  so  the  reflected 
wave  expression  includes  the  term  exp  [- j/3H0cos  (0)  ]  l:o  account 
for  any  resulting  phase  differences.  \ 

In  equations  4.34  and  4.35,  the  variable  0r  is  the  angle 
between  the  z  axis  and  the  radial  from  antenna  image  to  field 
point.  If  the  height  of  the  antenna  is  small  in  comparison  to 
the  far- field  distance,  0r  is  approximated  as  0  (the  angle 
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between  the  z  axis  and  the  radial  from  origin  co  field  point)  . 
The  Mathcad  applications  use  this  approxima'-  ‘.on  to  find  the 
reflection  coefficients  for  632  discrete  values  of  8  (316  for 
the  array  applications)  in  the  interval  -tt/2  <  8  <  7r/2  at  a 
specified  distance  (R)  from  the  antenna. 

The  direct  wave  component  of  the  space  wave  is  calculated 
for  discrete  values  of  8  with  the  free  space  electric  field 
equation  and  the  phase  terms  derived  from  the  antenna  height. 
The  reflected  wave  components  are  then  calculated  for  each  8 
|by  multiplying  the  direct  wave  components  by  their 

i  '  '  . 

i corresponding  reflection  coefficients  and  applicable  phase 

Iterms.  The  reflected  wave  components  are  then  added  to  the 

Idirect  wave  components  to  obtain  the  total  radiated  space  wave 

at  each  8 .  The  discrete  results  for  the  space  wave  are 

jnormalized  with  respect  to  the  value  with  maximum  magnitude, 

and  the  space  wave  radiation  pattern  is  depicted  by  plotting 

the  normalized  field  magnitudes  for  -tt/2  <  8  <  n/2. 

The  Mathcad  applications  plot  radiation  patterns  assuming 
a  planar  surface  below  the  antenna.  The  planar  earth  model 
provides  accurate  results  for  (1)  antenna  heights  within 
several  wavelengths  of  the  surface  and  (2)  wavelengths  much 
smaller  than  the  radius  of  the  earth.  The  antennas  addressed 
in  this  report  are  usually  mounted  well  within  several 
wavelengths  of  the  ground,  so  the  planar  assumption  in  not  a 
problem  for  the  majority  of  calculations.  If  a  configuration 
is  encountered  where  the  antenna  is  more  than  four  or  five 
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wavelengths  above  the  ground,  the  Mathcad  applications  still 
provide  a  fair  estimate  of  the  radiation  parameters.  However, 
if  highly  accurate  resu : cs  are  required,  a  computer  program 
which  accounts  for  the  curvature  of  the  earth  should  be  used. 

The  final  expression  for  the  radiated  space  wave  of  a 
vertical  dipole  antenna  is  [Ref  6:  p.  166] 
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This  is  the  equation  used  in  the  Elevated  Vertical  Dipole 
application  to  calculate  the  space  wave  radiation  pattern  and 
radiated  electric  field  distribution. 


E.  THE  SURFACE  WAVE 

The  radiated  electric  field  of  an  antenna  above  a  ground 
plane  includes  a  surface  wave  component  in  addition  to  the 
space  wave.  The  surface  wave  radiation  pattern  is  a  function 
of  the  distanc?  from  the  antenna  and  the  relative  dielectric 
constant  and  conductivity  of  the  ground  plane.  Since  the 
surf act  wave  attenuates  exponentially  with  increasing  range, 
the  surface  wave  radiation  pattern  varies  with  distance  from 
the  antenna.  While  the  space  wave  can  be  calculated  at 


different  distances  from  the  antenna,  its  radiation  pattern  is 
predominately  invariant  with  distance  from  the  antenna. 


Equations  for  the  radiated  electric  fields  of  a  vertical 
antenna  above  planar  earth  were  derived  by  Sommerfeld  in  1909 
[Ref  7] ,  but  the  final  form  of  these  equations  is  complicated 
and  difficult  to  evaluate.  Norton  simplified  these  equations 
to  a  form  suitable  for  engineering  applications  [Ref  8]  . 
Norton's  expression  for  the  surface  wave  of  an  elemental 
vertical  dipole  above  a  ground  plane  is  [Ref  5:  p.  644] 


(  Ocos(8)]\ 

j30$Idl  (l-rv)  F  I-? - - - I* 
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(4.37) 
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Its  derivation  is  not  explained  here,  but  references  6  and  7 
give  a  detailed  derivation.  In  equation  4.37,  R  is  distance 
to  the  field  point,  £  and  £  are  unit  vectors  perpendicular  and 
parallel  to  the  antenna  axis  respectively,  u  is  the  inverse  of 
the  complex  index  of  refraction  (n)  ,  and  F  is  the  surface  wave 
attenuation  factor.  The  function  F  introduces  a  surface  wave 
attenuation  dependent  upon  distance,  frequency,  and  the 
properties  of  the  earth  over  which  the  wave  is  travelling. 
The  attenuation  factor  is  given  by  [Ref  5:  pp.  644-651] 
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where  Pe  is  the  complex  numerical  distance  for  vertical 
polarization  and  the  erfc  term  is  the  complementary  error 
function  defined  by 

ezfc(jfPl)  =—  f  e'*2dx  (4.41) 


Mathcad  is  incapable  of  evaluating  the  complementary  error 
function  for  complex  numbers,  so  an  asymptotic  expansion  and 
infinite  series  approximation  [Ref  9]  was  written  into  the 
Mathcad  code  to  evaluate  the  erfc  terms. 

When  the  magnitude  of  j  •v/p7  is  less  than  2.18,  the 
Mathcad  applications  use 


erf(x>jy) =erf{x)  + 
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where 


fn{x,y)  =2x-2x  cosh (ny) cos (2xy)  +n  sinh (ny) sin (2xy) 
gn{x,y)  =2x  cosh (ny)  sin(2xy)  +n  sinh  (ny)  cos  (2xy) 
|e(x,y)  |  =  10"1S  | erf  (x+jy)  | 


an  infinite  series  approximation  which  evaluates  erfCjV'p]’), 
where  ’erf'  is  the  standard  error  function.  [Ref  9  p.  299]. 
The  solution  to  erfc(jV^7)  is  then  found  by 

erfc  (x+jy)  =1- erf  (x+jy)  (4.43) 
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When  the  magnitude  of  j Vp7  is  greater  than  2.18,  the 


Mathcad  code  uses 
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an  asymptotic  expansion  which  solves  for  j  VjrP,,  *erfc  ( jv/p7) 
[Ref  9  p.  298] .  Either  of  these  solutions  can  be  substituted 
directly  into  equation  4.38  to  find  the  surface  wave 
attenuation  factor. 

A  magnitude  of  2.18  is  used  as  the  transition  point  in  the 
Mathcad  applications  because  it  provides  the  smoothest 
transition  in  both  magnitude  and  phase  between  the  series 
approximation  and  asymptotic  expansion  results. 

For  antennas  which  radiate  an  electric  field  with  a  <b 
component,  the  surface  wave  attenuation  factor  for  horizontal 
polarization  (Fm)  must  be  found  using 


and  Pm * 
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where  PB  is  the  complex  numerical  distance  for  horizontal 
polarization.  The  variable  0,  in  equations  4.40  and  4.46  for 


P,  and  Pm,  respectively  'is  approximated  as  6  (just  as  it  was 
for  the  reflection  coefficients  in  equations  4.34  and  4.35). 
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The  numerical  distances  (Pe  and  Pra)  and  the  surface  wave 
attenuation  factors  (Fe  and  Fm)  are  evaluated  for  O°s0s9Q°  to 
determine  surface  wave  radiation  patterns.  Their  value  is  of 
most  interest,  however,  directly  at  the  surface  (0=90°)  where 
the  magnitude  of  F  is  called  the  ground  wave  attenuation 
factor.  Numerical  distance  is  proportional  to  distance  from 
the  antenna  and  the  square  of  the  frequency,  while  it  is 
approximately  inversely  proportional  to  ground  conductivity. 
When  the  magnitude  of  the  numerical  distance  is  less  than 
about  4.5,  empirical  formulas  show  that  the  ground  wave 
attenuation  factor  varies  exponentially  with  numerical 
distance.  For  numerical  distance  magnitudes  greater  than  4.5, 
empirical  formulas  show  that  the  ground  wave  attenuation 
factor  is  inversely  proportional  to  the  numerical  distance. 
Numerical  distance  magnitudes  are  typically  much  greater  than 
4.5  in  the  far- field.  This  implies  that  the  ground  wave 
(0=90°)  field  strength  varies  inversely  with  the  square  of  the 
distance  from  the  antenna  since  numerical  distance  is  directly 
proportional  to  distance,  from  the  antenna  and  equation  4.37 
shows  an  explicit  l/R  dependence  in  addition  to  the 
attenuation  factor. 

The  Mathcad  applications  determine  the  complex  numerical 
distances  for  the  same  discrete  values  of  0  described  in 
Chapter  4.D.  After  determining  the  complementary  error 
function  for  each  numerical  distance  result,  the  surface  wave 
attenuation  factor  is  calculated  for  each  discrete  0.  The 
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surface  wave  is  evaluated  for  each  0  using  the  equations  for 
the  electric  field  distribution,  and  the  discrete  results  of 
the  surface  wave  are  normalized  with  respect  to  the  value  with 
maximum  magnitude.  The  surface  wave  pattern  is  depicted  by 
plotting  the  normalized  magnitudes  for  -7r/2  <  9  <  tr/2.  The 
numerical  distances  (Pc  and  Pm)  at  the  surface  (0=90°)  are  also 
listed  with  the  predicted  radiation  parameters. 

With  all  of  the  terms  defined,  equation  4.37  is 
integrated  over  the  length  of  a  generic  finite  length  vertical 
dipole  to  obtain  the  radiated  surface  wave  distribution  of  an 
elevated  vertical  dipole  [Ref  6  p.166] 


Ee=j60I0 
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The  total  radiated  far- field  electric  field  distribution 
of  an  elevated  vertical  dipole  above  a  planar  earth  can  now  be 
expressed  by  combining  equations  4.36  and  4.37  to  write 


Eq  =j60I^  - 
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The  first  two  terms  inside  the  brackets  represent  the  space 
wave  and  the  third  term  represents  the  surface  wave. 

Chapter  IV  has  defined  the  necessary  terms  and  provided 
a  detailed  derivation  of  the  total  electric  field  distribution 
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for  an  elevated  vertical  dipole  antenna.  The  result  is  a 
radiated  electric  field  distribution  with  a  6  component  only. 
Most  antenna  types  radiate  an  electric  field  with  both  6  and 
<p  components.  While  a  detailed  derivation  of  an  electric 
field  <t>  component  was  not  presented  here,  the  necessary 
horizontal  polarization  terms  associated  with  deriving  <p 
components  were  discussed.  When  these  horizontal  terms  are 
substituted  for  their  vertical  counterparts,  the  method  used 
to  derive  the  electric  field  equation  for  the  4>  component  of 
a  radiated  space  wave  or  surface  wave  is  identical  to  that 
used  in  deriving  the  6  component  for  the  elevated  vertical 
dipole. 

The  remaining  chapters  of  this  thesis  explain  the  antenna 
configuration  and  the  calculations  of  each  individual  Mathcad 
application.  Only  the  final  expressions  for  the  radiated  far- 
field  electric  field  distributions  are  provided  since  the 
derivation  of  each  follows  the  same  procedure  as  explained  in 
this  chapter.  There  is  a  lot  of  text  repetition  among  the 
remaining  chapters  because  each  is  written  as  a  stand  alone 
users  guide  for  its  associated  Mathcad  application.  Any 
important  concepts  not  presented  in  previous  chapters  are 
discussed  as  required. 


The  orientation  of  the  elevated  vertical  dipole  is 
depicted  in  Figure  5.1,  where  1  is  the  length  of  the  dipole, 
Ho  is  the  height  of  the  feed  above  ground,  R  is  the  radial 
coordinate,  6  is  the  elevation  coordinate,  and  <t>  is  the 


azimuth  coordinate. 


dipole  antenna  for  its  corresponding  Mathcad  application. 


The  Mathcad  application  for  the  elevated  vertical  dipole 
requires  the  following  inputs: 


1 . length  of  the  dipole  (meters) 

fj . operational  frequency  (Hertz) 

H0 . antenna  feed  height  above  ground  plane  (meters) 

R . distance  from  antenna  (meters) 

er . relative  dielectric  constant  of  ground  plane 

a . conductivity  of  ground  plane 


The  length  of  the  dipole  is  used  to  determine  the  frequencies 
which  correspond  to  a  \/4,  A/2,  3\/4,  or  A  length  dipole.  The 
user  inputs  the  frequency  and  the  radiation  parameters 
discussed  in  Chapter  3  are  calculated.  Feed  height  ( Hq)  must 


be  at  least  equal  to  the  half-length  of  the  dipole  (h=l/2) , 
and  the  distance  from  the  antenna  (R)  must  meet  the  far- field 
requirements  of  Chapter  3.C. 

The  equation  for  the  radiated  electric  field  distribution 
of  the  elevated  vertical  dipole  was  derived  in  Chapter  4  as 
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The  first  two  terms  inside  the  brackets  represent  the  space 
wave  and  the  third  term  represents  the  surface  wave.  I0  is 
one  since  a  sinusoidal  current  input  with  a  maximum  of  unity 
is  assumed. 

The  requested  inputs  are  used  to  calculate  variables  at 
the  top  of  page  41  for  632  discrete  values  of  6  which  are 
equally  incremented  from  -tt/2  to  ir/2. 


h . half-length  of  the  dipole  (meters) 

Xj  . wavelength  of  the  operational  frequency  (meters) 

0 . free  space  wavenumber  for  operational  frequency 

Rd . distance  from  antenna  to  field  point  (meters) 

Rr . .  .distance  from  antenna  image  to  field  point  (meters) 

n . index  of  refraction 

Pc . complex  numerical  distance  (vertical  polarization) 

rv . vertical  reflection  coefficient 

Fe  . vertical  surface  wave  attenuation  factor 


The  calculated  variables  are  used  to  evaluate  the 
radiated  far- field  space  wave  and  surface  wave  for  each 
discrete  0.  The  space  wave  and  surface  wave  results  are 
combined  for  corresponding  values  of  0  to  obtain  the  total 
radiated  electric  field  distribution.  The  space  wave,  surface 
wave,  and  total  electric  field  results  are  then  normalized 
with  respect  to  the  maximum  field  intensity  of  each,  and  the 
normalized  magnitudes  are  plotted  for  each  value  of  0  to 
depict  the  radiation  patterns.  The  patterns  are  valid  for  any 
vertical  plane  containing  the  antenna  axis,  because  the 
radiated  electric  field  does  not  vary  with  changes  in  <t>. 

Equations  3.8  and  3.9  are  used  to  integrate  equation  5.1 

over  the  hemispherical  Gaussian  surface  above  the  ground  plane 

at  a  fixed  radius  (R)  from  the  antenna  to  find  total  average 

radiated  power  (Ptld)  .  With  the  discrete  values  of  the  electric 

field  and  total  average  radiated  power  determined,  the  Mathcad 

application  predicts  the  following  radiation  characteristics 

from  the  equations  in  Chapter  3: 

Rrad  . radiation  resistance  (Ohms) 

D0  . directivity 

EIRP . effective  radiated  isotropic  power  (Watts) 

An,,,  . maximum  theoretical  effective  area  (square  meters) 

lna,  . maximum  theoretical  effective  length  (meters) 

P . numerical  distance  (vertical  polarization,  0=90°) 

AnSflemw  .elevation  angle  of  maximum  directive  gain  (degrees) 
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These  results  are  consistent  with  expectations  for  this 
particular  configuration.  Appendix  A  contains  computer 
hardcopies  of  additional  example  calculations  for  the  vertical 
dipole  and  compares  predicted  radiation  parameters  to  those 
expected  based  on  previous  calculations. 
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THE  VERTICAL  MONOPOLE 


The  orientation  of  the  vertical  monopole  is  depicted  in 
Figure  6.1,  where  h  is  the  length  of  the  monopole,  the  antenna 
feed  is  at  ground  level,  R  is  the  radial  coordinate,  6  is  the 
elevation  coordinate,  and  <p  is  the  azimuth  coordinate. 


i 

z  axis 

•  FIELD  POINT 
( R,0,<p ) 

h. 

R 

/ 

/ 

y  axis 

$ 

x  axis 

Vertical  Monopole 

Orientation 

FIGURE  6.1j  Spatial  orientation  of  the  vertical  monopole 
antenna  for  its  corresponding  Mathcad  application. 


The  Mathcad  application  for  the  vertical  monopole 

requires  the  following  inputs: 

h . length  of  the  monopole  (meters) 

f5 . operational  frequency  (Hertz) 

R . distance  from  antenna  (meters) 

er . relative  dielectric  constant  of  ground  plane 

a. ...... .conductivity  of  ground  plane 

The  length  of  the  monopole  is  used  to  determine  frequencies 

which  correspond  to  a  X/8(  X/4,  3X/8,  or  X/2  long  monopole. 

The  user  inputs  the  frequency  for  which  the  radiation 

parameters  discussed  in  Chapter  3  are  calculated.  Distance 

from  the  antenna  (R)  must  meet  the  far- field  requirements  of 

Chapter  3 . C . 

The  monopole's  electric  field  distribution  is  derived  in 
the  same  way  as  the  vertical  dipole  equation  was  derived  in 
Chapter  4  and  is  given  by  [Ref  6:  p.152] 
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fl=sin[ph  cos (0)  ]  -cos (0)  sin(Ph) 


The  first  two  terms  inside  the  brackets  represent  the  space 
wave  and  the  second  two  terms  represent  the  surface  wave.  The 
vertical  monopole's  far- field  electric  field  has  a  6  component 
only  ,  just  as  the  vertical  dipole.  I0  is  one  since  a 
sinusoidal  current  input  with  a  maximum  of  unity  is  assumed. 
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The  requested  inputs  are  used  to  calculate  the  following 
variables  for  632  discrete  values  of  6  which  are  equally 
incremented  from  -tt/2  to  tt/2: 


X5 . wavelength  of  the  operational  frequency  (meters) 

18 . free  space  wavenumber  for  operational  frequency 

n. .......  index  of  refraction 

Pe . complex  numerical  distance  (vertical  polarization) 

rv . vertical  reflection  coefficient 

Fe  . vertical  surface  wave  attenuation  factor 


The  calculated  variables  are  used  to  evaluate  the 
radiated  far- field  space  wave  and  surface  wave  for  each 
discrete  6.  The  space  wave  and  surface  wave  results  are 
combined  for  corresponding  values  of  6  to  obtain  the  total 
radiated  electric  field  distribution.  The  space  wave,  surface 
wave,  and  total  electric  field  results  are  then  normalized 
with  respect  to  the  maximum  field  intensity  of  each,  and  the 
normalized  magnitudes  are  plotted  for  each  value  of  0  to 
depict  the  radiation  patterns.  The  patterns  are  valid  for  any 
vertical  plane  containing  the  antenna  axis,  because  the 
radiated  electric  field  does  not  vary  with  changes  in  <p. 

Equations  3.8  and  3.9  are  used  to  integrate  equation  6.1 
over  the  hemispherical  Gaussian  surface  above  the  ground  plane 
at  a  fixed  radius  (R)  from  the  antenna  to  find  the  total 
average  radiated  power  (Prid)  .  With  the  discrete  values  of  the 
electric  field  and  total  average  radiated  power  determined, 
the  Mathcad  application  predicts  the  radiation  characteristics 
listed  at  the  top  of  the  next  page  from  the  equations  in 
Chapter  3 . 
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Radiation  characteristics  determined  by  the  vertical 

monopole  application  are: 

Rra<i  . radiation  resistance  (Ohms) 

■  D0  . directivity 

EIRP . effective  radiated  isotropic  power  (Watts) 

A^  . maximum  theoretical  effective  area  (square  meters) 

l^  . maximum  theoretical  effective  length  (meters) 

Pe  .  numerical  distance  (vertical  polarization,  5=90°) 

Angle^  .  -  levation  angle  of  maximum  directive  gain  (degrees) 

As  an  example,  the  Mathcad  monopole  application  was 
executed  with  the  following  inputs: 


monopole  length 
frequency 

distance  from  the  antenna 
relative  dielectric  constant 
ground  conductivity 


0.5  meters 
150 -106  Hertz 
3000  meters 
15 

5  •  10'3 


Figures  6.2  and  6.3  depict  the  radiation  patterns  for  the 
space  wave  and  surface  wave  for  this  example. 


FIGURE  6.2s  Monopole  space  FIGURE  6.3:  Monopole  surface 
wave  radiation  pattern.  wave  radiation  pattern. 

The  following  radiation  parameters  were  predicted  by  the 
Mathcad  application  for  a  sinusoidal  current  input  of  one  Amp: 


Total  power  radiated 
Radiation  resistance 
Directivity 

Effective  isotropic  radiated  power 
Maximum  theoretical  effective  area 
Maximum  theoretical  effective  length 
Numerical  distance  (vertical) 
Elevation  of  directivity 


5.607  Watts 
11.215  Ohms 
3.246 

18.203  Watts 
1.033  sq  meters 
0.351  meters 
293.016 

26.815  degrees 


These  results  are  consistent  with  expectations  for  this 


particular  configuration.  Appendix  B  contains  computer 
hardcopies  of  additional  example  calculations  for  the  vertical 
monopole  and  compares  predicted  radiation  parameters  to  those 
expected  based  on  previous  calculations. 


The  orientation  of  the  horizontal  dipole  is  depicted  in 
gure  7.:,  where  1  is  length  of  the  dipole,  Hd  is  the  height 
the  feed  above  ground,  R  is  the  radial  coordinate,  0  is  the 
evation  coordinate,  and  4>  is  the  azimuth  coordinate. 


antenna  for  its  corresponding  Mathcad  application. 
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The  Mathcad  application  for  the  horizontal  dipole  requires 
the  following  inputs: 


1 . length  of  the  dipole  (meters) 

H0 . antenna  feed  height  above  ground  plane  (meters) 

f5 . operational  frequency  (Hertz) 

R . distance  from  antenna  (meters) 

e, . relative  dielectric  constant  of  ground  plane 

cr . conductivity  of  ground  plane 

d . elevation  angle  index  (from  Table  3.1) 


The  length  of  the  dipole  is  used  to  determine  frequencies 
which  correspond  to  a  X/4,  X/2,  3X/4,  or  X  length  dipole.  The 
user  inputs  the  frequency  for  which  the  radiation  parameters 
discussed  in  Chapter  3  are  computed.  Feed  height,  H0,  must  be 
greater  than  or  equal  to  zero,  and  distance  from  the  antenna 
(R)  must  meet  the  far- field  requirements  of  Chapter  3.C.  The 
elevation  angle  index  (d)  sets  the  elevation  coordinate,  6, 
for  which  a  horizontal  radiation  pattern  is  determined.  Table 
3.1  lists  possible  indices  and  their  corresponding  elevation 
angles  from  .285°  to.  88. 857°  in  increments  of  about  2°. 
Indices  between  those  listed  can  be  used  to  interpolate  a 
better  approximation  of  a  desired  elevation. 

The  horizontal  dipole's  radiated  electric  field  is  a 
combination  of  d  and  $  components.  In  addition,  unlike  the 
vertical  dipole,  its  radiation  pattern  varies  with  changes  in 
<t>.  The  total  radiated  electric  field  is  the  vector  sum  of  the 
d  and  $  components.  The  electric  field  for  the  horizontal 
dipole  is  obtained  in  a  manner  analogous  to  that  used  for  the 
vertical  dipole  discussed  in  Chapter  4.  The  horizontal 
dipole's  d  and  $  field  components  are  given  by  [Ref  6:  p.  144] 


V 
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The  first  two  terms  inside  the  brackets  of  each  equation  are 
the  space  wave  components  and  the  third  terms  are -the  surface 
wave  components.  I0  is  taken  to  be  one,  since  a  sinusoidal 
current  input  with  a  maximum  of  unity  is  assumed. 

The  requested  inputs  are  used  to  calculate  the  following 
variables  using  a  constant  <f>  of  0=0  and  0=7r/2  for  632  discrete 
values  of  6  which  are  equally  incremented  from  -7r/2  to  tt/2: 


h . half-length  of  the  dipole  (mecers) 

X3 . wavelength  of  the  operational  frequency  (meters) 

/3 . free  space  wavenumber  for  operational  frequency 

Rd . distance  from  antenna  to  field  point  (meters) 

Rr . distance  from  antenna  image  to  field  point  (meters) 

n . index  of  refraction 

. complex  numerical  distance  (vertical  polarization) 

Pm  . complex  numerical  distance  (horizontal  polarization) 

rv . vertical  reflection  coefficient 

rh . horizontal  reflection  coefficient 

Fe . vertical  surface  wave  attenuation  factor 

Fm  . horizontal  surface  wave  attenuation  factor 


The  calculated  variables  are  used  to  evaluate  the  far- 
field  space  wave  and  surface  wave  for  the  discrete  values  of 
6  for  0=0  and  0= ir/2 .  The  total  space  wave  and  surface  wave 
distributions  are  determined  from  vector  addition  of 
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corresponding  d  and  $  components.  The  space  wave  and  surface 
wave  results  are  then  combined  for  corresponding  discrete 
values  of  6  to  obtain  the  total  radiated  electric  field 
distribution  for  the  0=0  and  0=7r/2  vertical  planes.  The  space 
wave,  surface  wave,  and  total  electric  field  results  are  then 
normalized  with  respect  to  the  maximum  field  intensity  of 
each,  and  the  normalized  magnitudes  are  plotted  for  each 
discrete  d  to  depict  the  radiation  patterns.  The  Mathcad 
horizontal  dipole  application  computes  the  space  wave,  surface 
wave,  and  total  electric  field  radiation  patterns  and 

radiation  parameters  in  the  0=0  and  0=7r/2  vertical  planes. 

'  i 

The  variables  corresponding  to  the  selected  elevation 

I 

angle  index  (d)  are  used  to  evaluate  the  radiated  electric 
field  components  for  the  space  wave  and  surface  wave  at  the 

fixed  elevation  angle  (0d)  as  0  varies  from  0  to  2tt  in  632 

j  ■ 

equal  increments.  jThe  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  for  the  vertical  planes. 

To  find  the  contributions  of  the  #  and  $  electric  field 
components  to  the  total  average  radiated  power  (Pr,d)  ,  equations 
3.8  and  3.9  are  used  to  integrate  equations  7.1  and  7.2  over 
the  hemispherical  Gaussian  surface  above  the  ground  plane  at 
a  fixed  radius  (R)  from  the  antenna  With  the  discrete  values 
of  the  electric  field  and  total  average  radiated  power,  the 
Mathcad  application  predicts  the  following  radiation 
characteristics  from  the  equations  in  Chapter  3: 


52 


Rrid  . radiation  resistance  (Ohms) 

D0 . .directivity 

EIRP . effective  radiated  isotropic  power  (Watts) 

. maximum  theoretical  effective  area  (square  meters) 

lmjx  . maximum  theoretical  effective  length  (meters) 

P . numerical  distance  (vertical  polarization,  0=90°) 

Pm  . numerical  distance  (horizontal  polarization,  6  =  90°) 


Anglemjx  .elevation  angle  of  maximum  directive  gain  (degrees) 
The  directivity  (D0)  ,  effective  isotropic  radiated  power 
(EIRP) ,  maximum  theoretical  effective  area  (A^,)  ,  maximum 
theoretical  effective  length  (lm„) ,  and  the  elevation  angle  of 
maximum  directive  gain  (Angle^J  are  all  determined  for  both 


the  0=0  and  0=7r/2  vertical  planes. 

As  an  example,  the  Mathcad  horizontal  dipole  application 


was  executed  with  the  following  inputs: 


dipole  length 
frequency 

height  of  antenna  feed 
distance  from  the  antenna 
relative  dielectric  constant 
ground  conductivity 
elevation  angle  index 


1.0  meter 
150  *106  Hertz 
1.0  meter 
3000  meters 
15 
10*2 

535  (-27°) 


The  following  radiation  parameters  were  predicted  by  the 


Mathcad  application  for  a  sinusoidal  current  input  of  one  Amp: 


Total  power  radiated  (Watts) 

26.116 

Radiation  resistance  (Ohms) 

52.233 

Numerical  distance  (vertical) 

292.419 

Numerical  distance  (horizontal) 

66215.4 

0=0 

6- tt/2 

Directivity 

1.310 

7.215 

EIRP  (Watts) 

34.221 

188.435 

Max  eff  area  (sq  meters) 

0.417 

2.297 

Max  eff  length  (meters) 

0.481 

1.128 

Angle™*  (degrees) 

45.071 

28.241 

Figures  7.2  through  7.7  are  the  space  wave  and  surface 


wave  radiation  patterns  in  the  0=0  and  0=ir/2  vertical  planes 


and  the  designated  horizontal  plane  for  this  example. 
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FIGURE  7.6:  Horizontal  FIGURE  7.7:  Horizontal 
dipole  space  wave  radiation  dipole  surface  wave  radiation 
pattern  for  horizontal  plane.  pattern  for  horizontal  plane. 

These  results  are  consistent  with  expectations  for  this 
particular  configuration.  Appendix  C  contains  computer 
hardcopies  of  additional  example  calculations  for  the 
horizontal  dipole  and  compares  predicted  radiation  parameters 
to  those  expected  based  on  previous  calculations. 


The  Mathcad  application  for  the  arbitrarily  oriented 
dipole  requires  the  following  inputs: 


1.. . length  of  the  dipole  (meters) 

H0 . antenna  feed  height  above  ground  plane  (meters) 

a . .  .angle  of  antenna  axis  with  the  horizontal  (degrees) 

fj . operational  frequency  (Hertz) 

R . distance  from  antenna  (meters) 

er . relative  dielectric  constant  of  ground  plane 

c . .conductivity  of  ground  plane 

d . elevation  angle  index  (from  Table  3.1) 


The  length  of  the  dipole  is  used  to  determine  frequencies 
which  correspond  to  a  A/4,  A/2,  3 A/4,  or  A  length  dipole.  The 
user  inputs  the  frequency  for  which  the  radiation  parameters 
discussed  in  Chapter  3  are  computed.  Feed  height,  H0,  must  be 
greater  than  or  equal  to  (1/2) -sin(a) ,  and  the  distance  from 
the  antenna  (R)  must  meet  the  far- field  requirements  of 
Chapter  3.C.  The  elevation  angle  index  (d)  sets  the  elevation 
coordinate,  6,  for  which  a  horizontal  radiation  pattern  is 
determined.  Table  3.1  lists  possible  indices  and  their 
corresponding  elevation  angles  from  .285°  to  88.857°  in 
increments  of  about  2°.  Indices  between  those  listed  can  be 
used  to  interpolate  a  better  approximation  of  a  desired 
elevation. 

The  radiated  electric  field  for  the  arbitrarily  oriented 
dipole  is  a  combination  of  d  and  $  components,  and  its 
radiation  pattern  varies  with  changes  in  4>.  The  total 
electric  field  is  the  vector  sum  of  the  §  and  $  components. 
The  arbitrarily  oriented  dipole's  radiated  electric  field  is 
obtained  in  a  manner  analogous  to  that  used  for  the  vertical 
dipole  discussed  in  Chapter  4,  and  is  given  by 
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Equation  8.1  is  the  §  component  (vertically  polarized),  and 
equation  8.2  is  the  $  (horizontally  polarized)  component  [Ref 
6:  p.  174-175].  The  first  two  terms  inside  the  brackets  of 
each  equation  represent  the  space  wave  components  and  the 
third  terms  represent  the  surface  wave  components.  I0  is 
taken  to  be  one,  since  a  sinusoidal  current  input  with  a 
maximum  of  unity  is  assumed. 
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The  requested  inputs  are  used  to  calculate  the  following 
variables  using  a  constant  0  of  0=0  and  4>=tt/2  for  632  discrete 
values  of  6  which  are  equally  incremented  from  -ir/2  to  tr/2: 


h . half-length  of  the  dipole  (meters) 

X5 . wavelength  of  the  operational  frequency  (meters) 

0 . free  space  wavenumber  for  operational  frequency 

Rd . distance  from  antenna  to  field  point  (meters) 

Rr . distance  from  antenna  image  to  field  point  (meters) 

n . .index  of  refraction 

Pe . complex  numerical  distance  (vertical  polarization) 

Pm  . complex  numerical  distance  (horizontal  polarization) 

Tv . vertical  reflection  coefficient 

rh . horizontal  reflection  coefficient 

Fe . vertical  surface  wave  attenuation  factor 

Fm  . horizontal  surface  wave  attenuation  factor 


The  calculated  vari ah1  as  are  used  to  evaluate  the  far- 
field  space  wave  and  surface  wave  for  the  discrete  values  of 
6  for  0=0  and  0=7r/2 .  The  total  space  wave  and  surface  wave 
distributions  are  determined  from  vector  addition  of 
corresponding  d  and  $  components.  The  space  wave  and  surface 
wave  results  are  then  combined  for  corresponding  discrete 
values  of  6  to  obtain  the  total  radiated  electric  field 
distribution  for  the  0=0  and  0=ir/2  vertical  planes.  The  space 
wave,  surface  wave,  and  total  electric  field  results  are  then 
normalized  with  respect  to  the  maximum  field  intensity  of 
each,  and  the  normalized  magnitudes  are  plotted  for  each 
discrete  6  to  depict  the  radiation  patterns.  The  arbitrarily 
oriented  dipole  Mathcad  application  computes  the  space  wave, 
surface  wave,  and  total  electric  field  radiation  patterns  and 
radiation  parameters  in  the  0=0  and  0=7r/2  vertical  planes. 

The  variables  corresponding  to  the  selected  elevation 
angle  index  (d)  are  used  to  evaluate  the  radiated  electric 


59 


field  components  for  the  space  wave  and  surface  wave  at  the 
fixed  elevation  angle  (0d)  as  0  varies  from  0  to  2 7r  in  632 
equal  increments.  The  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  in  the  vertical  planes. 

To  find  the  contributions  of  the  §  and  $  electric  field 

components  to  the  total  average  radiated  power,  Prad,  equations 

3.8  and  3.9  are  used  to  integrate  equations  8.1  and  8.2  over 

the  hemispherical  Gaussian  surface  above  the  ground  plane  at 

a  fixed  radius  (R)  from  the  antenna.  The  sum  of  the  two 

integrals  is  the  antenna's  Prad.  With  the  discrete  values  of 

the  electric  field  and  total  average  radiated  power 

determined,  the  Mathcad  application  predicts  the  following 

radiation  characteristics  from  the  equations  in  Chapter  3: 

Rrad  . radiation  resistance  (Ohms) 

D0 . directivity 

EIRP . effective  radiated  isotropic  power  (Watts) 

A,^,  . maximum  theoretical  effective  area  (square  meters) 

. maximum  theoretical  effective  length  (meters) 

Pe . numerical  distance  (vertical  polarization,  0=90°) 

Pm  . numerical  distance  (horizontal  polarization,  8= 90°) 

Angle^  .elevation  angle  of  maximum  directive  gain  (degrees) 

The  directivity  (D0)  ,  effective  isotropic  radiated  power 

(EIRP) ,  maximum  theoretical  effective  area  (A^J  ,  maximum 

theoretical  effective  length  (l,^)  ,  and  the  elevation  angle  of 

maximum  directive  gain  (Angled  are  all  determined  for  both 

the  0=0  and  0=tt/2  vertical  planes. 
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The  Mathcad  application  for  the  inverted  L  antenna 
requires  the  following  inputs: 


H . length  of  vertical  antenna  segment  (meters) 

1 . .length  of  horizontal  antenna  segment  (meters) 

f5 . operational  frequency  (Hertz) 

R . distance  from  antenna  (meters) 

6 r . relative  dielectric  constant  of  ground  plane 

a . conductivity  of  ground  plane 

d. . elevation  angle  index  (from  Table  3.1) 


The  combined  length  of  the  vertical  and  horizontal  antenna 
segments  (H+l)  is  used  to  determine  frequencies  which 
correspond  to  a  X/8,  X/4,  3X/8,  or  X/2  long  inverted  L 
antenna.  The  user  inputs  the  frequency  for  which  the 
radiation  parameters  discussed  in  Chapter  3  are  computed. 
Distance  from  the  antenna  (R)  must  meet  the  far- field 
requirements  of  Chapter  3.C.  The  elevation  angle  index  (d) 
sets  the  elevation  coordinate,  0,  for  which  a  horizontal 
radiation  pattern  is  determined.  Table  3.1  lists  possible 
indices  and  their  corresponding  elevation  angles  from  .285°  to 
88.857°  in  increments  of  about  2°.  Indices  between  those 
listed  can  be  used  to  interpolate  a  better  approximation  of  a 
desired  elevation. 

The  inverted  L's  radiated  electric  field  is  a  varying 
combination  of  3  and  $  components,  and  its  radiation  pattern 
varies  with  changes  in  <t>.  The  total  radiated  electric  field 
is  the  vector  sum  of  the  $  and  $  components .  The  radiated 
electric  field  for  the  inverted  L  antenna  is  obtained  in  a 
manner  analogous  to  that  used  for  the  vertical  dipole 
discussed  in  Chapter  4.  The  radiated  fields  of  each  antenna 
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segment  are  considered  separately,  and  superposition  is  used 
to  determine  the  total . radiated  electric  field.  The  vertical 
se ament  radiates  only  a  $  field  component,  but  the  horizontal 
segment  radiates  both  0  and  components.  The  electric  field 
equations  for  the  inverted  L  are  [Ref  6:  p.  169-170} 


r  nr*  *"*’  Av*jB--  T  K~*Bv 
r  sin  (6)  v  sir.  (6) 
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Equation  9.1  is  the  9  component  for  the  vertical  segment, 
equation  9.2  is  the  d  (vertically  polarized)  component  for  the 
horizontal  segment,  and  equation  9.3  is  the  <j>  (horizontally 
polarized)  component  for  the  horizontal  segment.  The  first 
two  terms  inside  the  brackets  of  each  equation  are  the  space 
wave  components  and  any  additional  terms  are  the  surface  wave 
components.  I0  is  taken  to  be  one,  since  a  sinusoidal  current 
input  with  maximum  of  unity  is  assumed. 

The  requested  inputs  are  used  to  calculate  the  following 

variables  using  a  constant  0  of  0=0  and  0=7 r/2  for  632  discrete 

values  of  9  which  are  equally  incremented  from  -tr/2  to  tt/2: 

h . half-length  of  the  dipole  (meters) 

X3 . wavelength  of  the  operational  frequency  (meters) 

13 . free  space  wavenumber  for  operational  frequency 

R.J . distance  from  antenna  to  field  point  (meters) 

Rr . distance  from  antenna  image  to  field  point  (meters) 

n . index  of  refraction 

Pe . complex  numerical  distance  (vertical  polarization) 

Pm  . complex  numerical  distance  (horizontal  polarization) 

rv . vertical  reflection  coefficient 

rh . horizontal  reflection  coefficient 

F . vertical  surface  wave  attenuation  factor 

Fm  . horizontal  surface  wave  attenuation  factor 

The  calculated  variables  are  used  to  evaluate  the  far- 
field  space  wave  and  surface  wave  for  the  discrete  values  of 
9  for  0=0  and  0=7r/2 .  The  total  space  wave  and  surface  wave 
distributions  are  determined  from  vector  addition  of 
corresponding  9  and  $>  components.  The  space  wave  and  surface 
wave  results  are  then  combined  for  corresponding  discrete 
values  of  6  to  obtain  the  total  radiated  electric  field 
distribution  for  the  0=0  and  0=tt/2  vertical  planes.  The  space 
wave,  surfac;  wave,  and  total  electric  field  results  are  then 
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normalized  with  respect  to  the  maximum  field  intensity  of 
each,  and  the  normalized  magnitudes  are  plotted  for  each 
discrete  8  to  depict  the  radiation  patterns.  The  inverted  L 
Mathcad  application  computes  the  space  wave,  surface  wave,  and 
total  electric  field  radiation  patterns  and  radiation 
parameters  in  the  <t>- 0  and  <p=ir/2  vertical  planes. 

The  variables  corresponding  to  the  selected  elevation 
angle  index  (d)  are  used  to  evaluate  the  radiated  electric 
field  components  for  the  space  wave  and  surface  wave  at  the 
fixed  elevation  angle  (0d)  as  <fi  varies  from  0  to  2ir  in  632 
equal  increments.  The  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  for  the  vertical  planes. 

To  find  the  contributions  of  the  d  and  #  electric  field 

components  to  the  total  average  radiated  power,  Prad,  equations 

3.0  and  3.9  are  used  to  integrate  equations  9.1,  9.2,  and  9.3 

over  the  hemispherical  Gaussian  surface  above  the  ground  plane 

at  a  fixed  radius  (R)  from  the  antenna.  The  sum  of  the 

integrals  is  the  antenna's  Prid.  With  the  discrete  values  of 

the  electric  field  and  total  average  radiated  power 

determined,  the  Mathcad  application  predicts  the  following 

radiation  parameters  using  the  equations  from  Chapter  3: 

Rfld  . radiation  resistance  (Ohms) 

D0 . directivity 

EIRP . effective  isotropic  radiated  power  (Watts) 

Artai  .....maximum  theoretical  effective  area  (square  meters) 

lnu,  . maximum  theoretical  effective  length  (meters) 

Pe . numerical  distance  (vertical  polarization,  0=90°) 

Pm  . numerical  distance  (horizontal  polarization,  8= 90°) 

Angle^  .elevation  angle  of  maximum  directive  gain  (degrees) 
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The  directivity  (D0)  ,  effective  isotropic  radiated  power 
(EIRP) ,  maximum  theoretical  effective  area  (A^,)  ,  maximum 
theoretical  effective  length  (lmai)  ,  and  the  elevation  angle  of 
maximum  directive  gain  (Anglemax)  are  all  determined  for  both 
the  <t>= 0  and  <p=ir/2  vertical  planes. 
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The  Mathcad  application  for  the  sloping  long-wire  antenna 
requires  the  following  inputs: 


1 . length  of  the  antenna  (meters) 

a . ..angle  of  antenna  with  horizontal  (degrees) 

f5 . operational  frequency  (Hertz) 

R . distance  from  antenna  (meters) 

€ r . relative  dielectric  constant  of  ground  plane 

a . conductivity  of  ground  plane 

d . elevation  angle  index  (from  Table  3.1) 


The  user  inputs  the  frequency  (f})  for  which  the  radiation 
parameters  discussed  in  Chapter  3  are  computed.  The  distance 
from  the  antenna  (R)  must  meet  the  far- field  requirements  of 
Chapter  3.C.  The  elevation  angle  index  (d)  sets  the  elevation 
coordinate,  6,  for  which  a  horizontal  radiation  pattern  is 
determined.  Table  3.1  lists  possible  indices  and  their 
corresponding  elevation  angles  from  .285°  to  88.857°  in 
increments  of  about  2°.  Indices  between  those  listed  can  be 
used  to  interpolate  \  better  estimate  of  a  desired  elevation. 

The  electric  f-eld  equations  for  the  sloping  long-wire 
are  given  by  [Ref  6:  p.  180-181] 
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e  '-jP* 
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The  radiated  electric  field  for  the  sloping  long-wire  is  22 
obtained  in  a  manner  analogous  to  the  vertical  dipole  equation 
discussed  in  Chapter  4.  The  long-wire's  radiated  electric 
field  is  a  varying  combination  of  d  and  $  components,  and  its 
radiation  patterns  vary  with  changes  in  <t>.  The  total  radiated 
electric  field  is  the  vector  sum  of  the  d  and  #  components. 
The  first  two  terms  inside  the  brackets  of  each  equation  are 
the  space  wave  components  and  the  additional  terms  are  the 
surface  wave  components.  I0  is  taken  to  be  one  since  a 
sinusoidal  current  input  with  a  maximum  of  unity  is  assumed. 

Requested  inputs  are  used  to  compute  the  variables  on  the 

next  page  using  a  constant  of  <p= 0  and  <p=ir/2  for  632  discrete 

values  of  6  which  are  equally  incremented  from  -ir/2  to  ir/2. 

Xs . wavelength  of  the  operational  frequency  (meters) 

0 . free  space  wavenumber  for  operational  frequency 

n . index  of  refraction 

P, . complex  numerical  distance  (vertical  polarization) 

Pm  . complex  numerical  distance  (horizontal  polarization) 

rv . vertical  reflection  coefficient 

rh . horizontal  reflection  coefficient 

F, . vertical  surface  wave  attenuation  factor 

Fm  . horizontal  surface  wave  attenuation  factor 


4» 
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The  calculated  variables  are  used  to  evaluate  the  far- 
field  space  wave  and  surface  wave  for  the  discrete  values  of 
6  for  0=0  and  0=tt/2.  The  total  space  wave  and  surface  wave 
distributions  are  determined  from  vector  addition  of 
corresponding  d  and  $  components.  The  space  wave  and  surface 
wave  results  are  then  combined  for  corresponding  discrete 
values  of  6  to  obtain  the  total  radiated  electric  field 
distribution  for  the  0=0  and  0=ir/2  vertical  planes.  The  space 
wave,  surface  wave,  and  total  electric  field  results  are  then 
normalized  with  respect  to  the  maximum  field  intensity  of 
each,  and  the  normalized  magnitudes  are  plotted  for  each 
discrete  6  to  depict  the  radiation  patterns.  The  sloping 
long-wire  Mathcad  application  computes  the  space  wave,  surface 
wave,  and  total  electric  field  radiation  patterns  and 
radiation  parameters  in  the  0=0  and  0=ir/2  vertical  planes. 

The  variables  corresponding  to  the  selected  elevation 
angle  index  (d)  are  used  to  evaluate  the  radiated  electric 
field  components  for  the  space  wave  and  surface  wave  at  the 
fixed  elevation  angle  (0d)  as  0  varies  from  0  to  2ir  in  632 
equal  increments.  The  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  for  the  vertical  planes. 

To  find  the  contributions  of  the  d  and  $  electric  field 
components  to  the  total  average  radiated  power,  P^,  equations 
3.8  and  3.9  are  used  to  integrate  equations  10.1  and  10.2  over 
the  hemispherical  Gaussian  surface  above  the  ground  plane  at 
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a  fixed  radius  <R)  from  the  antenna.  The  sum  of  the  integrals 
is  the  antenna's  Prad.  With  the  discrete  values  of  the  electric 
field  and  total  average  radiated  power  determined,  the  Mathcad 
application  predicts  the  following  radiation  parameters  using 
the  equations  from  Chapter  3: 


Rrad  . radiation  resistance  (Ohms) 

D0 . directivity 

EIRP. ... .effective  isotropic  radiated  power  (Watts) 

Amax  . maximum  theoretical  effective  area  (square  meters) 

l^,  .....maximum  theoretical  effective  length  (meters) 

Pc . numerical  distance  (vertical  polarization,  fl=90°) 

Pm  . numerical  distance  (horizontal  polarization,  0=90°) 


Angle^  .elevation  angle  of  maximum  directive  gain  (degrees) 
The  directivity  (D0)  ,  effective  isotropiic  radiated  power 
(EIRP) ,  maximum  theoretical  effective  area  (A^J  ,  maximum 
theoretical  effective  length  (la*,)  ,  and  the  elevation  angle  of 
maximum  directive  gain  (Angle^)  are  all  determined  for  both 
the  <p=0  and  4>=ir/2  vertical  planes. 
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XI.  THte  TERMINATED  SLOPING  V  ANTENNA 
The  orientation  of  the  terminated  sloping  V  antenna  is 
depicted  on  the  next  page  in  Figure  11.1,  where  1  is  length  of 
the  antenna  wires  {numbered  1  and  2)  ,  H0  is  the  feed  height 
above  ground,  H,  is  the  height  of  the  termination  points  above 
ground,  7  is  the  angle  subtended  by  the  wires,  rj  is  the  angle 
between  the  x  axis  and  the  projection  of  the  antenna  legs  in 
the  x-y  plane,  a  is  the  angle  between  the  plane  which  contains 
the  antenna  and  the  x-y  plane,  R  is  the  radial  coordinate,  9 
is  the  elevation  coordinate,  and  <f>  is  the  azimuth  coordinate. 

The  Mathcad  application  for  the  sloping  V  antenna  requires 
the  following  inputs: 


1 . length  of  the  antenna  legs  (meters) 

Hq . feed  height  above  ground  (meters) 

H, . height  of  termination  points  above  ground  (meters) 

7 . angle  subtended  by  the  antenna  legs  (degrees) 

fj . operational  frequency  (Hertz) 

R . distance  from  antenna  (meters) 

er . relative  dielectric  constant  of  ground  plane 

a . conductivity  of  ground  plane 

d . elevation  angle  index  (from  Table  3.1) 


The  user  inputs  the  frequency  (f3)  for  which  the  radiation 
parameters  discussed  in  Chapter  3  are  computed.  The  feed 
height  (Hq)  and  termination  height  (H,)  can  be  any  value 
greacer  than  zero.  The  angle  a  is  found  from  Hq,  H,,  wire 
length  (1),  and  the  half-angle  between  the  wires  ( 7/2 )  by 

casin'1  — -1  —  (11.1) 
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FIGURE  11.1:  Spatial  orientation  of  the  terminated 
sloping  V  antenna  for  its  corresponding  Mathcad 
application. 


The  angle  t\  is  then  found  from  a  and  7/2  by 


ri=tan' 


Hi) 

cos (a) 


(11.2) 


Distance  from  the  antenna  (R)  must  meet  the  far- field 
requirements  of  Chapter  3.C.  The  elevation  angle  index  (d) 
sets  the  9  coordinate  for  which  a  horizontal  radiation  pattern 
is  determined.  Table  3.1  lists  possible  indices  and  their 
corresponding  elevation  angles  from  .285°  to  88.857°  in 


increments  of  about  2°.  Indices  between  those  listed  can  be 


used  to  interpolate  a  better  estimate  of  a  desired  elevation. 

The  radiated  electric  field  for  the  sloping  V  antenna  is 
obtained  in  a  manner  analogous  to  that  used  for  the  vertical 
dipole  discussed  in  Chapter  4.  The  sloping  V  electric  field 
equations  are  given  by  [Ref  6:  p.  185-186] 
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The  sloping  V's  radiated  electric  field  is  a  varying 
combination  of  9  and  $  components,  and  its  radiation  pattern 
varies  with  changes  in  <p.  The  total  radiated  electric  field 
is  the  vector  sum  of  the  d  and  $  components.  The  fields  of 
each  antenna  leg  are  considered  separately,  and  superposition 
is  used  to  determine  the  total  electric  field.  Equations  11.3 
and  11.4  are  the  radiated  9  and  $  components,  respectively  for 
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number  one,  and  equations  11.5  and  11.6  are  the  9  and  0 
components,  respectively  for  wire  number  two.  The  first  two 
terms  inside  each  set  of  brackets  represent  the  space  wave 
components  and  any  additional  terms  represent  the  surface  wave 
components.  I0  is  taken  to  be  one,  since  a  sinusoidal  current 
input  with  a  maximum  of  unity  is  assumed. 

The  requested  inputs  are  used  to  calculate  the  following 

variables  using  a  constant  0  of  0=0  and  0= rr/2  for  632  discrete 

values  of  6  which  are  equally  incremented  from  -ir/2  to  tt/2: 

\3 . wavelength  of  the  operational  frequency  (meters) 

/? . free  space  wavenumber  for  operational  frequency 

n . index  of  refraction 

Pc . complex  numerical  distance  (vertical  polarization) 

Pm  . complex  numerical  distance  (horizontal  polarization) 

rv . vertical  reflection  coefficient 

rh . horizontal  reflection  coefficient 

F . vertical  surface  wave  attenuation  factor 

Fm  . horizontal  surface  wave  attenuation  factor 

The  calculated  variables  are  used  to  evaluate  the  far- 
field  space  wave  and  surface  wave  for  the  discrete  values  of 
9  for  0=0  and  0= tt/2  .  The  total  space  wave  and  surface  wave 
distributions  are  determined  from  vector  addition  of 
corresponding  d  and  0  components.  The  space  wave  and  surface 
wave  results  are  then  combined  for  co; responding  discrete 
values  of  9  to  obtain  the  total  radiated  electric  field 
distribution  for  the  0=0  and  0=7r/2  vertical  planes.  The  space 
wave,  surface  wave,  and  total  electric  field  results  are  then 
normalized  with  respect  to  the  maximum  field  intensity  of 
each,  and  the  normalized  magnitudes  are  plotted  for  each 
discrete  6  to  depict  the  radiation  patterns.  The  sloping  V 
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Mathcad  application  computes  the  space  wave,  surface  wave,  and 
total  electric  field  radiation  patterns  and  radiation 
parameters  in  the  tf>=0  and  tt / 2  vertical  planes. 

The  variables  corresponding  to  the  selected  elevation 
angle  index  id)  are  used  to  evaluate  the  radiated  electric 
field  components  for  the  space  wave  and  surface  wave  at  the 
fixed  elevition  angle  (0 j)  as  <t>  varies  from  0  to  2 tt  in  632 
equal  increments.  The  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  for  the  vertical  planes. 

To  find  the  contributions  of  the  §  and  $  electric  field 

components  to  the  total  average  radiated  power,  Pri(1,  equations 

3.8  and  3.9  are  used  to  integrate  equations  11.3,  11.4,  11.5, 

and  11.6  over  the  hemispherical  Gaussian  surface  above  the 

ground  plane  at  a  fixed  radius  (R)  from  the  antenna.  The  sum 

of  the  integrals  is  the  antenna's  Prtd.  With  the  discrete 

values  of  the  electric  field  and  total  average  radiated  power 

determined,  the  Mathcad  application  predicts  the  following 

radiation  parameters  using  the  equations  from  Chapter  3: 

R,*  . radiation  resistance  (Ohms) 

D0 . directivity 

EIRP . effective  isotropic  radiated  power  (Watts) 

A**,  . maximum  theoretical  effective  area  (square  meters) 

l*,,  .....maximum  theoretical  effective  length  (meters) 

P,  . numerical  distance  (vertical  polarization,  8-90°) 

P„  . numerical  distance  (horizontal  polarization,  0=> 90°) 

Angle,™,  .elevation  angle  of  maximum  directive  gain  (degrees) 

The  directivity  (D0)  ,  effective  isotropic  radiated  power 

(EIRP),  maximum  theoretical  effective  area  (A^,)  ,  maximum 
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theoretical  effective  length  (lm„),  and  the  elevation  angle  of 
maximum  directive  gain  (Anglemj<)  are  all  determined  for  both 
the  <t>=0  and  <£  =  ? r/2  vertical  planes. 
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FIGURE  12. li  Spatial  orientation  of  the  side- loaded 
vertical  half  rhombic  antenna  for  its  Mathcad  application. 
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The  Mathcad  application  for  the  vertical  half  rhombic 
antenna  requires  the  following  inputs: 


1 . length  of  the  antenna  legs  (meters) 

a . angle  of  antenna  wires  with  horizontal  (degrees) 

f5 . operational  frequency  (Hertz) 

R . distance  from  antenna  (meters) 

€r . relative  dielectric  constant  of  ground  plane 

a... . conductivity  of  ground  plane 

d . elevation  angle  index  (from  Table  3.1) 


The  user  inputs  the  frequency  (f5)  for  which  the  radiation 
parameters  discussed  in  Chapter  3  are  computed.  Distance  from 
the  antenna  (R)  must  meet  the  far -field  requirements  of 
Chapter  3.C.  The  elevation  angle  index  (d)  sets  the  6 
coordinate  for  which  a  horizontal  radiation  pattern  is 
determined.  Table  3.1  lists  possible  indices  and  their 
corresponding  elevation  angles  from  .285°  to  88.857°  in 
increments  of  about  2°.  Indices  between  those  listed  can  be 
used  to  interpolate  a  better  approximation  of  a  desired 
elevation. 

The  vertical  half  rhombic's  radiated  electric  field  is  a 
combination  of  0  and  components,  and  its  radiation  pattern 
varies  with  changes  in  <f>.  The  total  radiated  electric  field 
is  the  vector  sum  of  the  9  and  $  components.  The  radiated 
electric  field  for  the  vertical  half  rhombic  antenna  is 
obtained  in  a  manner  analogous  to  that  used  for  the  vertical 
dipole  discussed  in  Chapter  4.  The  radiated  fields  of  each 
antenna  leg  are  considered  separately,  and  superposition  is 
used  to  determine  the  total  electric  field.  The  electric 
field  equations  for  the  half  rhombic  are  [Ref  6:  p.  198-200] 
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where  U,  and  U2  are  given  by 

t/^l-cos  (7^  )  y2  =  l-cos(72) 

cos  (7:)  =cos  (0)  sin  (a)  +sin  (0)  cos  (a)  cos  (4>) 

cos  (72)  =sin(0)  cos  (a)  cos  (<J>)  -cos  (0)  sin  (a) 

Equations  12.1  and  12.2  are  the  §  (vertically  polarized)  and 

$  (horizontally  polarized)  components,  respectively  for  leg 

number  one,  and  equations  12.3  and  12.4  are  the  0  and  $ 

components,  respectively  for  leg  number  two.  The  first  two 

terms  inside  each  set  of  brackets  are  the  space  wave 

components  and  any  additional  terms  are  the  surface  wave 

components.  I0  is  taken  to  be  one,  since  a  sinusoidal  current 

input  with  a  maximum  of  unity  is  assumed.  The  requested 

inputs  are  used  to  calculate  the  following  variables  using  a 

constant  4>  of  </>= 0  and  <£=7r/2  for  632  discrete  values  of  0  which 

are  equally  incremented  from  -ir/2  to  7t/2: 

X 5 . wavelength  of  the  operational  frequency  (meters) 

,3 . free  space  wavenumber  for  operational  frequency 

n . index  of  refraction 

P, . complex  numerical  distance  (vertical  polarization) 

Pm  . complex  numerical  distance  (horizontal  polarization) 

rv . vertical  reflection  coefficient 

rh . horizontal  reflection  Coefficient 

F, . vertical  surface  wave  attenuation  factor 

Fm  . horizontal  surface  wave  attenuation  factor 

I 

The  calculated  variables  are  used  to  evaluate  the  far- 

\ 

field  space  wave  and  surface  wave  for  the  discrete  values  of 
0  for  <£»0  and  <£=7r/2.  The  total  space  wave  and  surface  wave 
distributions  are  determined  from  vector  addition  of 
corresponding  §  and  $  components.  The  space  wave  and  surface 
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wave  results  are  then  combined  for  corresponding  discrete 
values  of  $  to  obtain  the  total  radiated  electric  field 
distribution  for  the  <p=0  and  0=7r/2  vertical  planes.  The  space 
wave,  surface  wave,  and  total  electric  field  results  are  then 
normalized  with  respect  to  the  maximum  field  intensity  of 
each,  and  the  normalized  magnitudes  are  plotted  for  each 
discrete  6  to  depict  the  radiation  patterns.  The  vertical 
half  rhombic  Mathcad  application  computes  the  space  wave, 
surface  wave,  and  total  electric  field  radiation  patterns  and 
radiation  parameters  in  the  <f>=  0  and  <f>=n/2  vertical  planes. 

The  variables  corresponding  to  the  selected  elevation 
angle  index  (d)  are  used  to  evaluate  the  radiated  electric 
field  components  for  the  space  wave  and  surface  wave  at  the 
fixed  elevation  angle  (0d)  as  <fi  varies  from  0  to  2n  in  632 
equal  increments.  The  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  for  the  vertical  planes. 

To  find  the  contributions  of  the  9  and  $  electric  field 
components  to  the  total  average  radiated  power,  P„d,  equations 
3.8  and  3.9  are  used  to  integrate  equations  12.1,  12.2,  12.3, 
and  12.4  over  the  hemispherical  Gaussian  surface  above  the 
ground  plane  at  a  fixed  radius  (R)  from  the  antenna.  The  sum 
of  the  integrals  is  the  antenna's  PrKl.  With  the  discrete 
values  of  the  electric  field  and  total  average  radiated  power 
determined,  the  Mathcad  application  predicts  the  following 
radiation  parameters  using  the  equations  from  Chapter  3: 
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. radiation  resistance  (Ohms) 

D0  . . directivity 

EIRP . effective  isotropic  radiated  power  (Watts) 

A™,  . maximum  theoretical  effective  area  (square  meters) 

lma,  . maximum  theoretical  effective  length  (meters) 

P . numerical  distance  (vertical  polarization,  0  =  90°) 

Pm  . numerical  distance  (horizontal  polarization,  0=90°) 


Angle,^,  .elevation  angle  of  maximum  directive  gain  (degrees) 
The  directivity  (D0)  ,  effective  isotropic  radiated  power 
(EIRP) ,  maximum  theoretical  effective  area  (A^,)  ,  maximum 
theoretical  effective  length  (I™,),  and  the  elevation  angle  of 
maximum  directive  gain  (Angle™,)  are  all  determined  for  both 
the  0=0  and  0=tt/2  vertical  planes. 
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XIII.  THE  TERMINATED  SLOPING  OR  HORIZONTAL  RHOMBIC  ANTENNA 

The  orientation  of  the  terminated  sloping  or  horizontal 
rhombic  antenna  is  depicted  on  the  next  page  in  Figure  13.1, 
where  1  i3  the  length  of  the  antenna  wires  (numbered  1  to  4), 
Ho  is  the  feed  height  above  ground,  H2  is  the  height  of  the 
termination  point  above  ground,  y  is  the  angle  subtended  by 
the  feed  wires  (1  and  2),  jj  is  the  angle  between  the  x  axis 
and  the  projection  of  the  feed  wires  in  the  x-y  plane,  a  is 
the  angle  between  the  plane  which  contains  the  antenna  and  the 
x-y  plane,  R  is  the  radial  coordinate,  6  is  the  elevation 
coordinate,  and  0  is  the  azimuth  coordinate. 

The  Mathcad  application  for  the  rhombic  antenna  requires 
the  following  inputs: 


1 . length  of  the  antenna  legs  (meters) 

Hq . feed  height  above  ground  (meters) 

H2 . height  of  termination  point  above  ground  (meters) 

y . angle  subtended  by  the  antenna  legs  (degrees) 

f3 . operational  frequency  (Hertz) 

R . distance  from  antenna  (meters) 

€, . relative  dielectric  constant  of  ground  plane 

a. ..... . .conductivity  of  ground  plane 

d . elevation  angle  index  (from  Table  3.1) 


The  feed  he  .ght  (Hq)  and  termination  height  (H2)  can  be  any 
value  greater  than  zero.  The  angle  a  is  found  from  Hq,  H;, 
wire  length  (1),  and  the  half-angle  between  the  wires  iy/2)  by 
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or  horizontal  rhombic  antenna  for  its  Mathcad  application. 


The  angle  »?  is  then  found  from  a  and  y/2  by 


tj  =  tan 


-U 


tan(i) 


(13.2) 


cos (a) 

The  user  inputs  the  frequency  (f5)  for  which  the  radiation 
parameters  discussed  in  Chapter  3  are  computed.  Distance  from 
the  antenna  (R)  must  meet  the  far- field  requirements  of 
Chapter  3.C.  The  elevation  angle  index  (d)  sets  the  6 
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coordinate  for  which  a  horizontal  radiation  pattern  is 
determined.  Table  3.1  lists  possible  indices  and  their 
corresponding  elevation  angles  from  .285°  to  88.857°  in 
increments  of  about  2°.  Indices  between  those  listed  can  be 
used  to  interpolate  a  better  estimate  of  a  desired  elevation. 

The  radiated  electric  field  for  the  rhombic  antenna  is 
obtained  in  a  manner  analogous  to  that  for  the  vertical  dipole 
discussed  in  Chapter  4.  The  radiated  field  of  each  wire  is 
considered  separately,  and  superposition  is  used  to  determine 
the  total  electric  field.  The  rhombic's  electric  field 
equations  are  [Ref  6:  p.  183-191] 
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The  rhombic's  radiated  electric  field  is  a  varying  combination 
of  9  and  $  components,  and  its  radiation  pattern  varies  with 
changes  in  <f>.  The  total  radiated  electric  field  is  the  vector 
sum  of  the  5  and  $  components.  Equations  13.3  and  13.4  are 
the  d  (vertically  polarized)  and  $  (horizontally  polarized) 
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components,  respectively  for  wire  one,  and  equations  13.5  and 
13.6  are  the  8  and  $  components,  respectively  for  wire  number 
two.  Equations  13.7  and  13.8  are  the  analogous  equations  for 
wire  three,  and  13.9  and  13.10  are  those  for  wire  four.  The 
first  two  terms  inside  each  set  of  brackets  represent  are  the 
space  wave  components,  and  any  additional  terms  represent  the 
surface  wave  components.  I0  is  taken  to  be  one,  since  a 
sinusoidal  current  input  with  a  maximum  of  unity  is  assumed. 

The  requested  inputs  are  used  to  calculate  the  following 
variables  using  a  constant  0  of  0=0  and  0= rr/2  for  632  discrete 
values  of  8  which  are  equally  incremented  from  -tt/2  tc  rr/2: 

X3 . wavelength  of  the  operational  frequency  (meters) 

. free  space  wavenumber  for  operational  frequency 

n . index  of  refraction 

Pe . complex  numerical  distance  (vertical  polarization) 

Pm  . complex  numerical  distance  (horizontal  polarization) 

rv . vertical  reflection  coefficient 

rh  ..... .horizontal  reflection  coefficient 

Fe . vertical  surface  wave  attenuation  factor 

Fm  . horizontal  surface  wave  attenuation  factor 

The  calculated  variables  are  used  to  evaluate  the  far- 
field  space  wave  ard  surface  wave  for  the  discrete  values  of 
6  for  0=0  and  0=7r/2.  The  total  space  wave  and  surface  wave 
distributions  are  determined  from  vector  addition  of 
corresponding  §  and  $  components.  The  space  wave  and  surface 
wave  results  are  then  combined  for  corresponding  discrete 
values  of  8  to  obtain  the  total  radiated  electric  field 
distribution  for  the  0=0  and  0=7r/2  vertical  planes.  The  space 
wave,  surface  wave,  and  total  electric  field  results  are  then 
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normalized  with  respect  to  the  maximum  field  intensity  of 
each,  and  the  normalized  magnitudes  are  plotted  for  each 
discrete  6  to  depict  the  radiation  patterns.  The  rhombic 
antenna  Mathcad  application  computes  the  space  wave,  surface 
wave,  and  total  electric  field  radiation  patterns  and 
radiation  parameters  in  the  <p=0  and  <£= 7r/2  vertical  planes. 

The  variables  corresponding  to  the  selected  elevation 
angle  index  (d)  are  used  to  evaluate  the  radiated  electric 
field  components  for  the  space  wave  and  surface  wave  at  the 
selected  elevation  angle  (0d)  as  <p  varies  from  0  to  2rr  in  632 
equal  increments.  The  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  for  the  vertical  planes. 

To  find  the  contributions  of  the  d  and  #  electric  field 
components  to  the  total  average  radiated  power,  Prad,  equations 
3.8  and  3.9  are  used  to  integrate  equations  13.3  to  13.10  over 
the  hemispherical  Gaussian  surface  above  the  ground  plane  at 
a  fixed  radius  (R)  from  the  antenna.  The  sum  of  the  integrals 
is  the  antenna's  P^,.  With  the  discrete  values  of  the  electric 
field  and  total  average  radiated  power  determined,  the  Mathcad 
application  predicts  the  following  radiation  parameters  using 
the  equations  from  Chapter  3: 


Rrld  . radiation' resistance  (Ohms) 

D0 . directivity 

EIRP . effective  isotropic  radiated  power  (Watts) 

Piau  . maximum  theoretical  effective  area  (square  meters) 

Inuu,  . maximum  theoretical  effective  length  (meters) 

Pe . numerical  distance  (vertical  polarization,  6=30°) 

P„  . numerical  distance  (horizontal  polarization,  6=90°) 


Angle^  .elevation  angle  of  maximum  directive  gain  (degrees) 
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The  directivity  (D0)  ,  effective  isotropic  radiated  power 
(EIRP) ,  maximum  theoretical  effective  area  (A^,)  ,  maximum 
theoretical  effective  length  (lm„) ,  and  the  elevation  angle  of 
maximum  directive  gain  (Angle,.,,)  are  all  determined  for  both 
the  4>=Q  and  4>=tt/2  vertical  planes. 


92 


XIV.  THE  TERMINATED  SLOPING  DOUBLE  RHOMBOID  ANTENNA 

The  orientation  of  the  terminated  sloping  double  rhomboid 
antenna  is  depicted  on  the  next  page  in  Figure  14.1,  where  1, 
is  the  length  of  the  four  short  antenna  wires  (A,  D,  E,  and 
H)  ,  1:  is  the  length  of  the  four  long  antenna  wires  (B,  C,  F, 
and  G)  ,  lj  is  the  length  of  the  radials  which  connect  the  feed 
to  the  two  termination  points,  Hq  is  the  feed  height  above 
ground,  H3  is  the  height  of  the  termination  points  above 
ground,  y  is  the  angle  between  by  the  radials  which  connect 
the  feed  to  the  termination  points,  6  is  the  angle  between  the 
termination  radials  and  the  short  feed  wires,  77  is  the  angle 
between  the  termination  radials  and  the  long  feed  wires,  a  is 
the  angle  between  the  plane  which  contains  the  antenna  and  the 
x-y  plane,  R  is  the  radial  coordinate,  6  is  the  elevation 
coordinate,  and  <f>  is  the  azimuth  coordinate. 

The  Mathcad  application  for  the  double  rhomboid  antenna 
requires  the  following  inputs: 


1,  . length  of  the  short  antenna  wires  (meters) 

1,  . length  of  the  long  antenna  wires  (meters) 

13  . length  of  the  termination  radials  (meters) 

Hq . feed  height  above  ground  (meters)  \ 

H3 . height  of  termination  points  above  ground  (meters) 

7 . angle  subtended  by  termination  radials  (^egrees) 

fj . operational  frequency  (Hertz)  \ 

R . ...distance  from  antenna  (meters)  \ 

e, . relative  dielectric  constant  of  ground  plane 

a . conductivity  of  ground  plane 

d . elevation  angle  index  (from  Table  3.1) 


The  feed  height  (Hq)  and  termination  height  (H3)  can  be  any 
value  greater  than  zero.  The  length  of  the  radials  to  the 
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TOP  VIEW 


B 


1 


Sloping  Double  Rhomboid  Orientation 


FIGURE  14.1s  Spatial  orientation  of  the  terminated  sloping 
double  rhomboid  antenna  for  its  Mathcad  application. 
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termination  points  must  be  less  than  the  combined  lengths  of 
a  long  and  short  antenna  wire  (13  <  1,  +  1,)  .  The  angle  r;  is 
found  from  the  law  of  sines  and  cosines  by 


Ji,  . 

tl=sin'i|-j±sm  cos 


(i:)  2+  (i3) 2-  (1,) 2 

2 1A 


(14.1) 


The  angle  5  is  also  found  from  the  law  of  cosines  by 


,  (11) 2+ (13) 2- (12) 

8fs  — —  TTT7 - 


(14.2) 


The  angle  a  is  then  found  from  H0,  H3,  1,,  12,  and  the  half¬ 


angle  between  the  wires  ( y/2 )  by 


(14.3) 

The  user  inputs  the  frequency  (f3)  for  which  the  radiation 
parameters  discussed  in  Chapter  3  are  computed.  Distance  from 
the  antenna  (R)  must  meet  the  far- field  requirements  of 
Chapter  3.C,  and  the  elevation  angle  index  (d)  sets  the  6 
coordinate  for  which  a  horizontal  radiation  pattern  is 
determined.  Table  '3.1  lists  possible  indices  and  their 
corresponding  elevation  angles  from  .285°  to  88.857°  in 
increments  of  about  2°.  Indices  between  those  listed  can  be 
used  to  interpolate  a  better  estimate  of  a  desired  elevation. 

The  radiated  electric  field  for  the  double  rhomboid  is 
obtained  in  a  manner  analogous  to  that  used  for  the  vertical 
dipole  discussed  in  Chapter  4.  The  electric  field  equations 
for  the  double  rhomboid  are  given  by  [Ref  6:  p.  238-244] 


a=sm* 


H2-H0 


11-cos|-|+6|  +  12-cos|-^+ti 


95 


^•.a-3  0  Io 
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H  2 
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(14.4) 
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-sin (a) sin (0) 
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■cos (0 


e-JP* 
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(14.5) 


,  /-,  r  %  p  i-e~jgJ^5  -jap^cosio 

«  rr  '  ■*■  1  A'  *  m  --  ^ 


£4 


£4 


Fe.s=30Ic 
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(7, 


(14.6) 


r  o-^2Phocos(0)  a  M  pi  \  r  l-e"'7'p'!jC,‘  „ -j^3«„cos(0) 
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where  U,,  U:,  U3,  U4,  UJf  U6,  U7,  and  U8  are  given  by 

J71  =  l-cos{71)  t/,  =  l-cos  (Y2) 

C/3  =  1-cos  (T3)  f/4  =  1-cos  (T4) 

C/5  =  1-cos(75)  f/6  =  l-cos  (Y6) 

Cf7  =  l-cos  (T7)  lr8  =  l-cos  (Y8)  where 

cos  (7^  =cos  (0)sin(a).  +  sin(0)cos(o)sin(<j>-o-6) 
cos  (72)  =cos  (0)  sin  (a )  +sin (0) cos (a )  sin  ( <J>- or  -r) ) 
cos  (73)  =cos  (0)  sin(a)  +  sin{0)  cos  (a)  sin(4>+a+ri) 
coc(74)  =cos  (0)  sin  (a)  +sin(0)  cos  (a)  sin(4>+a+6) 
cos(Y5)=sin(0)cos(a)sin(<{>-a-6)  -cos  (0)  sin  (a) 
cos  (Y6)  =sin(0)  cos  (a)  sin  (<j>-a-r| )  -cos  (0)  sin  (a) 
cos  ( Y7)  =sin  (0)  cos  (a)  sin  (<}>+a+T| )  -cos  (0)  sin  (a) 
cos  (Y8)  =sin(0)  cos  (a)  sin  (4>+e  +  8)  -cos  (0)  sin  (a) 


The  double  rhomboid's  radiated  electric  field  is  a  combination 
of  §  and  $  components,  and  its  radiation  pattern  varies  with 
changes  in  <t> .  The  total  electric  field  is  the  vector  sum  of 
the  d  and  $  components.  The  radiated  fields  of  each  antenna 
wire  are  considered  separately,  and  superposition  is  used  to 


determine  the  total  electric  field.  Equations  14.4  and  14.5 
are  the  d  (vertically  polarized)  and  $  (horizontally 


polarized)  components,  respectively  for  wire  A,  and  equations 
14.6  and  14.7  are  the  d  and  $  components,  respectively  for 
wire  B.  Equations  14.8  and  14.9  are  the  analogous  equations 
for  wire  C,  14.10  and  14.11  those  for  wire  D,  14.12  and  14.13 


those  for  wire  E,  14.14  and  14.15  those  for  wire  F,  14.16  and 


14.17  those  for  wire  G,  and  14.18  and  14.19  those  for  wire  H. 
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The  first  two  terms  inside  each  set  of  brackets  represent  the 
space  wave  components  and  any  additional  terms  represent  the 
surface  wave  components.  I0  is  taken  to  be  one  since  a 
sinusoidal  current  input  with  a  maximum  of  unity  is  assumed. 

The  requested  inputs  are  used  to  calculate  the  following 

variables  using  a  constant  0  of  0=0  and  0=tt/2  for  632  discrete 

values  of  6  which  are  equally  incremented  from  -7r/2  to  r/2: 

X5 . wavelength  of  the  operational  frequency  (meters) 

/3 . free  space  wavenumber  for  operational  frequency 

n . inuex  of  refraction 

Pe . complex  numerical  distance  (vertical  polarization) 

P,,,  . complex  numerical  distance  (horizontal  polarization) 

rv . vertical  reflection  coefficient 

rh . horizontal  reflection  coefficient 

Fe . vertical  surface  wave  attenuation  factor 

Fm  . horizontal  surface  wave  attenuation  factor 

The  calculated  variables  are  used  to  evaluate  the  far- 
field  space  wave  and  surface  wave  for  the  discrete  values  of 
6  with  0=0  and  0=7r/2 .  The  total  space  wave  and  surface  wave 
distributions  are  determined  from  vector  addition  of 
corresponding  d  and  $  components.  The  space  wave  and  surface 
wave  results  are  then  combined  for  corresponding  discrete 
values  of  6  to  obtain  the  total  radiated  electric  field 
distribution  for  the  0=0  and  0=m/2  vertical  planes.  The  space 
wave,  surface  wave,  and  total  electric  field  results  are  then 
normalized  with  respect  to  the  maximum  field  intensity  of 
each,  and  the  normalized  magnitudes  are  plotted  for  each 
discrete  6  to  depict  the  radiation  patterns.  The  Mathcad 
double  rhomboid  antenna  application  computes  the  space  wave, 
surface  wave,  and  total  electric  field  radiation  patterns  and 
radiation  parameters  in  the  0=0  and  0=tt/2  vertical  planes. 
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The  variables  corresponding  to  the  selected  elevation 
angle  index  (d)  are  used  to  evaluate  the  radiated  electric 
field  components  for  the  space  wave  and  surface  wave  at  the 
fixed  elevation  angle  (Cj)  as  0  varies  from  0.  to  27 r  in  632 
equal  increments.  The  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  for  the  vertical  planes. 

To  find  the  contributions  of  the  d  and  $  electric  field 

components  to  the  total  average  radiated  power,  PraJ,  equations 

3.8  and  3.9  are  used  to  integrate  equations  14.4  to  14.19  over 

the  hemispherical  Gaussian  surface  above  the  ground  plane  at 

a  fixed  radius  (R)  from  the  antenna.  The  sum  of  the  integrals 

is  t.ie  antenna  '3  Prad-  with  the  discrete  values  of  the  electric 

field  and  total  average  radiated  power  determined,  the  Mathcad 

application  predicts  the  following  radiation  parameters  using 

the  equations  from  Chapter  3: 

Rrad  . radiation  resistance  (Ohms) 

D0 . directivity 

EIRP . effective  isotropic  radiated  power  (Watts) 

A^,  . maximum  theoretical  effective  area  (square  meters) 

lou,  . maximum  theoretical  effective  length  (meters) 

Pe . numerical  distance  (vertical  polarization,  6= 90°) 

Pm  . numerical  distance  (horizontal  polarization,  0=90°) 

Angle,^  .elevation  angle  of  maximum  directive  gain  (degrees) 

The  directivity  (D0)  ,  effective  isotropic  radiated  power 

(EIRP) ,  maximum  theoretical  effective  area  (A^,)  ,  maximum 

theoretical  effective  length  (1^) ,  and  the  elevation  angle  of 

maximum  directive  gain  (Angle^)  are  all  determined  for  both 

the  <£=0  and  <p=ir/2  vertical  planes. 
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XV. 


THE  VERTICALLY  POLARIZED  LOG- PERIODIC  DIPOLE  ARRAY 


The  orientation  of  the  vertically  polarized  log-periodic 
dipole  array  is  depicted  in  Figure  15.1,  where  N  is  the  number 
of  dipoles,  n  is  an  index  equal  to  0,1,2..N-1,  ln  is  the 
length  of  the  n^+l  dipole,  d„  is  the  distance  between  the  n^+l 
and  ^+2  dipole  elements,  Ho  is  the  height  of  the  antenna  feed 
above  ground,  H„  is  the  height  of  the  n^+l  dipole  center,  ^  is 
the  angle  between  the  vertical  axis  and  the  array's  center 
axis,  a,  is  the  angle  between  the  center  axis  and  the  line 
which  connects  the  upper  tips  of  the  elements,  a,  is  the  angle 
between  the  center  axis  and  the  line  which  connects  the  lower 
tips  of  the  elements,  a3  is  the  angle  between  the  horizontal 
and  the  line  connecting  the  lower  tips  of  the  elements,  R  is 
the  radial  coordinate,  6  is  the  elevation  coordinate,  and  0  is 
the  azimuth  coordinate. 

The  Mathcad  application  for  the  vertically  polarized  log- 
periodic  array  requires  the  following  inputs: 


N . number  of  dipole  elements 

10  . length  of  the  shortest  dipole  (meters) 

1,  . length  of  the  second  shortest  dipole  (meters) 

do  . . . distance  between  first  two  elements  (meters) 

rado . radius  of  the  shortest  dipole  (meters) 

Ho . feed  height  above  ground  plane  (meters) 

V . ...angle  from  vertical  axis  to  array  axis  (degrees) 

fj . operational  frequency  (Hertz) 

R . distance  from  array  (meters) 

€, . relative  dielectric  constant  of  ground  plane 

<r . conductivity  of  ground  plane 

w. . . elevation  angle  index  (from  Table  3.2) 

ADM . transmission  line  characteristic  admittance 

TIMP . termination  impedance  connected  to  longest  dipole 


The  user  inputs  the  frequency  for  which  the  radiation 
parameters  discussed  in  Chapter  3  are  computed.  The  lengths 
of  the  first  two  elements  (10  and  1,)  are  used  to  calculate  the 
log-periodic  relationship  (r)  given  by  [Ref  10:  p.  317] 

t  =  A  (n=0 , 1 , 2 _ N-2.N-1)  (15.1) 

in*  i 


The  remaining  element  lengths  (12  to  1N.,)  are  then  found  from 


(15.2) 


As  noted  in  equation  15.1,  indices  start  at  zero  instead  of  1, 
so  an  index  simply  refers  to  the  index- plus -one  position  in 
terms  of  successive  values  for  a  given  parameter.  This  is  to 
maintain  continuity  with  Mathcad  which  also  indexes  from  zero. 
The  radius  of  the  first  element  (rado)  is  used  to  find  the 
remaining  radii  by  substituting  rad„  and  rad„+,  for  ln  and  ln+, 
in  equation  15.2.  The  angle  between  the  vertical  axis  and  the 
array  axis  W)  can  be  between  0°  and  90°.  The  angles  oi:  and 
a3  are  found  by  substituting 


a2  +  «3  =  f-7 
into 

dn  _ _ 1-T _ 

in,!  sin(a2+a3)  -tan(a3)  cos  (a2  +  a3) 
which  yields 


(15.3) 

(15.4) 


a3=tan'1| 


d„cos(|-T) 
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The  angle  a,  is  then  found  from  [Ref  10:  pp.  351-354] 


sin(ax)  sin(a2) 
cos  (a1  +  a2  +  a3)  cos(a3) 


(15.5) 


a1  =  tan'1 


sin (cc2)  cos  (ot2  +  a3) 
cos(a3)  +  sin(a2)sin(a2  +  a3) 


The  remaining  separations  ( d,  to  dN.2)  are  then  found  from 


equation  15.4.  The  dipole  half-lengths  (hj  are  found  from 
ln/2,  the  feed  height  (H0)  must  be  at  least  equal  to  the  half- 
length  of  the  longest  dipole  (hN-1)  ,  and  distance  from  the 
array  (R)  must  meet  the  far- field  requirements  of  Chapter  3.C. 
The  elevation  angle  index  (w)  sets  the  6  coordinate  for  which 
a  horizontal  radiation  pattern  is  determined.  Table  3.2  lists 
possible  indices  and  their  corresponding  elevation  angles  from 
.57°  to  89.1°  in  increments  of  about  2.3°.  Indices  between 


those  listed  can  be  used  to  interpolate  a  better  approximation 
of  a  desired  elevation.  The  termination  impedance  (TIMP)  is 
assumed  to  be  connected  to  the  array's  center  axis  opposite 
the  feed  at  a  distance  of  hN.,/2  from  the  longest  element. 
Typical  termination  impedances  are  the  transmission  line's 
matched  impedance  (complex  conjugate)  or  a  short  circuit 
(zero) .  The  characteristic  admittance  (ADM)  is  simply  the 
inverse  of  the  transmission  line's  characteristic  impedance. 

If  computed  lengths,  separations  and  radii  are  not 
desired,  they  can  be  entered  manually  as  described  on  the  page 
two  of  the  application.  Consult  a  Mathcad  manual  to  be 
certain  the  entries  are  made  correctly. 
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The  open  circuit  impedance  matrices  must  be  calculated  to 
determine  the  current  distribution  among  the  dipole  elements. 
The  matrix  of  open  circuit  mutual  impedances  between  the 
dipole  elements  is  represented  as 


IZ)  = 


*o.o 

*o.i  ' 

*0,2  •  •  • 

*  *0,  IN- 1) 

*1.0 

*1.1. 

*1.2  •  •  • 

*1.  IN- 1) 

*2.0 

*2.1 

*2,2  *  ‘  • 

*  *2 ,  1  N-l) 

■ 

■ 

.... 

*  1 

. 

. 

.  .  .  . 

•  • 

.*(*-!>  ,o 

*(W- 1)  .  1 

*<W-1),2  ‘  *  V 

•  *(W-1> .  (N-l). 

(15.6) 


Each  main  diagonal  term  (Zu)  is  the  self  -  impedance  of  the  i^+l 
dipole.  The  off-diagonal  terms  (Zjk)  are  the  mutual  impedances 
between  the  i^+l  and  k^+l  elements.  Subscripts  i  and  k  are 
matrix  indices  from  0,1,2. .. .N-l. 

The  matrix  of  open  circuit  mutual  impedances  between  the 
actual  dipoles  and  their  images  is  represented  as 


oH 

o 

**0.1 

**0 ,2  *  • 

•  ’  **l(J\r-l) 

**1.0 

**1.1 

**1.2  •  • 

•  *  **1,(AM) 

[ zi ]  = 

**2.0 

• 

**2.1, 

•  **2,2  '  * 

•  •  • 

•  '  **2 ,  'N-l) 

•  •  1 

• 

**<AM),0 

**(»- 1),1 

•  •  • 

**<W- 1),2  *  • 

•  •  • 

.  .  ZI{N.  1)((lV_i, 

The  main  diagonal  terms  (ZIM)  are  the  mutual  impedances 
between  the  i^+l  element  and  its  image,  and  the  off-diagonal 
terms  (ZIjJt)  are  the  mutual  impedances  between  the  i^+l  element 
and  the  image  of  the  k^+l  element. 
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The  main  diagonal  self  -  impedance  terms  of  equation  15.6 
( Z, , )  are  found  from  [Ref  6:  pp.  205-206] 


Z,  ,  =  - 


60 


[Q{Ua)  -2Q(U1)] 


(15.8) 


1,1  l-cos(2P  hi) 

[Q(V0)  -2Q(V1)]+2[Q(U'0)  -QU\)  -i?(V1)]  +  2!?(t//0)[l+cos(2ph,)]  j 

where  Q(x)  =Ci  ( x )  -jsi  (x)  =  j  9osLY\  dy-j-J  (15.9) 

C/0  =  P  [^2radl+4h|-2hJ  U/0=>/2Pradi  V0  =  P  [/2rad|+4h|+2hij 

C4  =  p  [y/iradj+hj-hj  V^P  [/iTadJ+hJ+Aj] 


For  ^=90°,  the  off-diagonal  mutual  impedance  terms  of 
equation  15.6  (Zitl)  are  given  by  [Ref  11] 


Zi.k~Zk.i: 


60 


(15.10) 


cos  [P  ihrhk)  ]  -cos  [p  (h^h*)  ] 

[e-jP  ( uQ )  -o ( ux )  -0 ( U2 )  ]+e 'hi'hk) [0 ( V0 )  -P ( Vi )  -0 ( V2 )  ] 

+  e -^P (u'o-Q{Vl)-Q(U2)]+ej* 'v‘v [0 ( W0 }  -<? ( )  -0 ( V2 ) ] 
+2Q($Dik)  (cos  [p  +cos  [p  v^+A*)  ] )] 


where  Q(x)  =Ci(x)  -j‘Si(x) 

=  ?  cos  (y)  .  j.f  sin  (y)  . 

“>0 

t 

->« 

£/0  =  p  [lDl^(h^hk)2-(h^hk)\ 

v0=P  [/Oi,k*ihi*hk)i*(hi*hk)\ 

u'0* P  yD!.k+(hrhk)»-{hrhk)] 

^0=P 

^  =  P  [plk+hl-h] 

*i  =  P  [^jc+Ai+Aj 

tf2  =  P  [pl.k+ht-h,] 

V2  =  P  [^ia.*+A^Aj 
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For  ¥  less  than  90°,  the  dipoles  are  in  echelon.  The 
expression  for  the  mutual  impedance  terms  of  equation  15.6 


(Zlk)  for  an  echelon  configuration  is  given  by  [Ref  12] 

Zitj=15{e'j^ihi-SDi-k)[0(UQ) ?(^)]  ,(15.11) 
+  e ihitSD^[o  ( u'0 )  -Q  ( U2 )  ] *  e Jp  (  ^ )  .Q  (  v2 )  ] 


rjt(Hr2hk-SDitk) 


-jp  <hi*2hk+SDiik) 


[L QiUt)-Q{U2)]  +  e 


3jf><hr2hk-SD1'k) 


Jp(/»i*2hfc-SDi,*) 


[0(v3)  -Of  VJ] 
[Q(v,)-q(v2)) 


+2cos  ( Phi)  e'ipsDi*[£)(V2)  ~Q{  W2)  ]+2cos  ( pAj)  ej*SDl\Q(Y2)  -Q{YJ\ 
+2cos(PAi)[e*'P(2VSDj*,[0(W,2)  -£)(W3)]+ei|5t2^SDi  *)p(y2)  -0(r3)]  J  } 

X  X 

where  Q{x)  =Ci  ( x )  -j-(Si  (x)  =J ---s  dy-j dy 

m  y  o  y 

u0  =  P  [v iDlk+(SDi'k-hi)  2+  (SDitk-hi)  ] 

Vq  =  P  [pl.k+  (SDi't-hJ  2-  (SDi't-hi)  ] 

^o-P  [ff^-(®i,^)| 

V'o-P  [v^iV 

A“P  [v^i.*+ S.SD^-h^hj,)  2^SDiik-h^hk) } 

Vi=p  [K*+  VA*> 2-  VV ] 

^2aP  [\^iV  <SUlffc+V**> 2-  (S0i>Jk+VA*>] 

V2-P  [^V  (SDi  k+hi+hk)  2+  {SDik+hi+hk)  ] 
y3  =  P  [^a,*+  (SI^-A^A*)  2  +  (5DUk-hi+2hk)  ] 

V3  =  P  [}jDfi k+  ( SI?,  k-hi+2 hk)  2-  (SD^-A^A*)] 

C/4  =  p  *+  (SO,  k+hi  +  2hk)  2-  ( SDi  k*h^2h} 


V4  =  P  [jDlk+  {SDi  k+hi*2hk)  2  +  (SDi  k+hi  +  2hk) } 

^i  =  P  [\JDlk+SDf'k-SDi  k] 

VP 

^2  =  P  [^E7TiD~^-(5Di<Jc+v] 

^2  =  P  [v/A2,^  (SDiik+hk)  2+  (SDiik+hk) }  ' 

V3  =  p  (SDIiJc+2h*)  2-  (SD^+2^)] 

^3  »  P  [\/  (^i. *♦  < SD, ,  **2/1*)  2  +  ( SDiik+2hk)  ] 

The  Dj-k  terms  in  equations  15.10  and  15.11  are  the  perpen¬ 
dicular  distances  between  the  i^+l  and  k^+l  elements  given  by 


Jc-l 


n»i 


sin  (T) 


(15.12) 


The  SDjk  terms  for  equation  15.11  are  the  vertical  distance 

i 

from  the  center  of  the  i*+l  dipole  to  the  z  coordinate  of  the 
nearest  end  of  the  k^+l  dipole  and  are  given  by 


Jc-l 


SDilk~SDkti^da-  cos (V)-hj 


(15 


13) 


For  ground  planes  with  high  conductivities  relative  to 
the  operating  frequency,  it  is  necessary  to  compute  equation 
15.7,  [ZI] ,  for  the  mutual  impedances  between  the  elements  and 
their  images.  The  main  diagonal  terms  (ZIlt)  are  the  mutual 
impedances  between  the  i^+l  dipole  and  its  own  image.  These 
Zu  terms  are  calculated  from  the  mutual  impedance  equation  for 
a  collinear  configuration  given  by  [Ref  12] 
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+eJP*VMi.i>[0(  v>o)  -Q(v2)]+e-j*(hi'hiii)  ln|  ^'-~~ 


+  e^'.1/*A<I.,>[C(y3) -Q{Ux)  ]  +  e 


-jfiiht-hij'i)  ■ 


hij.j 

hii.i+hi 


[QWJ-QlVjy 


-riOh^hii 


ejp(3Ai*hii,i)^^y.j  ~q{v2)  j+e';P'3hi*'hIi'i)ln  h^£— 

+2cos  ( phj)  e^hii,i[Q[Y2)  -£>(  1^) ]+2cos  ( (JAj)  e'J?Aii,1ln  - 
+2cos  ( p/jj)  ej*(2hrl>iii)[Q(Y2)  -Q(  r3)] 
*2ccS,pAlle^-.-,,lln[^^i]j 

where  OU)  -Ci  <*>  -j-Si  <*>  .f£S^Xldy-j-U Hi 
Uo-P  toi.i,*  (ilij.j-Aj)  2*  (Wi.c-Aj>  ] 


hi+2hk)  z+  (Ai 


hi+2hk) 2-  (Ai 


t/4  =  p  yDlk+  {hi^+h^hj 2-  [hiik+hi+2hk)  ] 

V4  =  P  yDl^(hilik  +  hi  +  2hk)  2+  (hiiik+hs2hk)  ] 

*1=P  [pl^hU't-hi,' J  vP 

V2  =  p 

*2«P  [v,£>i.*+  (hiik+hk)  2+  {hiiik+hk)  ] 

V3  =  p  yDf'k+  ( hilik*2hk ) 2-  ( hiiik+2hk )  ] 
i'3  =  P  [/ (D-ik*  (hii§k*2hk)  2+  (Aiijit+2AJt)] 

The  symbol  hi(j  in  equation  15.14  is  the  height  of  the  im+l 
dipole  above  the  nearest  end  of  its  own  image,  and  with  i=k  is 
obtained  from 

(15.15) 

The  off-diagonal  mutual  impedances  between  the  dipoles  and  the 
images  (ZIa)  are  found  using  equation  15.13,  except  the  Q(x) 
arguments  are  those  given  for  equation  15.14,  but  with  SD,k 
replaced  by  hilvk  from  equation  15.15. 

The  function  Q(x)  given  by  equation  15.9  is  defined  in 
terms  of  sine  [Si (x) ]  and  cosine  [Ci (x) ]  integrals.  Mathcad 
is  incapable  of  evaluating  Si(x)  and  Ci(x)  directly,  so  a 
series  expansion  and  polynomial  approximation  are  written  into 
the  Mathcad  code  for  both  the  sine  and  cosine  integrals  to 
evaluate  the  functions  Q(x).  The  arguments  for  Q(x)  are 
always  real,  so  it  is  unnecessary  to  find  Si(x)  and  Ci  (x) 
expressions  valid  for  complex  numbers.  For  arguments  less 
than  one,  Si(x)  and  Ci(x)  are  evaluated  by  series  expansions 
given  by  rRef  9:  p.  232] 
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In  addition  to  the  mutual  impedance  matrices,  it  is 
necessary  to  calculate  the  short  circuit  admittance  matrix  to 
determine  the  elemental  current  distribution.  The  short 
circuit  admittance  matrix  is  given  by  [Ref  10:  p.  321] 


^2,1  ^2 , 2  *2,3 


[Y]  = 


ri'o.a  *Vi  0  0 

*1,0  *1,1  *1.2  0 

0 
0 

0 
0 


0 

0 

0 


*3,2  *3,3  *3 , 4 


0 

0 


0 

0 


0 

0 


*(W-2)  ,  (W-3)  *(W-2),ftf-2)  *fW-2),(W-l) 

0  *<//-!)  ,  (AT-2)  *<JV-1)  (W-l) 


[15.22) 


where  the  diagonal  terms  are  given  by 
*o (o  =  *"J *o cot  ( P c?q )  and  Y^,^  T-jYQcot(^dN.1) 

and  for  0  <  n  <  N- 1 : 

*n,n*  •  •  •  *(W-  2)  ,  (N-2)~~j^0  (c°t  (  P^n-i)  +COS  (  P<d„) 

and  the  off  diagonal  terms  are  given  by 

*<«-!),  n=*n.  (n-1)  —  J  *qCSC  ( 0da.i) 

The  term  Y0  is  the  characteristic  admittance  of  the  trans¬ 
mission  line  given  by  the  input  (ADM)  ,  and  Y'T  is  given  by 


W  -y 
z  T  1 0 


COS 

P-^i 

k _ ? _ 

+jY0Z7sin|p^| 

YgZjjeosfp  - 

Vi] 

2  J 

+jsin 

(15.23) 


where  Zx  is  the  termination  impedance  input  (TIMP) . 

The  ground  conductivity  relative  to  the  operating 
frequency  (<rrel)  is  then  computed  by  [Ref  6:  p.  640] 


°rel  = 


o>ere0 


18000  o 

^Mhz  er 


(15.24) 
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to  determine  which  equation  to  use  for  the  array's  current 
distribution.  A  surface  with  a  relative  conductivity  greater 
than  20  is  considered  to  be  highly  conductive  relative  to  the 
transmitted  frequency  and  the  current  distribution  is  given  by 
[Ref  10:  p.  342] 

[IB]  ={  CC/]  +  [Y]  [Z]+m  [ZI]  }"*  [IREF]  (15.25) 

where  [U]  is  the  N  X  N  identity  matrix.  The  entries  of  the 
[IB]  matrix  are  the  base  currents  of  the  dipole  elements.  The 
matrix  [IREF]  represents  the  input  current  to  which  the  [IB] 
distribution  is  referenced.  Since  the  current  input  is  only 
to  the  feed  at  the  center  of  the  first  element,  and  since  the 
input  current  is  sinusbidal  with  a  maximum  of  unity,  [IREF]  is 
a  1  X  N  matrix  with  a  one  as  the  first  entry  and  N-l  zeroes  in 
the  remaining  positions. 

For  surfaces  with  a  relative  conductivity  of  less  than 

! 

I 

20,  the  mutual  impedances  due  to  the  images  are  neglected  and 
the  current  distribution  is  given  by 

[I£]={[U]  +[y]  [Z]  I'M  IREF]  (15.26) 


The  electric  field  equations  are  referenced  to  the 
current  maxima  values  for  the  dipole  elements,  [I] ,  which  are 
obtained  from  the  [IB]  matrix  using  [Ref  6:  p.  209] 


sin(Phi) 


(15.27) 
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The  mutual  impedance  calculations  are  not  valid  when  any 
element  length  is  an  exact  integer  multiple  of  the  wavelength. 
When  this  occurs,  there  is  a  singularity  error  in  the  mutual 
impedance  calculations,  and  it  is  necessary  to  vary  the 
frequency  such  that  no  element  is  an  exact  integer  multiple  of 
the  wavelength.  The  change  required  is  only  one  or  two 
percent,  and  the  predicted  radiation  parameters  are  still  a 
good  estimate  of  tnose  for  the  original  frequency. 

The  electric  field  for  the  vertical  log-periodic  array  is 
obtained  for  each  individual  dipole  in  a  manner  analogous  to 
that  used  for  the  vertical  dipole  discussed  in  Chapter  4.  The 
equations  for  the  individual  elements  are  combined  into  a 
single  expression  for  the  array's  total  electric  field  by  the 
array  factor.  The  array  factor  for  the  vertical  log-periodic 
array  is  [Ref  6  p.  222] 

AM. 

A+jB=£  Iiej*H'cosi9)[cos  [p^cos  (6)  ]  -cos  (ph,)]  •  (15.28) 

i=Q 

1  +  e-^pw'cos(e)|rv+(l-rv)Fff  |sin2  (0)  _^Bpmcos(e)))  • 

e  JPlTjSoc  («,♦«,)  [cos  (8)  sin(aj»B5)  ♦sin (8)  cos  (a2*«3)  sin  ($)  ] 


where  Y,  is  the  y  coordinate  of  the  i^+1  dipole  given  by 


rtm 


\n-0  ) 


sin  (7) 


(15.29) 
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The  first  two  terms  inside  the  brackets  of  equation  15.28 
represent  the  space  wave,  the  third  term  represents  the 
surface  wave,  and  I,  is  the  current  term  for  each  element  based 
on  a  sinusoidal  current  input  with  a  maximum  of  unity.  The 
array  factor  is  combined  with  the  element  factor  to  find  the 
expression  for  the  total  radiated  electric  field  given  by 


je  Oe-'9*  A+iB 
0  R  sin (6) 


(15.30) 


The  requested  inputs  are  used  to  calculate  the  following 
variables  using  a  constant  0  of  $=0  and  0=7t/2  for  312  discrete 
values  of  6  which  are  equally  incremented  from  -ir/2  to  tt/2  : 

h, . half-length  of  the  i^+l  dipole  (meters) 

H,  . height  of  the  i4+l  dipole  (meters) 

Y . .  coordinate  of  the  i^+l  dipole 

I,  . . current  for  the  ilh+l  dipole 

X5 . wavelength  of  the  operational  frequency  (meters) 

(3.  ......  free  space  wavenumber  for  operational  frequency 

n . index  of  refraction 

Pt . complex  numerical  distance  for  vertical  polarization 

r\ . vertical  reflection  coefficient 

F, . vertical  surface  wave  attenuation  factor 

The  calculated  variables  are  used  to  evaluate  the  far- 
field  space  wave  and  surface  wave  for  the  discrete  values  of 
9  with  0=0  and  0=jt/2.  The  space  wave  and  surface  wave  results 
are  combined  for  corresponding  values  of  9  to  obtain  the  total 
radiated  electric  field  distribution  for  the  0=0  and  0=7t/2 
vertical  planes.  The  space  wave,  surface  wave,  and  total 
electric  field  results  are  then  normalized  with  respect  to  the 
maximum  field  intensity  of  each,  and  the  normalized  magnitudes 
are  plotted  for  each  discrete  9  to  depict  the  radiation 
patterns.  The  Mathcad  vertical  log-periodic  dipole  array 


117 


( 


application  computes  the  space  wave,  surface  wave,  and  total 
electric  field  radiation  patterns  and  radiation  parameters  in 
the  0=0  and  <p=ir/2  vertical  planes. 

The  variables  corresponding  to  the  selected  elevation 
angle  index  (w)  are  used  to  evaluate  the  radiated  electric 
field  components  for  the  space  wave  and  surface  wave  at  the 
fixed  elevation  angle  (9W)  as  0  varies  from  0  to  2n  in  312 
equal  increments.  The  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  for  the  vertical  planes. 

Equations  3.8  and  3.9  are  used  to  integrate  equation 

15.30  over  the  hemispherical  Gaussian  surface  above  the  grcand 

plane  at  a  fixed  radius  (R)  from  the  array  to  find  total 

average  radiated  power  (Prad)  .  With  the  discrete  values  of  the 

electric  field  and  total  average  radiated  power  determined, 

the  Mathcad  application  predicts  the  following  radiation 

characteristics  from  the  equations  in  Chapter  3. 

Rrad  . radiation  resistance  (Ohms) 

D0 . directivity 

EIRP . effective  radiated  isotropic  power  (Watts) 

A^  . maximum  theoretical  effective  area  (square  meters) 

l^u,  . maximum  theoretical  effective  length  (meters) 

P . numerical  distance  (vertical  polarization,  0=90°) 

Pm  . numerical  distance  (horizontal  polarization,  0=90°) 

Angles  .elevation  angle  of  maximum  directive  gain  (degrees) 

The  directivity  (D0) ,  effective  isotropic  radiated  power 

(EIRP) ,  maximum  theoretical  effective  area  (A^)  ,  maximum 

theoretical  effective  length  (1^)  ,  and  the  elevation  angle  of 

maximum  directive  gain  (Angle^)  are  all  determined  for  both 

the  0=0  and  <p=n/2  vertical  planes. 


V 
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As  an  example,  the  Mathcad  vertical  log-periodic  array- 
application  was  executed  with  the  following  inputs: 


number  of  elements 

length  of  first  dipole 

length  of  second  dipole 

first  separation  distance 

radius  of  first  element  0. 

height  of  antenna  feed 

frequency 

distance  from  the  antenna 
relative  dielectric  constant 
ground  conductivity- 
characteristic  admittance 
termination  impedance 
angle  from  vertical  to  array  axis 
elevation  angle  index 


12 

3.246  meters 
3.732  meters 
1.096  meters 
00325  meters 
8.0  meters 
10  -106  Hertz 
3000  meters 
4 

1-10‘J 

1/450  Mhos 
0  Ohms 
90° 

285  (-17°) 


Figures  15.2  through  15.7  are  the  space  wave  and  surface  wave 
radiation  patterns  in  the  <£«0  and  4)=*tt/2  vertical  planes,  and 
the  designated  horizontal  plane  for  this  example. 


plane. 


vertical  plane. 
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The  following  radiation  parameters  were  predicted  by  the 
Mathcad  application  for  a  sinusoidal  current  input  of  one  Amp. 


Total  power  radiated  (Watts) 
Radiation  resistance  (Ohms) 
Numerical  distance  (vertical) 

37.967 

75.934 

57.126 

6=it/2 

Directivity 

3.415 

12.176 

EIRP  (Watts) 

129.656 

462.272 

Max  eff  area  (sq  meters) 

244.578 

872.012 

Max  eff  length  (meters) 

14.038 

26.506 

Anglem„  (degrees) 

21.71 

21.14 

These  results  are  consistent  with  expectations  for  this 
particular  configuration.  Appendix  D  contains  computer 
hardcopies  of  additional  example  calculations  for  the  vertical 
log-periodic  dipole  array  and  compares  predicted  radiation 
parameters  to  those  expected  based  on  previous  calculations. 
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XVI.  THE  HORIZONTALLY  POLARIZED  LOG- PERIODIC  DIPOLE  ARRAY 


The  orientation  of  the  horizontally  polarized  log- 
periodic  dipole  array  is  depicted  in  Figure  16.1,  where  N  is 
the  number  of  dipoles,  n  is  an  index  equal  to  0,1,2..N-1,  ln 
is  the  length  of  the  n^+l  dipole,  d„  is  the  distance  between 
the  n^+l  and  ^+2  dipole  elements,  H0  is  the  height  of  the 
antenna  feed  above  ground,  H„  is  the  height  of  the  n^+l  dipole, 
^  is  the  angle  between  the  vertical  axis  and  the  array's 
center  axis,  a  is  the  angle  between  the  array's  center  axis 
and  the  lines  which  connect  the  tips  of  the  elements,  R  is  the 
radial  coordinate,  6  is  the  elevation  coordinate,  and  4>  is  the 
azimuth  coordinate. 

The  Mathcad  application  for  the  horizontally  polarized 
log-periodic  array  requires  the  following  inputs: 


N . number  of  dipole  elements 

10 . length  of  the  shortest  dipole  (meters) 

1,  . length  of  the  second  shortest  dipole  (meters) 

do . distance  between  first  two  elements  (meters) 

rado . radius  of  the  shortest  dipole  (meters) 

Hq . feed  height  above  ground  plane  (meters) 

¥ . angle  from  vertical  axis  to  array  axis  (degrees) 

f3 . operational  frequency  (Hertz) 

R . distance  from  array  (meters) 

er . relative  dielectric  constant  of  ground  plane 

a . conductivity  of  ground  plane 

w . elevation  angle  index  (from  Table  3.2) 

ADM . transmission  line  characteristic  admittance 

TIMP . termination  impedance  connected  to  longest  dipole 


The  user  inputs  the  frequency  for  which  the  radiation 
parameters  discussed  in  Chapter  3  are  computed.  The  lengths 


of  the  first  two  elements  (10  and  lt)  are  used  to  calculate  the 
log-periodic  relationship  (r)  given  by  [Ref  10:  p.  317] 

t  =  (n=0, 1, 2 _ N-2.N-1)  (16.1) 

■‘■a*  l 


The  remaining  element  lengths  (12  to  1N.,)  are  then  found  from 


(16.2) 


As  noted  in  equation  16.1,  indices  start  at  zero  instead  of  1, 
so  the  index  simply  refers  to  the  index- plus -one  position  in 
terms  of  successive  values  for  a  given  parameter.  This  is  to 
maintain  continuity  with  Mathcad  which  also  indexes  from  zero. 
The  dipole  half-lengths  (hj  are  found  from  ln/2 .  The  radius 
of  the  first  element  (radj  is  used  to  find  remaining  radii  by 
substituting  rad„  and  rad,,*,  for  ln  and  ln+1  into  equation  16.2. 
The  angle  between  the  vertical  axis  and  array  axis  (¥)  can  be 
between  0°  and  90°,  and  a  is  found  from  [Ref  10:  p.  319] 


— 2-=(l-t)cot(o)  (16.3) 

1 

which  can  be  written  as 

g  =  tan-1[(1"Cj-hntl1  (16.4) 

dn 

The  remaining  separations  (d,  to  dN.2)  are  found  from  equation 
16.3.  Feed  height  (H*,)  can  be  greater  than  or  equal  to  zero. 
Distance  from  the  array  (R)  must  meet  the  far- field  require¬ 
ments  of  Chapter  3.C.  The  elevation  angle  index  (w)  sets  the 
6  coordinate  for  which  a  horizontal  radiation  pattern  is 
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determined.  Table  3.2  lists  possible  indices  and  their  cor¬ 
responding  elevation  angles  from  .57°  to  89.1°  in  increments 
of  about  2.3°.  Indices  between  those  listed  can  be  used  to 
interpolate  a  better  estimate  of  a  desired  elevation.  The 
termination  impedance  (TIMP)  is  assumed  to  be  connected  to  the 
array's  center  axis  opposite  the  feed  at  a  distance  of  hN.,/2 
from  the  last  element.  Typical  termination  impedances  are  the 
transmission  line's  matched  impedance  (complex  conjugate)  or 
a  short  circuit  (zero) .  The  characteristic  admittance  (ADM) 
is  just  the  inverse  of  the  transmission  line's  characteristic 
impedance.  If  computed  lengths,  separations  and  radii  are  not 
desired,  they  can  be  entered  manually  on  page  two  of  the 
application.  Consult  a  Mathcad  manual  to  be  certain  the 
entries  are  made  correctly. 


The 

open 

circuit 

impedance 

matrices  are 

calculated  to 

determine 

the  current  - 

distribution  among  the  dipole  elements. 

The  matrix  of 

mutual 

impedances 

between  dipole  elements  is 

represented  as 

Zo.o 

*0.1 

*0,2  •  -  * 

•  •  Z0,  (W-l) 

. . . . . 

*1.0 

*i.i 

*1.2  *  • 

•  •  *1,  (W-l) 

t  Z]  = 

• 

• 

•  •  • 

•  •  • 

(16.5) 

i.*(W-l),0  *<iM),  1 

*W-1),2  *  • 

•  •  • 

•  •  *<N-1)  ,  (W-l). 

Each  main  diagonal  term  (Zu)  is  the  self-  impedance  of  the 
i^+l  dipole.  The  off-diagonal  terms  (Zik)  are  the  mutual 
impedances  between  the  i^+l  and  k*+l  elements.  Subscripts  i 
and  k  are  matrix  indices  (0,1,2 _ N-l) . 
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The  matrix  of  open  circuit  mutual  impedances  between  the 
actual  dipoles  and  their  images  is  represented  as 


[ZI]  = 


ZJ0i0 

ZI0.i 

^0.2  • 

.  .  .  ZIq,  (W-l) 

ZI  1.0 

^1.2  * 

•  • 

.  .  . 

«I1MI  .0 

ZI(N-l) ,  1 

•  • 

^"<W-1),2  • 

•  •  *  'II(N-l)  .  (AT-1). 

(16.6) 


The  main  diagonal  terms  (ZI,,)  are  the  mutuaj.  impedances 
between  the  i^+l  element  and  its  image,  and  the  off-diagonal 
terms  (ZIiik)  are  the  mutual  impedances  between  the  i^+l  element 
and  the  image  of  the  klh+l  element. 

The  main  diagonal  self  - impedance  terms  of  equation  16.5 
are  found  from  [Ref  6:  pp.  205-206] 


60 

i,i  1-cos  (2phi) 


[<?(£/«,)  -2Q(Uj)] 


(i6.7: 


+  ej2*h‘  [(?(  V0)  '-2{?(  Vj)  J+2[C?(f//0)  -Q(UX)  -j?(  V^)  ]+2j?(  C/^)  [1+cos  (2phi)  ]  ] 

X  X 

where  Q(x)  =Ci  ( x )  -j-Si  (x)  =  f  co3/>r)  dy-j-f  sin^-  dy  (16.8) 

i  y  o  y 

U0= P  [j2rad}+4hj-2h^  C7/0=v/2Pradi  V0=p  [^2 rad] + 4 hf+2h^ 

U^V^ftxadj+hj-h^  VP  [j2radi*hj*h1] 

I 

The  of ^-diagonal  mutual  impedance  terms  of  equation  16.5 
are  given  by  [Ref  11] 

60 


*i.*“3*> 


(16.9) 


cos  [p  (bi-h*)  ]  -cos  [p  {ht+hk)  ] 

[»  -jp  ( V**  p  (  Uq  )  _  Q  { q )  _  q  ( ^  j  j+ eJt  (*i  [Q(Vq)  -Q(VX)  -Q(V2)] 

+  e~i*'h™[Q{U,9-Q[VJ)-Q{Ui)yei*{bi'h*]p{V,0)  -0(C/x)-i?(y2)] 
+2C’(P^,JC)(cos  [p  (hx-hk)  ]  +cos  [p  (hi+hk)  ] )] 
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X  x 

where  Q(x)  -Ci  ( x )  -j-Si  (x)  =  J SBllZLdy  -  j-f 


sin(vO 


dy 


tf0=p  [fin  +(Ai+hJC)2-(v*,c>] 

v0=p  (h^h 

t/'o-p  yDik+(hrhk)>~(hrhk)] 

^0=p  yDik*(hrhk)^{hr 

U.=V  [/Dl.t+hi-hi] 

^=p  yDit+hf+hi] 

[M  »♦*/-*»]' 

^2=P  [fi>Fk+hk+hk] 

The  Djk  terms  are  the  distances  between  the  i^+l  and  k^+l 
dipole  elements  and  are  given  by 


(16.10) 

n-i 

For  ground  planes  with  high  conductivities  relative  to 
the  operating  frequency,  it  is  necessary  to  compute  equation 
16.6,  [ZI] ,  for  the  mutual  impedances  between  the  elements  and 
their  images.  The  mutual  impedance  terms  (ZIi-k)  are  given  by 
equation  16.9  with  the  Q(x)  arguments 


u0= P  ydilk+(hi+hk)2-{hi+hk)\ 

v0  =  p  [ycriJ2.Jt+(A1+/iJt)*+(Ai+AJf)j 

u'0= P  [^/diJ2.^(hi-V2-(V^)] 

vVP  [^/dii2,Jc+(hi-ii*)2-*-(hi-hi)] 

^=P  [Idil^hj-h,] 

^i=P  [^i.*+*i2+*J 

C/2=p  [v/cfii.  jc+ 

V2=P  [jdilk+hZ+hk] 

The  dii-k  entries  are  the  distance  from  the  i*h+l  dipole  to  the 
image  of  the  k^+l  dipole  and  are  given  by 


2+  (^i,jcsin(T) ) 2 


(16.11) 
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The  function  Q(x)  given  by  equation  16.8  is  defined  in 
terms  of  sine  [Si(x}]  and  cosine  [Ci(x)]  integrals.  Mathcad 
is  incapable  of  evaluating  Si (x)  and  Ci (x)  directly,  so  a 
series  expansion  and  polynomial  approximation  are  written  into 
the  Mathcad  code  for  both  the  sine  and  cosine  integrals  to 
evaluate  the  functions  Q{x).  The  arguments  for  Q(x)  are 
always  real,  so  it  is  unnecessary  to  find  Si(x)  and  Ci(x) 
expressions  valid  for  complex  numbers. 

For  arguments  greater  than  one,  two  auxiliary  functions, 
f(x)  and  g(x)  are  determined  by  [Ref  9:  p.  233] 


x8+a1x6+a2x4+a3x2+a4 ' 
xB+b1x6  +h2x4  +b2x2 + b4 1 


+e(x) 


where 

ax=38. 027264 
a2=265. 187033 
a3=335. 677320 
a.  =38 . 102495 


je  (x)  |  <  5-10*7 

b3=40. 021433 
b2=322. 624911 
jb3  =  57  0.236280 
b,  =157 . 105423 


g(x) 


x° + a,  x” + a,x  *  +  a,x  * + a 


y  x8+bjX6+b2x4+b3x2+b4 


±  +e(x> 


where 

a3=42. 242855 
a2=302. 757865 
a3=352. 018498 
a.=21. 821899 


|e(x)  |  <  3-10‘7 

b2  =  48  . 196927 
b2=482 .485984 
b3=1114. 978885 
b.=449 .690326 


(16.12) 


(16.13) 


\ 

4 


/ 
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and  Si(x)  and  Ci{x)  are  defined  in  terms  of  these  auxiliary 
functions  by  [Ref  9:  p.  232] 


Si  (x)  =-^  -f  (x)  cos  (x)  ~g(x)  sin  (x) 
Ci  (x)  =f(x)  sin  (x)  -g(x)  cos  ix) 


(16 . 14) 


(16.15) 


For  arguments  less  than  one,  Si(x)  and  Ci(x)  are 
evaluated  by  series  expansions  given  by  [Ref  9:  p.  232] 


where  y=0 . 5772156649 


(16.16) 


(16.17) 


It  is  also  necessary'  to  calculate  the  short  circuit 
admittance  matrix  to  find  the  elemental  current  distribution. 
The  admittance  matrix  is  given  by  (Ref  10:  p.  321] 


*1.0  *1,1 


[y]=  0 


*o,i  0  0 

*1.1  ^1.2  0 
*2.1  *^2, 2  *2.3 


**3 , 2  *3.3  ^3.4 


(16 . 18) 


*W 


Ur- 3)  1  (W-2)  ,  (W-2)  1  (N-2)  ,  (W-l) 


'  (W-l) ,  {«- 2) 


(AT-1)  (W-l) 


where  the  diagonal  terms  are  given  by 
*o,o  =  ~.7  *oco*-  ( and  *^/^i) ,  tw-i)  =  *r/r”-7  *oc°t  ( Pdw_ 
and  for  0  <  n  <  N- 1: 

i'n.n - *V2> .  (W-2)  —  J  ^  (cot  (  Pdn.i)  +<=OS  (  Mn) 

and  the  off  diagonal  terms  are  given  by 
*01-1)  .n=Yn.  lB-l)u-JYoC8cWda_l) 
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The  term  Y0  is  the  transmission  line's  characteristic  admit¬ 
tance  obtained  from  the  input  (ADM),  and  Y'T  is  given  by 


Y1  -r- 


cos 

(  h 
p  n 

^ _ 2_J 

*jY0z,sin 

(  h 

p  T1 

k  2 

Y0ZjCOS 

(f>- 

Vxl 

2  J 

+  jsin 

l"  rj 

(16 . 19) 


The  ground  conductivity  relative  to  the  operating 
frequency  (anl)  is  then  computed  by  [Ref  G:  p.  640] 


’tel' 


coere0 


18000  o 
f Mhz  er 


(16.20) 


A  surface  with  a  relative  conductivity  greater  than  20  is 
considered  to  be  highly  conductive  relative  to  the  frequent/ 
and  the  current  distribution  is  given  by  [Ref  10:  p.  342] 

[IB]  «{  [U]  +  [Y]  [Z]  +  [Y]  [ZI]  }-*  [IREF]  (16 .21) 

where  [U]  is  the  N  X  N  identity  matrix.  The  entries  of  the 
[IB]  matrix  are  the  base  currents  of  the  dipole  elements.  The 
matrix  [IREF]  represents  the  input  current  to  which  the  [IB] 
distribution  is  referenced.  Since  the  current  is  input  only 
to  the  feed  at  the  center  of  the  first  element,  and  since  the 
input  current  is  sinusoidal  with  a  maximum  of  unity,  [IREF]  is 
a  1  X  N  matrix  with  a  one  as  the  first  entry  and  N-l  zeroes  in 
the  remaining  positions. 

For  surfaces  with  a  relative  conductivity  of  less  than 
20,  the  mutual  impedances  due  to  the  image  elements  are 
neglected  and  the  current  distribution  is  given  by 

[IB]»{[TJ]  +  [Y]  [Z]}-‘[IREF]  (16.22) 
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The  electric  field  equations  are  referenced  to  the 
current  maxima  values  for  the  dipole  elements,  [I] ,  which  are 
obtained  from  the  [IB]  matrix  using  [Ref  6  p.209] 


■IB. 


sin  i ph.) 


(16.23) 


The  mutual  impedance  calculations  are  not  valid  when  any 
element  length  is  an  exact  integer  multiple  of  the  wavelength. 
When  this  occurs,  there  is  a  singularity  error  in  the  mutual 
impedance  calculations,  and  it  is  necessary  to  vary  the 
frequency  such  that  no  element  is  exactly  an  integer  multiple 
of  the  wavelength.  The  change  required  is  only  a  percent  or 
two,  and  the  predicted  radiation  parameters  are  still  a  good 
estimate  of  those  for  the  original  frequency. 

The  electric  field  for  the  horizontal  log-periodic  array 
is  obtained  for  each  individual  dipole  in  a  manner  analogous 
to  that  for  the  vertical  dipole  discussed  in  Chapter  4.  The 
equations  for  the  individual  elements  are  combined  into  a 
single  expression  for  the  array's  total  electric  field  by  the 
array  factors.  The  array  factors  for  the  horizontal  log- 
periodic  array  are  [Ref  6:  p.  214] 


•V- 1 


IVcosfPh^cos  (0)  cos  (7)  +sin  (0)  sin  (7)  sin  (cj>)  )]-cos  ( P hi)  ] 

(16.24) 

|sin2  (0)  -  (flTeos  (Qj  j  gjiY,csctW)  !cos<8)cos(T)  *sin(8)  sin(T) 


e-ji0w1c3si8i  +  j  1_p  j  p  c-j2iH1coB(Q)  y(n2-sin2  (0) 
v  v  *  n2cos(0) 


’ ) sin ] 
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N-l 

S^^jT  Ii[cos[Phi(cos  (0)  cos  (7)  +sin  (0)  sin  (7)  sin  (4>)  )]-cos  ( p/2^ )  ] 

i  =  0 

[l+rAe‘-7iPWiC0S<e)  +  (l-rh)  Fm  e -j^P^iCos  <0)  j.  (16.25) 

Jp  l'jCsc(T)  [cos  (8) cos  (?)  *sin  (8)  sin  (?)  sin ($)  1 

where  Y,  is  the  y  coordinate  of  the  i^+l  dipole  given  by 


(i- 1 

£*. 

\n»0  / 


sin  (7) 


(16.26) 


The  first  two  terms  inside  the  brackets  of  equations 
16.24  and  16.25  represent  the  space  wave,  the  third  term 
represents  the  surface  wave,  and  I;  is  the  current  term  for 
each  element  based  on  a  sinusoidal  current  input  with  a 
maximum  of  unity.  The  array  factors  are  combined  with  the 
element  factors  to  find  the  total  radiated  electric  field 
equation,  and  the  expression  for  the  radiated  electric  field 
distribution  of  the  horizontal  log-periodic  dipole  array  is 
given  by  [Ref  6:  p.  213] 


E  =  ~J'60e~jPJ?  cos  (<}>)  CCS  (0)  s 
8  R  1-sin2  (0)  cos2  (<}>)  9 

E  -  J60e"JtR _ sin(4>) _  s 

*  R  1-sin2  (0) cos2 (0)  * 


(16.27) 

(16.28) 


The  requested  inputs  are  used  to  calculate  the  following 
variables  using  a  constant  0  of  0=0  and  0=tt/2  for  312  discrete 
values  of  6  which  are  equally  incremented  from  -ir/2  to  rr/2 : 


v 
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hj . half-length  of  the  i^+i  dipole  (meters) 

Hj . height  of  the  i^+i  dipole  (meters) 

Y, . y  coordinate  of  the  i^+l  dipole 

I, . current  for  the  i^+l  dipole 

X3 . wavelength  of  the  operational  frequency  (meters) 

/3 . free  space  wavenumber  for  operational  frequency 

n . index  of  refraction 

Pe  ......  complex  numerical  distance  (vertical  polarization) 

Pm  . complex  numerical  distance  (horizontal  polarization) 

rv  ..... .vertical  reflection  coefficient 

rh . horizontal  reflection  coefficient 

Fe . vertical  surface  wave  attenuation  factor 

Fm  . horizontal  surface  wave  attenuation  factor 


The 

calculated 

variables 

are 

used  to 

evaluate 

the 

radiated 

far- field 

space  wave 

and 

surface 

wave  for 

the 

discrete  values  of  6  with  0=0  and  0=7r/ 2.  The  space  wave  and 


surface  wave  results  are  combined  for  corresponding  values  of 
6  to  obtain  the  total  radiated  electric  field  distribution  for 
the  0=0  and  0=7r/2  vertical  planes.  The  space  wave,  surface 
wave,  and  total  electric  field  results  are  then  normalized 
with  respect  to  the  maximum  field  intensity  of  each,  and  the 
normalized  magnitudes  are  plotted  for  each  discrete  0  to 
depict  the  radiation  patterns.  The  Mathcad  horizontal  log- 
periodic  dipole  array  application  computes  the  space  wave, 
surface  wave,  and  total  electric  field  radiation  patterns  and 
radiation  parameters  in  the  0=0  and  0=7r/2  vertical  planes. 

The  variables  corresponding  to  the  selected  elevation 
angle  index  (w)  are  used  to  evaluate  the  radiated  electric 
field  components  for  the  space  wave  and  surface  wave  at  the 
fixed  elevation  angle  (0W)  as  0  varies  from  0  to  2n  in  312 
equal  increments.  The  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  for  the  vertical  planes. 


Equations  3.8  and  3.9  are  used  to  integrate  equations 

16.27  and  16.28  over  the  hemispherical  Gaussian  surface  above 

the  ground  plana  at  a  fixed  radius  (R)  from  the  array  to  find 

total  average  radiated  power  (Pfad)  .  With  the  discrete  values 

of  the  electric  field  and  total  average  radiated  power 

determined,  the  Mathcad  application  predicts  the  following 

radiation  characteristics  from  the  equations  in  Chapter  3: 

Rrad  . radiation  resistance  (Ohms) 

D0 . directivity 

EIRP . effective  radiated  isotropic  power  (Watts) 

A^,  . maximum  theoretical  effective  area  (square  meters) 

1^,  .....maximum  theoretical  effective  length  (meters) 

P . numerical  distance  (vertical  polarization,  0  = 90°) 

Pm  . numerical  distance  (horizontal  polarization,  0=90°) 

Angle^,  .ele  *ation  angle  of  maximum  directive  gain  (degrees) 

The  directivity  (D0) ,  effective  isotropic  radiated  power 

(EIRP) ,  maximum  theoretical  effective  area  (A^)  ,  maximum 

theoretical  effective  length  (l^J ,  and  elevation  angle  of  the 

directivity  (Angle^J  are  determined  for  both  vertical  planes. 

As  an  example,  the  Mathcad  horizontal  log -periodic  array 
application  was  executed  with  the  following  inputs: 


number  of  elements 
length  of  first  dipole 
length  of  second  dipole 
first  separation  distance 
radius  of  first  element 
height  of  antenna  feed 
frequency 

distance  from  the  antenna 
relative  dielectric  constant 
ground  conductivity 
characteristic  admittance 
termination  impedance 
angle  from  vertical  to  array  axis 
elevation  angle  index 


12 

3.246  meters 
3.732  meters 
1.096  meters 
0.00325  meters 
8.0  meters 
20 -106  Hertz 
3000  meters 
4 

5  •  10'3 
1/450 

0  Ohms 
90° 

285  (-17°) 
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XVII.  THE  HORIZONTAL  YAGI-UDA  ARRAY 

The  orientation  of  the  horizontal  Yagi-Uda  array  is 
depicted  in  Figure  17.1,  where  1NS  is  the  length  of  the 
current  driven  element,  NR  is  the  number  of  parasitic 
reflector  elements,  ND  is  the  number  of  parasitic  director 
elements,  n  is  an  index  equal  to  0, 1, 2 . .NR+ND,  and  ln  is  the 
length  of  the  n^+l  element  starting  with  the  outermost 
reflector  and  counting  consecutively  toward  the  opposite  end. 
The  terms  are  the  separation  distances  between  the  n*+l  and 
n®+2  elements,  Hq  is  the  height  of  all  antenna  elements  above 
ground,  R  is  the  radial  coordinate,  6  is  the  elevation 
coordinate,  and  </>  is  the  azimuth  coordinate. 

The  Mathcad  application  for  the  horizontal  Yagi-Uda  array 

requires  the  following  inputs: 

NR . number  of  reflector  elements 

ND . number  of  director  elements 

ln . length  of  the  n^+l  element  (meters) 

dg . distance  between  n^+l  and  n®+2  elements  (meters) 

rad„  .....  radius  of  the  n^+l  element  (meters) 

Hq  .... ...antenna  height  above  ground  plane  (meters) 

fj . operational  frequency  (Hertz) 

R . distance  from  array  (meters) 

e , . relative  dielectric  constant  of  ground  plane 

a . conductivity  of  ground  plane 

w . elevation  angle  index  (from  Table  3.2) 

Unlike  the  log-periodic  applications,  all  element 
lengths,  radii,  and  separation  distances  must  be  input 
manually  on  page  two  of  the  application.  Consult  a  Mathcad 
manual  to  be  certain  the  entries  are  made  correctly.  The  user 
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inputs  the  frequency  for  which  the  radiation  parameters 
discussed  in  Chapter  3  are  calculated.  As  noted  in  the  first 
paragraph,  indices  start  at  zero  instead  of  1,  so  the  index 
simply  refers  to  the  index- plus -one  position  in  terms  of 
successive  values  for  a  given  parameter.  This  is  to  maintain 
continuity  with  Mathcad  which  also  indexes  from  zero.  The 
dipole  half-lengths  (hj  are  found  from  ln/2.  Antenna  height 
(Hq)  can  be  any  value  greater  than  or  equal  to  zero,  and  the 
distance  from  the  array  (R)  must  meet  the  far- field 
requirements  of  Chapter  3.C.  The  elevation  angle  index  (w) 
sets  the  6  coordinate  for  which  a  horizontal  radiation  pattern 
is  determined.  Table  3.2  lists  possible  indices  and  their 
corresponding  elevation  angles  from  .57°  to  &9.'°  ..u 
increments  of  about  2.3°.  Indices  between  those  listed  can  be 
used  to  interpolate  a  better  approximation  of  a  desired 
elevation. 

The  open  circuit  impedance  matrices  are  calculated  to 
determine  the  current  distribution  among  the  dipole  elements. 
The  matrix  of  mutual  impedances  between  the  dipole  elements  is 
represented  as 
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Each  main  diagonal  term  (Zltl)  is  the  self  impedance  of  the  i*+i 
element.  The  off-diagonal  terms  (Zik)  are  the  mutual 
impedances  between  the  im+l  and  k“+l  elements.  Subscripts  i 
and  k  are  matrix  indices  (0, 1, 2 . . . .NR+ND) . 

The  matrix  of  open  circuit  mutual  impedances  between  the 
actual  elements  and  their  images  is  represented  as 


ZIo.  0 

ZIo,  i 

ZIo,  2  *  •  • 

•  ZI 0i  iN.x) 

ZIx.o 

zi,.. 

ZI  1,2  •  * 

•  zi,t  (N.X) 

[ZI]  = 

ZI2l  o 

ZI2.1 

zi2  2  ... 

•  zi2llf,1) 

•  > 

ZI(N-l)  ,0 

ZI{N- 1)  ,1 

ZI (N- 1)  .2  ’  •  • 

•  • 

■  ZIfff.  ,i,(N-l) 

The  main  diagonal  terms  (ZI^)  are  the  mutual  impedances 
between  the  i^+l  element  and  its  image,  and  the  off-diagonal 
terms  (ZIU)  are  the  mutual  impedances  between  the  i^+l  element 
and  the  image  of  the  k^+l  element. 

The  main  diagonal  self -impedance  terms  of  equation  17.1 
are  found  from  [Ref  6:  pp.  205-206] 

[Q(U0)  -2Q{UX))  (17.3) 


z  60 

1,1  1-cos  (2p/ii) 


+  ej2*h‘  [£>(  VJ,)  -2Q{VX)  J+2[(?(  £7^)  -Q{UX)  -Q{Vl)  ]+2!?(  O^q)  [l+cos  (2Phi)  ]  ] 

X  X 

where  Q(x)  =Ci  ( x )  -j-Si  ( x )  =  J  cos-^)~  dy-j-f  s-n  dy  (17.4) 

m  y  0  y 

£70  =  P  [v^2zad|+4hi-2hi]  U'0=y/2$radi  V0=P  [v/2radI2  +  4hI2+2hi] 


Ux- P  yTzadf+hj-hi 


V^-p  [^2radf +h2+hi] 
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The  orf -diagonal  mutual  impedance  terms  of  equation  17.1 
are  given  by  CRef  11] 


2  =2  .  = _ — — _ 

i.k  k.i  Cos  [p  (h2-hk)  ]  -cos  [p  (h1+hk)  ] 


(17.5) 


-Qitfj  -o(u2)]^ej^,hrhk)[0(  /0)  -Q(Vl)  -Q{V2)) 
+  e-jp'Wjc,(C7/o_£,(yi)  -cXU2)]+e^'A*''h*'[C)(y/0)  ~Q(U2)  -Q(V2)) 
+20(  PDiJc)  (cos  [p  (hrhk)  ]  +cos  [  P  (h^hk)  ])] 


X  X 

where  Q(x)  =Ci  (x)  -j-Si  ( x )  =  f  — s  ■  dy  -  j-f  dy 

t  y  {  y 


The  Dik  terms  are  the  distances  between  the  i^+l  and  k*+i 
dipole  elements  and  are  given  by 

Jc-l 

Dl.)CDk.l=Y,dn  (17.6) 

The  mutual  impedances  between  the  actual  elements  and 
their  images  must  also  be  computed  to  determine  the  elemental 
current  distribution.  The  mutual  impedance  terms  of  equation 
17.2  (ZIiik)  are  given  by  equation  17.5,  except  the  Di-k  terms  in 
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the  Q(x)  arguments  are  replaced  by  dij>k  terms  which  represent 
the  distance  from  the  i^+l  element  to  the  image  of  the  k^+l 
element  and  are  given  by  [Ref  6:  pp.  205-208] 

diUk=sjH^Dlk  (17.7) 


The  function  Q(x)  given  by  equations  17.4  is  defined  in 
terms  of  sine  [Si (x) ]  and  cosine  [Ci(x)]  integrals.  Mathcad 
is  incapable  of  evaluating  Si(x)  and  Ci(x)  directly,  so  a 
series  expansion  and  polynomial  approximation  are  written  into 
the  Mathcad  code  for  both  the  sine  and  cosine  integrals  to 
evaluate  the  functions  Q(x).  The  arguments  for  Q(x)  are 
always  real,  so  it  is  unnecessary  to  find  Si  (x)  and  Ci(x) 
expressions  valid  for  complex  numbers. 

For  arguments  less  than  one,  Si(x)  and  Ci(x)  are 
evaluated  by  series  expansions  given  by  [Ref  9:  p.  232] 


Si  (x)  =£ 

73=0 


(-1)”  x2"*1 
(2/i+D  (2/3+1)  ! 


Ci  (x)  =v+ln(x)  +£ 

n-1 


(~l)n  x2n 
(2/3)  (2/3)  .' 


where  y=0 . 5772156649 


(17.8) 

(17.9) 


For  arguments  greater  than  one,  Si(x)  and  Ci  (x)  are 
expressed  in  terms  of  auxiliary  functions,  f  (x)  and  g(x), 
given  by  Ref  9:  p.232] 

Si  (x)  =-j-f  (x)  cos  (x)  -g(x)  sin(x)  (17 . 10) 

Ci  (x)  =f (x) sin(x) -g(x) cos (x)  (17 . 11) 
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The  auxiliary  functions  f  (x)  and  g(x)  are  evaluated  by  the 
polynomial  approximations  [Ref  9:  p.  233] 


f(x)=-± 


x8+a1x6+a2x4  +  a3x2  +  a4) 


x[  x6+b1x6+b2x4+b2x2+bi  J 
where  |e(x)|  <  5-10"7 


+e  fx) 


a1=38 . 027264 


£^=40.021433 


a2  =  265 . 187033 


£>2  =  322. 624911 


a3=335. 677320 


jb3  =  570. 236280 


a4=38. 102495 

/  -,8 


b4=157 .105423 


g{x) 


x°+a1xD+a2x‘+a3x^+a4 


y  xB+b1x6+b2x4+b2x2+bi) 


+e(x) 


where 


a,  =42. 242855 


|e(x)  |  <  3*1  O'7 
^=48.196927 


a2  =  302. 757865 
a3  =  352 . 018498 
a4=21. 821899 


2?2=482. 485984 


jb3  =  1114 .978885 


(17 .12) 


(17 .13) 


jb4  =449 .690326 

When  the  mutual  impedances  are  found,  the  array's  base 
current  distribution  is  computed  from  [Ref  10  p.  258-259] 

[IS]  ={[Z]  +rh'0[ZI]yl[VREF]  (17.14) 

where 


^h,0~^h  I  0-0 = 


1-n 


1+n 


and  n  is  the  complex  index  of  refraction.  The  entries  of  the 
[IB]  matrix  are  the  base  currents  of  the  dipole  elements,  and 
[VREF]  is  the  voltage  matrix  to  which  the  [IB]  distribution  is 
referenced.  Since  the  Yagi-Uda  has  only  one  driven  element 
(corresponding  to  the  NR*  index) ,  and  since  a  sinusoidal 
voltage  response  with  a  maximum  of  unity  is  assumed  across  the 


input  terminals,  [VREF]  is  a  1  X  NR+ND  matrix  with  a  1  as  the 
NR*  entry  and  NR+ND- 1  zeroes  in  the  remaining  positions.  The 
electric  field  equations  are  referenced  to  the  current  maxima 
values  for  the  dipole  elements,  [I] ,  which  are  determined  frcm 
the  [IB]  matrix  using  [Ref  6  p.  209] 


IB, 

sin(p  hj) 


(17.15) 


The  mutual  impedance  calculations  are  not  valid  when  any 
element  length  :  s  an  exact  integer  multiple  of  the  wavelength. 
When  this  occurs,  there  is  a  singularity  error  in  the  mutual 
impedance  calculations,  and  it  is  necessary  to  vary  the 
frequency  such  that  no  element  is  exactly  an  integer  multiple 
of  the  wavelength.  The  change  required  is  only  a  percent  or 
two,  and  the  predicted  radiation  parameters  are  still  a  good 
estimate  of  those  for  the  original  frequency. 

The  electric  field  for  the  horizontal  Yagi-Uda  array  is 
obtained  for  each  element  in  a  manner  analogous  to  that  used 
for  the  vertical  dipole  in  Chapter  4.  The  equations  for  the 
individual  elements  are  simplified  into  two  expressions  for 
the  array's  electric  field  components  by  the  array  factor. 
The  Yagi-Uda' s  array  factor  is  given  by  [Ref  6:  p.  204] 


s  =  E  i’i[cos[Phi(sin(0)  sin(<}>)  )]-cos  (P^)  ]  e■^1'isin<0,sin^♦,  (17.16: 
1*0 

where  Yi  is  the  y  coordinate  for  the  ieh+ 1  element  given  by 

i-x 

^i=E  dn  (17.17) 

73-0 
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The  array  factor  is  multiplied  by  the  element  factors  to 
find  the  equations  for  the  total  radiated  electric  field 
distribution  of  the  horizontal  Yagi-Uda  array, 

o  — - _ e  ''pWy-'os " 81  ( S)  •  (17.18) 

r  P  i-sir.M0)  cos*  (4>) 

'coS(0)(l-r.,e--J^coi  0>(l-rjFee-;‘^cos,9)- 


v  ~  ~  s (sin1  (0)  -  <ni-8i??±*Lcos  (0) 
n  *  n  * 


sin  < «t>) 


c  -  _ "¥• _  g 

P  i-sin2  (6)  cos*  (0) 


,’PHncos  (8) 


(5)-  (17.19) 


The  first  two  terms  inside  the  brackets  of  each  equation  are 
the  space  wave,  the  third  terms  are  the  surface  wave,  and  I, 
is  the  current  term  for  each  element  based  on  a  sinusoidal 
voltage  across  the  input  terminals  with  a  maximum  of  unity. 

The  requested  inputs  are  used  to  calculate  the  following 
variables  using  a  constant  <t>  of  0  and  <p-ir/2  for  312  discrete 
values  of  6  which  are  equally  incremented  from  - ir/2  to  tt/2  : 


h, . half-length  of  the  i*+l  element  (meters) 

Y, . y  coordinate  of  the  iUl+l  element 

I, . current  for  the  i'h+l  element 

X, . wavelength  of  the  operational  frequency  (meters) 

& . free  space  wavenumber  for  operational  frequency 

n . index  of  refraction  \ 

Pc . complex  numerical  distance  (Vertical  polarization) 

Pm  . complex  numerical  distance  (horizontal  polarization) 

r, . vertical  reflection  coefficient 

I\ . horizontal  reflection  coef f icjient 

F, . vertical  surface  wave  attenuation  factor 

F„  . horizontal  surface  wave  attenuation  factor 
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The  calculated  variables  are  used  to  evaluate  the  far- 
field  space  wave  and  surface  wave  for  the  discrete  values  of 
6  with  0=0  and  0=7r/2.  The  space  wave  and  surface  wave  results 
are  combined  for  corresponding  values  of  6  to  obtain  the  total 
radiated  electric  field  distribution  for  the  0=0  and  0=7r/2 
vertical  planes.  The  space  wave,  surface  wave,  and  total 
electric  field  results  are  then  normalized  with  respect  to  the 
maximum  field  intensity  of  each,  and  the  normalized  magnitudes 
are  plotted  for  each  discrete  6  to  depict  the  radiation 
patterns.  The  Mathcad  horizontal  Yagi-Uda  array  application 
computes  the  space  wave,  surface  wave,  and  total  electric 
field  radiation  patterns  and  radiation  parameters  in  the  0=0 
and  0=7r/2  vertical  planes. 

The  variables  corresponding  to  the  selected  elevation 
angle  index  (w)  are  used  to  evaluate  the  radiated  electric 
field  components  for  the  space  wave  and  surface  wave  at  the 
fixed  elevation  angle  ($v)  as  0  varies  from  0  to  2t r  in  312 
equal  increments.  The  horizontal  radiation  patterns  are  then 
plotted  for  the  space  wave,  surface  wave,  and  total  radiated 
electric  field  just  as  those  for  the  vertical  planes. 

Equations  3.8  and  3.9  are  used  to  integrate  equations 
17.18  and  17.19  over  the  hemispherical  Gaussian  surface  above 
the  ground  plane  at  a  fixed  radius  (R)  from  the  array  to  find 
total  average  radiated  power  (P^)  .  Since  the  Yagi-Uda  current 
distribution  is  referenced  to  a  sinusoidal  voltage  input  with 
a  maximum  of  unity,  the  radiation  resistance  can  not  be 
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calculated  from  equation  3.10  as  in  previous  applications. 
Instead,  radiation  resistance  must  be  found  with  respect  to 
the  unity  voltage  input  using 


rad 


=_M 


2  R 


or 


rad 


rad 


,JZl? 


2  P 


(17.20) 


rad 


With  the  discrete  values  of  the  electric  field  and  total 
average  radiated  power,  the  Mathcad  application  predicts  the 
following  radiation  characteristics  from  Chapter  3: 


Rrad  . radiation  resistance  (Ohms) 

D0 . directivity 

EIRP . effective  radiated  isotropic  power  (Watts) 

A^  . maximum  theoretical  effective  area  (square  meters) 

l,^  . maximum  theoretical  effective  length  (meters) 

P . numerical  distance  (vertical  polarization,  0=90°) 

Pm  . numerical  distance  (horizontal  polarization,  0=90°) 


Angle^  .elevation  angle  of  maximum  directive  gain  (degrees) 

The  directivity  (D0)  ,  effective  isotropic  radiated  power 
(EIRP) ,  maximum  theoretical  effective  area  (A^J  ,  maximum 
theoretical  effective  length  (1^)  ,  and  elevation  angle  of  the 
directivity  (Angle,^)  are  determined  for  both  vertical  planes. 

As  an  example,  the  Mathcad  horizontal  Yagi-Uda  array 
application  was  executed  with  the  following  inputs: 


number  of  reflectors 
number  of  directors 
element  lengths 
element  separations 
element  radii  0.001, 

height  of  antenna 
frequency 

distance  from  the  antenna 
relative  dielectric  constant 
ground  conductivity 
elevation  angle  index 


1 

1 

15.6,  15,  and  14 
7.5  and  6.5 
0.001,  and  0.001 
8.0 
10 -106 
3000 
72 
4 

285 


meters 

meters 

meters 

meters 

Hertz 

meters 


(-17°) 
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Yagi-Uda  space  wave  radiation 
pattern  for  horizontal  plane. 


radiation  pattern  for 
horizontal  plane. 


The  predicted  radiation  parameters  based  on  a  one  volt 
response  across  the  input  terminals  are: 


Total  power  radiated  (Watts) 
Radiation  resistance  (Ohms) 
Numerical  distance  (vertical) 
Numerical  distance  (horizontal) 


0.00482 

103.760 

0.04363 

2.26-106 


<6*  0 

Directivity 

1.679 

14.631 

EIRP  (Watts) 

0.0081 

0.0705 

Max  eff  area  (sq  meters) 

120.26 

1047.9 

Max  eff  length  (meters) 

11.506 

33.966 

Angles  (degrees) 

89.7 

41.714 

These  results  are  consistent  with  expectations  for  this 


particular  configuration.  Appendix  F  contains  computer 
hardcopies  of  additional  example  calculations  for  the 


horizontal  Yagi-Uda  array  and  compares  predicted  radiation 
parameters  to  those  expected  based  on  previous  calculations. 


XVIII.  REMARKS  AND  CONCLUSION 


This  thesis  is  the  culmination  of  nine  months  of  research 
and  computer  programming .  No  new  electric  field  equations 
were  derived  exclusively  for  this  project.  It  was  directed  by 
NAVMARINTCEN  prior  to  commencing  that  existing  electric  field 
equations  would  provide  the  basis  for  the  radiation  parameter 
predictions.  The  expressions  used  in  this  thesis  are  all 
previously  derived  by  pioneers  in  antenna  radiation  theory 
such  as  Sommerfeld,  Norton,  Cox,  King,  Ma,  and  Walters.  The 
equations  are  compiled  here  solely  for  the  purpose  of 
crediting  the  source  references  and  to  describe  how  the 
Mathcad  code  was  assembled. 

The  accuracy  of  the  predicted  radiation  characteristics 
is  totally  dependent  upon  the  extent  to  which  the  equations 
used  realistically  model  the  actual  radiated  electric  fields 
of  the  antennas.  It  is  near  impossible  to  obtain  analytic 
results  which  accurately  model  an  antenna's  radiation 
characteristics  in  all  cases.  There  are  too  many  operational 
and  environmental  variables  to  obtain  one  general,  all-purpose 
expression.  Even  if  all  possible  variables  could  be  accounted 
for  in  a  single  expression,  computer  processing  times  would  be 
unacceptably  long.  The  equations  used  herein  are  simple 
enough  to  be  evaluated  by  Mathcad  in  a  reasonable  length  of 
time,  but  the  radiation  parameter  predictions  are  accurate 
enough  to  be  useful  for  antenna  analysis. 
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As  addressed  in  appendices  A-F,  the  equations  presented 
in  previous  chapters  and  evaluated  by  the  Mathcad  applications 
provide  radiation  parameters  consistent  with  expectations  for 
the  inputs  which  have  been  executed  to  date.  All  antenna 
configurations  which  have  been  computed  exhibit  radiation 
characteristics  consistent  with  other  computational  programs 
and  empirically  obtained  patterns  and  parameters.  There  is  no 
way  to  test  the  accuracy  of  the  Mathcad  predictions  other  than 
empirical  measurements  for  each  antenna  configuration,  clearly 
a  task  beyond  the  scope  and  purpose  of  this  report.  However, 
adequate  analysis  of  the  Mathcad  results  has  been  provided  to 
demonstrate  that  the  applications  provide  very  good  estimates 
of  antenna  radiation  parameters  as  a  base  level  analysis  tool. 

Most  of  the  applications  have  computations  times  under 
ten  minutes.  However,  even  with  the  simplified  equations  used 
by  Mathcad,  some  of  the  applications  can  take  upwards  of  two 
hours  to  compute  on  a  33  MHz  80386  PC.  The  applications  for 
the  rhombic,  double  rhomboid,  and  array  (with  more  than  four 
to  five  elements)  antennas  can  take  an  exorbitantly  long  time 
for  power  calculations.  For  this  reason,  the  applications 
should  be  run  on  the  fastest  PC  available;  a  50  MHz  80486  is 
preferable.  Computation  times  can  also  be  reduced  by 
computing  only  the  radiation  patterns  until  it  is  determined 
which  frequencies  are  of  greatest  interest.  Then  the  average 
power  calculations  for  those  frequencies  can  be  executed. 
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The  Mathcad  applications  were  written  for  NAVMARINTCEN  to 
provide  antenna  analysts  with  the  capability  to  predict 
radiation  parameters  based  solely  on  antenna  physical 
dimensions  and  ground  properties.  The  Mathcad  results 
computed  thus  far  are  consistent  with  expectations  and  are 
most  likely  providing  accurate  predictions  of  actual  antenna 
radiation  characteristics.  Only  extended  use  of  the 
applications  and  empirical  confirmation  of  results  will  prove 
out  the  accuracy  of  the  Mathcad  code.  Results  obtained  thus 
far  certainly  justify  continued  use  of  the  Mathcad  code  as  an 
antenna  analysis  tool. 


APPENDIX  A; 

VERTICAL  DIPOLE  ARRAY  COMPUTER  OUTPUT 

This  appendix  contains  computer  hardcopies  from  the  Mathcad 
vertical  dipole  application  which  show  the  input  values  and 
predicted  radiation  characteristics  for  two  sample  calculations. 
The  configuration  is  a  half-wave  dipole  at  one  quarter  wavelength 
above  soil  (er=10  and  cr=10'2)  for  the  first  example  and  the  same 
configuration  above  seawater  (er=72  and  <r=4)  for  the  second. 
Reference  6  [p.  91]  provides  the  radiation  patterns  and  gain 
predictions  from  several  sources  for  the  configuration  in  the 
first  example. 

The  elevation  of  maximum  directive  gain  is  slightly  higher 
for  the  Mathcad  output  than  for  those  given  in  reference  6,  but 
the  overall  radiation  patterns  are  very  similar.  The  4.17  (6.2 
dB)  value  of  directivity  for  Mathcad  is  quite  a  bit  higher  than 
the  maximum  gain  values  of  about  zero  dB  in  reference  6 . 
However,  since  a  half-wa\ »  dipole  has  a  free-space  directivity  of 
1.64  (2.1  dB)  ,  one  would  expect  the  actual  directivity  to  be 
closer  to  the  Mathcad  prediction  because  of  the  effect  of  the 
reflected  wave  (constructive  and  destructive  interference)  on  the 
space  wave,  and  the  lower  total  average  radiated  power  resulting 
from  ground  plane  losses.  The  seawater  example  yields  results 
consistent  with  expectations.  With  respect  to  the  soil  example, 
the  seawater  example's  directivity  is  slightly  higher  due  to  a 
stronger  reflected  wave,  and  the  surface  wave  is  stronger  at 
grazing  angles  (0-90°)  due  to  higher  conductivity. 


VERTICAL  DIPOLE 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated 
with  vertical  thin-wire  dipole  antennas  (diameter  «  wavelength).  The  antenna  is 
mounted  vertically  along  the  z-axis  in  a  rectangular  coordinate  system  with  the  feed 
at  the  ante,  na's  center  at  a  set  height  above  the  surface.  Required  inputs  are  the 
antenna  length,  feed  height  above  the  surface,  transmitted  frequency,  distance  from 
the  antenna,  the  conductivity  and  dielectric  constant  of  the  surface  below  the 
antenna.  The  planar  earth  model  is  assumed  in  predicting  radiation  patterns. 
Predicted  operating  frequencies  assume  that  the  antenna  is  a  quarter*  wavelength, 
half-wavelength,  three-quarter-wavelength,  or  full-wavelength  dipole.  A  sinusoidal 
current  input  with  a  maximum  of  unity  is  assumed.  All  radiation  patterns  are 
normalized  with  respect  to  the  maximum  electric  field  intensity  transmitted  by  the 
antenna.  Plotted  radiation  patterns  are  valid  for  any  vertical  plane  passing  through 
the  antenna  axis.  The  electric  field  magnitudes  to  which  the  Radiation  Patterns  are 
normalized  are  displayed  below  their  respective  plots.  Polarization  is  vertical  for  all 
vertical  dipoles. 


Input  the  Dipole 
length  in  meters 


Wavelengths  and  Frequencies 


Possible  Operating  Frequencies  (Hertz) 


^  Dipole: 

4 

ft  =  1.5*  107 

X  , 

-  Dipole: 

h  =  3- 107 

Dipole: 

4 

f3  =  4.5- 107 

X  Dipole: 

r4  =  6- io7 

Input  the  operating 
Frequency  (Hertz) 

f5  30-  i06 

Input  the  Distance 
from  the  Antenna 
(meters) 

R  :=  3000 

Input  the  ground 

cr  =  10 

Input  the  ground 

o  ■«  10' 2 

Dielectric  Constant 

Conductivity 

Input  the  Height  of 
the  Antenna  Feed 


(meters) 
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2.5 


Radiation  Patterns  valid  for  any  Vertical  Plane  passing  through  the  Antenna 
Space  Wave  Radiation  Pattern 


Maximum  Space  Wave  Electric  Field 

Intensity  (Volts  per  meter)  max(MagEl)  =  0.01733 


153 


Numerical  Distance  for 
Vertical  Polarization 


PeQ|  =  74.95924 


Elevation  Angle  of  Maximum 
Power  Gain  (Degrees) 


Angle  3  = 


19.39778 

19.39778 


VERTICAL  DIPOLE 


Tliis  application  calculates  far  field  radiation  patterns  and  parameters  associated 
with  vertical  thin-wire  dipole  antennas  (diameter  «  wavelength).  The  antenna  is 
mounted  vertically  along  the  z-axis  in  a  rectangular  coordinate  system  with  the  feed 
at  the  antenna's  center  at  a  set  height  above  the  surface.  Required  inputs  are  the 
antenna  length,  feed  height  above  the  surface,  transmitted  frequency,  distance  from 
the  antenna,  the  conductivity  and  dielectric  constant  of  the  surface  below  the 
antenna.  The  planar  earth  model  is  assumed  in  predicting  radiation  patterns. 
Predicted  operating  frequencies  assume  that  the  antenna  is  a  quarter-  wavelength, 
half-wavelength,  three-quarter-wavelength,  or  full-wavelength  dipole.  A  sinusoidal 
current  input  with  a  maximum  of  unity  is  assumed.  All  radiation  patterns  are 
normalized  with  respect  to  the  maximum  electric  field  intensity  transmitted  by  the 
antenna.  Plotted  radiation  patterns  are  valid  for  any  vertical  plane  passing  through 
the  antenna  axis.  The  electric  field  magnitudes  to  which  the  Radiation  Patterns  are 
normalized  are  displayed  below  their  respective  plots.  Polarization  is  vertical  for  all 
vertical  dipoles. 

Input  the  Dipole  1  =  5  h  -  - 

length  in  meters  2 


Wavelengths  and  Frequencies 


A,  =  41 

X2  21 

X  -  41 
*3  '  t 

X4  1 

-n 

ii 

_ri« 

f  *  ± 

2  *2 

'|-f 

A3 

a4 

Possible  Operating  Frequencies  (Hertz) 

~  Dipole: 

4 

ft  =  1.5- 107 

^  Dipole: 

f2  =  3- 10 

Dipole: 

4 

f3  -  4.5- 107 

X  Dipole: 

f4  =  6- 10 

Input  the  operating 
Frequency  (Hertz) 


f5  :=  30- 10$ 


Input  the  Distance 

from  the  Antenna  r  =  3000 

(meters) 


Input  the  ground 
Dielectric  Constant 


Input  the  ground 
Conductivity 


Input  the  Height  of 
the  Antenna  Feed 
(meters) 
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Radiation  Patterns  valid  for  any  Vertical  Plane  passing  through  the  Antenna 
Space  Wave  Radiation  Pattern 


Maximum  Space  Wave  Electric  Field 

Intensity  (Volts  per  meter)  max(MagEl)  =  0.03394 


161 


i 


power  ••=  30- 


* 

2 


.'0 


cos(/9- h  cos(  £) )  _  cos(/3  h) 
sin(  ?) 


'  HQ-coi(a)  -|-p-(*  * 


R  -  H0  cos(O 


(|E3,I)2 

SqE3|  :=  — - — — 

M  1  2(120x) 


_.  ..  4-x-R z  max(sqE3) 

Directivity  =  - — 3 — - 

power 


Radres  =  2 ‘power 


Total  Power  Radiated  (Watts) 


Radiation  Resistance  (Ohms) 


Directivity  (maximum  Power  Gain 
assuming  100%  Antenna  efficiency) 


Effective  Isotropic  Radiated 
Power  (EIRP)  (Watts) 


Maximum  Effective  Area 
(Along  Redial  of  Directivity) 
(square  meters) 


Maximum  Effective  Length 
(Along  Radial  of  Directivity) 
(meters) 


power  =  37.64371 

Radres  =  75.28741 

Directivity  =  4.59967 

Directivity  power  =  173.14876 


/A«')  •  Directivity 

-  =  36.60304 

4‘X 


Radres-  (Ag)  •  Directivity 
-  =  5.40735 


480-x‘ 


R  *  H0 
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Numerical  Distance  for 
Vertical  Polarization 


=  0.39252 


Elevation  Angle  of  Maximum 
Power  Gain  (Degrees) 


Ang!e3 


/  8.55784  V 
\  8.55784/ 


APPENDIX  B: 

VERTICAL  MONOPOLE  ARRAY  COMPUTER  OUTPUT 


This  appendix  contains  computer  hardcopies  from  the  Mathcad 
vertical  monopole  application  which'  show  the  input  values  and 
predicted  radiation  characteristics  for  two  sample  calculations. 
The  configuration  is  a  quarter- wave  monopole  above  soil  (er=lO  and 
a=10‘2)  for  the  first  example  and  the  same  configuration  above 
seawater  (er=72  and  <r= 4)  for  the  second.  Reference  6  [p.  89] 
provides  the  radiation  patterns  and  gain  predictions  from  several 
sources  for  the  configuration  in  the  first  example. 

The  elevation  of  maximum  directive  gain  is  slightly  higher 
for  the  Mathcad  output  than  for  those  given  in  reference  6,  but 
the  overall  radiation  patterns  are  very  similar.  The  3.24  (5.1 
dB)  value  of  directivity  for  Mathcad  is  quite  a  bit  higher  than 
the  maximum  gain  values  of  about  zero  d3  in  reference  6 . 
However,  since  a  half-wave  dipole  has  a  free-space  directivity  of 
1.64  (2.1  dB) ,  the  actual  directivity  of  a  quarter-wave  monopole 
above  a  ground  plane  should  be  closer  to  the  Mathcad  prediction 
because  of  the  effect  of  the  reflected  wave  (constructive  and 
destructive  interference)  on  the  space  wave  and  the  lowdr  total 
average  radiated  power  resulting  from  ground  plane  losses.  The 
seawater  example  yields  results  consistent  with  expectations. 
With  respect  to  the  soil  example,  the  seawater  example's 
directivity  is  slightly  higher  due  to  a  stronger  reflected!  wave, 
and  the  surface  wave  is  stronger  at  grazing  angles  (0*» 90°)  due  to 
higher  conductivity. 


\ 
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VERTICAL  MONOPOLE 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated 
with  vertical  thin-wire  monopole  antennas  (diameter  «  wavelength).  The  antenna 
is  mounted  vertically  along  the  z-axis  in  a  rectangular  coordinate  system  with  the 
feed  at  the  origin.  Required  inputs  are  the  antenna  length,  transmitted  frequency, 
distance  from  the  antenna,  the  conductivity  and  dielectric  constant  of  the  surface 
below  the  antenna.  The  planar  earth  model  is  assumed  in  predicting  radiation 
patterns.  Predicted  operating  frequencies  assume  that  the  antenna  is  an  eighth- 
wavelength,  quarter-wavelength,  three-eighths-wavelength,  or  haif-wavetength 
monopole.  A  sinusoidal  current  input  with  a  maximum  of  unity  is  assumed.  All 
radiation  patterns  are  normalized  with  respect  to  the  maximum  eiectic  field  intensity 
transmitted  by  the  antenna.  The  electric  field  magnitudes  to  which  the  patterns  are 
normalized  are  displayed  below  (heir  respective  plots.  Radiation  patterns  are  valid 
fo.  any  vertical  plane  passing  through  the  z-axis,  because  the  radiation  pattern  is 
symmetrical  with  respect  to  phi.  Polarization  is  vertical  for  all  vertical  monopoles. 

I 


Input  the  Monopole 
length  in  meters 


h  «  7.5 


1  •  2-h 


Wavelengths  and  Frequencies 


*1 

■«  41 

X2 

=  21 

X,  •  i! 

x4 

II 

f1 

it 

f2 

tl 

'3  ■  f 

A3 

f4 

li 

Possible  Operating  Frequencies  (Hertz) 


g  Monopole:  f|  =  5- 10^  ^  Monopole:  f<>  =  1*10* 


Monopole:  fg  =  1.5*10* 


^  Monopole:  f^  =  2*10* 


Input  the  operating  jq*  10^  I°Put  the  Distance  from  ^  ^goo 


Frequency  (Hertz) 


the  Antenna  (meters) 


Input  the  ground  a  ..  jq  Input  the  ground  g  _  jq- 


Dieiectric  Constant 


Conductivity: 


/ 


power 


T5 

R2 


-'0 


cos(0-hcos(f))  -  cos(jS-h)  ... 


+  J 


sin(/3  h  cos(n)  - 

+  (-(cos(f)sin(j8-h))) 


sb(f) 


[  cos(/J-h 

•cos(f))  -  cos(/3-h)  ... 

9  | 

+ 

- 

i* 

sin(/3-h-cos(C))  |  1 

n-cos(f)  -  J 

nz  -  sin(0 

1  +  (-(cos(O-sin(0-h)) 

n*-cos(  O 


4-  Ja 


sin(C)1 


sin(f) 


sqE3,  = 


(lE3i|)2 

2(120*) 


Directivity  = 


4  *  R  -max(sqE3) 


power 


Radres  -  2*  power 


Total  Power  Radiated  (Watts)  power  =  5.74454 

Radiation  Resistance  (Ohms)  Radres  =  1 1.48908 


Directivity  (maximum  Power  Gain 
assuming  100%  Antenna  Efficiency) 


Effective  Isotropic  Radiated 
Power  (EIRP)  (Watts) 


Maximum  Effective  Area 
(Along  Radial  of  Directivity) 
(square  meters) 


Directivity  =  3.23898 


Directivity- power  =  18.60648 


(Aj)  2  Directivity 

4* 


=  231.97519 
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/ 

/ 


Maximum  EITeclive  Length 
(Along  Radial  of  Directivity) 
(meters) 


A 

Radres-(Atj) 4-  Directivity 


480-x2 


=  5.31775 


Numerical  Distance  for  ,  , 

Vertical  Polarization  |*eo|  ~  14.91116 


Elevation  Angle  of  Maximum 
Power  Gain  (Degrees) 


Angle3  = 


25.67353  \ 
25.67353/ 


172 


VERTICAL  MONOPOLE 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated 
with  vertical  thin-wire  monopole  antennas  (diameter  «  wavelength).  The  antenna 
is  mounted  vertically  along  the  z-axis  in  a  rectangular  coordinate  system  with  the 
feed  at  the  origin.  Required  inputs  are  the  antenna  length,  transmitted  frequency, 
distance  from  the  antenna,  the  conductivity  and  dielectric  constant  of  the  surface 
below  the  antenna.  The  planar  earth  model  is  assumed  in  predicting  radiation 
patterns.  Predicted  operating  frequencies  assume  that  the  antenna  is  an  eighth- 
wavelength.  quarter-wavelength,  three-eighths-wavelength.  or  half- wavelength 
monopole.  A  sinusoidal  current  input  with  a  maximum  of  unity  is  assumed.  All 
radiation  patterns  are  normalized  whh  respect  to  the  maximum  eleclic  field  intensity 
transmitted  by  the  antenna.  The  electric  field  magnitudes  to  which  the  patterns  are 
normalized  are  displayed  below  their  respective  plots.  Radiation  patterns  are  valid 
for  any  vertical  plane  passing  through  the  z-axis,  because  the  radiation  pattern  is 
symmetrical  whL  respect  to  phi.  Polarization  is  vertical  for  all  vertical  monopoles. 


Input  the  Monopole 
length  in  meters 


h  -  7.5 


I  •  2  h 


Wavelengths  and  Frequencies 


X,  =  41 


f<  ■  r, 


X2  =  21 


f  :  JL 

X2 


r 


4-1 

3  *  3 

c 


Possible  Operating  Frequencies  (Hertz) 


x4  --  1 


^  Monopole:  fj  =  5*l0* 
Monopole:  f j  =  i .5- 10^ 


-  Monopole:  f2  =  1  •  107 

4  ■ 


-  Monopole:  f|  =  2- 10* 


Input  the  operating  r  ,  10- 10^  Input  the  Distance  from  « 
Frequency  (Hertz)  5  the  Antenna  (meters) 


3000 


Input  the  ground  ff  ^  Input  the  ground 
Dielectric  Constant  Conductivity: 
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/ 

/ 


Combined  Space  and  Surface  Wave  Pattern 


-1  -0.3  0  0.5  1 


.  h3j 

Maximum  Radiated  Electric  Field 

Intensity  (Volis  per  meter)  max(MagE3)  =0.01976 


176 


177 


/ 


Maximum  Effective  Length 
(Along  Radial  of  Directivity) 
(meters) 


Numerical  Distance  for 
Vertical  Polarization 


i - ; - 

Radres-^X^  4  Directivity 
480  *2 


PeQ|  =  0.04363 


Elevation  Angle  of  Maximum 
Power  Gain  (Degrees) 


9.4327 


178 


APPENDIX  C: 

HORIZONTAL  DIPOLE  ARRAY  COMPUTER  OUTPUT 


This  appendix  contains  computer  hardcopies  from  the  Mathcad 
horizontal  dipole  application  which  show  the  input  values  and 
predicted  radiation  characteristics  for  four  sample  calculations. 
The  configuration  is  a  half-wave  dipole  at  one-  quarter,  one-half, 
and  three-quarters  wavelengths  above  soil  {er=lO  and  cr=10'3)  for  the 
first  three  examples,  respectively,  and  a  half-wave  dipole  at  one- 
quarter  wavelength  above  seawater  {er=72  and  a- 4)  for  the  fourth. 
Reference  6  [pp.  86-88]  provides  the  radiation  patterns  and  gain 
predictions  from  several  sources  for  the  configurations  in  the 
first  three  examples. 

The  radiation  patterns  and  maximum  directive  gain  computed  by 
the  first  three  Mathcad  examples  are  almost  identical  to  the  those 
given  in  reference  6.  The  directivity  predictions  are  slightly 
higher  for  Mathcad,  but  the  overall  similarity  between  the 
predicted  radiation  characteristics  is  noteworthy.  The  seawater 
example  also  yields  results  consistent  with  expectations.  With 
respect  to  the  soil  example,  the  seawater  example's  directivity  is 
slightly  higher  due  to  a  stronger  reflected  wave,  and  the  surface 
wave  is  slightly  stronger  at  grazing  angles  (0-90°)  due  to  higher 
conductivity.  For  horizontal  polarization,  the  higher  conductivity 
surface  below  the  antenna  does  not  result  in  a  greatly  enhanced 
surface  wave  and  increased  directivity  as  it  does  for  vertical 
polarization. 
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HORIZONTAL  DIPOLE 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated  with 
horizontal  thin-wire  dipole  antennas  (diameter  «  wavelength).  The  antenna  is 
mounted  above  and  parallel  to  the  x-axis  in  a  rectangular  coordinate  system.  The 
feed  is  at  the  center  of  the  antenna  at  a  set  height  directly  above  the  origin. 

Required  inputs  are  (he  antenna  length,  feed  height  above  the  surface,  transmitted 
frequency,  distance  from  the  antenna,  the  conductivity  and  dielectric  constant  of  the 
surface  below  the  antenna.  The  planar  earth  model  is  assumed  in  predicting 
radiation  patterns.  Predicted  operating  frequencies  assume  that  the  antenna  is  a 
quarter- wavelength,  half-wavelength,  three-quarter -wavelength,  or  full  wave¬ 
length  dipole.  A  sinusoidal  current  input  with  a  maximum  of  unity  is  assumed.  All 
radiation  patterns  are  normalized  with  respect  to  the  maximum  electric  field 
intensity  transmitted  by  the  antenna  in  the  plane  of  interest.  The  electric  field 
magnitudes  to  which  the  patterns  are  normalized  are  displayed  below  their 
respective  plots.  The  radiation  patterns  are  plotted  for  the  phi=0  and  phi=pi/2 
vertical  planes  parallel  and  perpendicular  to  the  x-axis  respectively.  A  horizontal 
radiation  pattern  is  plotted  at  an  elevation  selected  by  the  index  from  the  elevation 
angle  index  table.  Polarization  is  a  combination  of  vertical  and  horizontal 
depending  upon  spatial  orientation  relative  to  the  dipole. 

Input  the  dipole  l  =  15  h  =  - 

length  in  meters  2 


Input  the  index  of  the  elevation  angle  for 
which  to  calculate  the  horizontal  radiation 
pattern  (from  the  angle  index  table) 


Wavelengths  and  Frequencies 


X,  .  41 

Xg  21 

u 

II 

uj| 

X4  :• 

f'  ■  k 

U 

M 

«ri» 

f4  = 

Possible  Operating  Frequencies  (Hertz) 

X  , 

-  Dipole: 

4 

f,  =  5-10* 

^  Dipole: 

f2  = 

3X  TV  . 

— —  Dipole: 

4 

f3  =  1.5*  107 

X  Dipole: 

f4  = 

180 


Input  the  operating  j  jq-  10^ 
Frequency  (Hertz)  5 


X5  •  j- 


f  - 


2 •* 


Input  the  Height  of  ^  =76  lflPut  the  ground 

Antenna  Feed  (meters)  0  '  Dielectric  Constant 


ct  ■=  30 


Input  the  Distance  ^  ..  jqqq  Input  the  ground  j.jq- 

from  Antenna  (meters)  Conductivity:  " 


Index  of  Refraction 


Rdj  =  R  -  Hq-cos(0j) 
Rr,  =  R  +  Hq-cos^j^ 


n  = 


.  18000-o 

er  -  j- 


(f510'6) 


Complex  Numerical  Distance 
for  Vertical  Polarization 


Pe, 


-j/3Rrj 
2-sin^«?jj  1 


cos 


W)  * 1 


n2  -  singly 


Complex  Numerical  Distance  D_  j  -  Rrj  f~^  /  ,  'j 

for  Horizontal  Polarization  *  -  .  /«\2  [  cos(  l)  +  Ja  “  \sul(  ji 


•sin^i)' 


i)2) 


Vertical  Reflection  Coefficient  T Vj 


Horizontal  Reflection  Coefficient 


(n2-cos(0,))  -  ( 

J”2  -  sta(*l)2) 

(n2  cos(0,))  «.  ( 

J“2  -  sin(®l)2) 

h ..  cos(9i)  -  (• 

n2  -  sin(0()2j 

cos(9|)  *  (. 

»2  -  5in(9|)2) 

181 


Radiation  Patterns  in  Phi=Pi /2  Plane  (Perpendicular  to  Dipole) 
Space  Wave  Radiation  Pattern  (Pfci=Pi/2) 


Max  E-Field  Intensity  (Volts  per  meter)  max(MagElP)  =  0.03632 
Surface  Wave  Radiation  Pattern  (Phi=pi/2) 


-1  -0.3  0  0.3  1 

,  x2Pj 


Max  E-Field  Intensity  (Volts  per  meter)  max(MagE2P)  =  1.03856*10  ® 
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Combined  Space  and  Surface  Wave  Radiation  Pattern  (Phi=pi/2) 
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/ 

/ 

/ 

/ 

/' 


Radiation  Patterns  in  Phi=0  Plane  (Parallel  to  Dipole) 
Space  Wave  Radiation  Pattern  (Phi=0) 


-1  -0.5  0  0.5  1 


HreC^.  HtOj 

Max  E-Field  Intensity  (Vohs  per  meter)  max(MagElO)  =  0.03594 
Surface  Wave  Radiation  Pattern  (Phi=0) 


-1  -0.5  0  0.5  1 

,  m20j 


Max  E-Field  Intensity  (Volts  per  meter)  max(MagE20)  =  3.84453-10'* 
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Combined  Space  and  Surface  Wave  Radiation  Pattern  (Phi=0) 
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Total  Power  Radiated  (Watts)  power  =  37.80286 


Radiation  Resistance  (versus 
Maximum  Antenna  Current) 
(Ohms) 


Radres  =  75.60572 


Directivity  (or  Maximum  Power  Gain  assuming  100%  Antenna  Efficiency) 

Phi  -  pi/2  Plane  Phi  =  0  Plane 

DirectivityP  =  5.23417  DirectivityO  =  5.12535 

Effective  Isotropic  Radiated  Power  (EIRP)  (W atts) 

Phi  -  pi/2  Plane  Phi  =  0  Plane 

DirectivityP- power  =  197.86674  DirectivityO- power  =  193.75275 


Maximum  Effective  Area  (Along  Radial  of  Directivity) 
(square  meters) 


Phi  =  pi/2  Plane 


Phi  =  0  Plane 


(Aj)  ^  DirectivityP 

4-x 


=  374.87006 


(Xj)  ^  PTectivityO 
4-x 


=  367.07586 
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Maximum  Effective  Length  (Along  Radial  of  Directivity) 

(meters) 

Phi  =  pi/2  Plane 

Phi  =  0  Plane 

Cadres-  2*  DirectivhyP 

Radres-  (  Xe') 2  DirectivityO 
17  34133  2-  v  ' 

480  *2 

,  480-x2 

Numerical  Distances 

Vertical  Polarization 

Horizontal  Polarization 

|PeQ|  =  5.04617 

|PmQ|  =  1.92562- 104 

Elevation  Angle  of  Directivity  (Maximum  Gain) 
above  the  Horizon  (Degrees) 

Phi  =  pi/2  Plane 

Phi  =  0  Plane 

AngleP=  (  65,895f) 

\  65.8954/ 

A»gl,0=(  89  85737  ) 

\  89.85737/ 

HORIZONTAL  DIPOLE 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated  with 
horizontal  thin-wire  dipole  antennas  (diameter  «  wavelength).  The  antenna  is 
mounted  above  and  parallel  to  the  x-axis  in  a  rectangular  coordinate  system.  The 
feed  is  at  the  center  of  the  antenna  at  a  set  height  directly  above  the  origin. 

Required  inputs  ase  the  antenna  length,  feed  height  above  the  surface,  transmitted 
frequency,  distance  from  the  antenna,  the  conductivity  aad  dielectric  constant  of  the 
surface  below  the  antenna.  The  planar  earth  model  is  assumed  in  predicting 
radiation  patterns.  Predicted  operating  frequencies  assume  that  the  antenna  is  a 
quarter-wavelength,  balf-wavelength.  three-quarter-wavelength,  or  full  wave¬ 
length  dipole.  A  sinusoidal  current  input  with  a  maximum  of  unity  is  assumed.  All 
radiation  patterns  are  normalized  with  respect  to  the  maximum  electric  field 
intensity  transmitted  by  the  antenna  in  the  plane  of  interest  The  electric  field 
magnitudes  to  which  the  patterns  are  normalized  are  displayed  below  their 
respective  plots.  The  radiation  patterns  are  plotted  for  the  phi=0  and  phi=pi/2 
vertfc*'  planes  parallel  and  perpendicular  to  the  x-axis  respectively.  A  horizontal 
radir .  n  pattern  is  plotted  at  an  elevation  selected  by  the  index  from  the  elevation 
ang  e  index  table.  Polarization  is  a  combination  of  vertical  and  horizontal 
depending  upon  spatial  orientation  relative  to  the  dipole. 


Input  the  dipole 
length  in  meters 


Input  the  index  of  the  elevation  angle  for 
which  to  calculate  the  horizontal  radiation 
pattern  (from  the  angle  index  table) 

Wavelengths  and  Frequencies 


Xt  .  4  1 

X2  2  1 

X  4,1 

h  *  t 

X^  =  1 

r  c 

f1  \ 

A1 

<2  *  t 
a2 

h  ■  r 

A3 

*  i 

Possible  Operating  Frequencies  (Hertz) 

7  Dipole: 

r,  =  i-io7 

^  Dipole: 

f2  *  2- I07 

3X  IV  . 

—  Dipole: 

4 

f,  =  3-  to7 

X  Dipole: 

r4  =  4* io7 
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Input  the  operating 
Frequency  (Hertz) 


Input  the  Height  of  U  i  76  ground 

Antenna  Feed  (meters)  ®  Dielectric  Constant 


Input  the  Distance  ^  ^qqq  Input  the  ground  g  _ 
from  Antenna  (meters)  Conductivity: 


Index  of  Refraction 


Complex  Numerical  Distance 
for  Vertical  Polarization 


Complex  Numerical  Distance 
for  Horizontal  Polarization 


Vertical  Reflection  Coefficient 


Horizontal  Reflection  Coefficient 
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Radiation  Patterns  in  Phi=Pi/2  Plane  (Perpendicular  to  Dipole) 
Space  Wave  Radiation  Pattern  (Phi=Pi/2) 


-1  -0.5  0  0.5  1 

ureL  ,«tP. 

Max  E-Field  Intensity  (Voltr.  per  meter)  max(MagElP)  =  0.03737 
Surface  Wave  Radiation  Pattern  (Phi=pi/2) 


h2Pj 

Max  E-Field  Intensity  (Vchs  per  meter)  max(MagE2P)  =  8.03309*10 
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Radiation  Patterns  in  Horizontal  Plane 
Space  Wave  Radiation  Pattern 


-1  -0.3  0  C.3  1 


,  HtE, 

Max  E-Field  Intensity  (Volts  per  meter)  max(MagElH)  =  0.03734 
Surface  Wave  Radiation  Pattern  (Horizontal  Plane) 


-1  -0.3  0  0.5  1 

we(^ ,  h2Hj 


Max  E-Field  Intensity  (Volts  per  meter)  max  ( MagE2H )  =  1 .62046*  10 


Total  Power  Radiated  (Watts) 


power  =  27.68975 


Radiation  Resistance  (versus 
Maximum  Antenna  Current) 
(Ohms) 


Radres  =  55.3795 


Directivity  (or  Maximum  Power  Gain  assuming  100%  Antenna  Efficiency) 

Phi  =  pi/2  Plane  Phi  -  0  Plane 

Directi  vityP  =  7.56408  DireclivityO  =  1.32021 

Effective  Isotropic  Radiated  Power  (EIRP)  (Watts) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 

DirectivityP- power  =  209.44744  DirectivityO- power  =  36.55637 
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Maximum  Effective  Length  (Along  Radial  of  Directivity) 

(meters) 

Phi  -  pi/2  Plane  Phi  *  0  Plane 


Radres-  (Aj)  2  DirectivilyP 
j  480- *2 


=  8.92079 


Radres*  2 •  DirectivhyO 


480*x2 


Numerical  Distances 


Vertical  Polarization  Horizontal  Polarization 

|  PeQ|  =  15.283  |  Pm0|  =  2.48%*  104 


Elevation  Angle  of  Directivity  (Maximum  Gain) 
above  (he  Horizon  (Degrees) 

Phi  =  pi/2  Plane  Phi -0  Plane 


AngteP  = 


/  28.24089  \ 
\  28.24089/ 


AngleO  = 


/  41.64818  \ 
\  41.64818/ 


3.7269 
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HORIZONTAL  DIPOLE 

This  application  calculates  far  field  radiation  patterns  and  parameters  associated  with 
horizontal  thin-wire  dipole  antennas  (diameter  «  wavelength).  The  antenna  is 
mounted  above  and  parallel  to  the  x-axis  in  a  rectangular  coordinate  system.  The 
feed  is  at  the  center  of  the  antenna  at  a  set  height  directly  above  the  origin. 

Required  inputs  are  the  antenna  length,  feed  height  above  the  surface,  transmitted 
frequency,  distance  from  (he  antenna,  the  conductivity  and  dielectric  constant  of  the 
surface  below  the  antenna.  The  planar  earth  model  is  assumed  in  predicting 
radiation  patterns.  Predicted  operating  frequencies  assume  that  the  antenna  is  a 
quarter-wavelength,  half-wavelength,  three-quarter-wavelength,  or  full  wave¬ 
length  dipole.  A  sinusoidal  current  input  with  a  maximum  of  unity  is  assumed.  Ail 
radiation  patterns  are  normalized  with  respect  to  the  maximum  electric  field 
intensity  transmitted  by  the  antenna  in  the  plane  of  interest.  The  electric  field 
magnitudes  to  which  the  patterns  are  normalized  are  displayed  below  their 
respective  plots.  The  radiation  palterns  are  plotted  for  the  phi=0  and  phi=pi/2 
vertical  planes  parallel  and  perpendicular  to  the  x-axis  respectively.  A  horizontal 
radiation  pattern  is  plotted  at  an  elevation  selected  by  the  index  from  the  elevation 
angle  index  table.  Polarization  is  a  combination  of  vertical  and  horizontal 
depending  upon  spatial  orientation  relative  to  the  dipole. 

Input  the  dipole  1*5  h  -  - 

length  in  meters  /  2 

Input  the  index  of  the  elevation  angle  for 

v/hich  to  calculate  the  horizontal  radiation  d  =  535 

pattern  (from  the  angle  index  table) 


Wavelengths  and  Frequencies 


X,  =  41 

X2  =  2  1 

x3 

41 

3 

x4  ,  1 

f<  =  f, 

h  *  f 

a2 

f3 

c 

h 

•■n 

ii 

-ri" 

Possible  Operating  Frequencies  (Hertz) 

A  , 

-  Dipole: 

4 

ft  =  1.5- 107 

X 

2 

Dipole: 

f2  =  3- 10 

3X  TV  . 

— —  Dipole: 

4 

f3  =  4.5- 107 

X 

Dipole: 

f4  =  6- 10 

197 


Input  the  operating 
Frequency  (Hertz) 


f5  =  30-  10S 


Input  ihe  Height  of  ^  -  76  ^P0*  the  ground 

Antenna  F'ed  (meters)  0  '  ‘  Dielectric  Constant 


Input  the  Distance 
from  Antenna  (meters) 


R  ■-  3000  Ittpul  the  ground 
Conductivity: 


o  =  3* I0~  2 


Index  of  Refraction 


Rdj  =  R  -  Hq-cos^j) 
RTj  :a  R  +  Hg-  COS^jj 


Complex  Numerical  Distance 
for  Vertical  Polarization 


Complex  Numerical  Distance  p 
for  Horizontal  Polarization 


Vertical  Reflection  Coefficient 


Horizontal  Reflection  Coefficient 


2-sin  2 


r.,  ..  (n2<COS(9l))  -  (j 

fn2  -  sin ( dyj  2 j 

(«*•«»(»,))  *  (j 

n2  -  sin (*«)*) 

rh  ..  CM(*i)  -  (i" 

|2  -  sta(5i)2) 

2 


Radiation  Patterns  is  Phi=0  Plane  (Parallel  to  Dipole) 
Space  Wave  Radiation  Pattern  (Phi=0) 


Hrei^ ,  m(0| 

Max  E-Field  Intensity  (Vohs  per  meter)  max(MagElO)  =*  0.03422 
Surface  Wave  Radiation  Pattern  (Phi=0) 


-1  -0.3  0  0.3  I 

Hrti^ ,  m20| 

Max  E-Field  Intensity  (Volts  per  meter)  max(  MagE20)  =>  9.38091  •  10  * 
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Radiation  Patterns  in  Horizontal  Plane 
Space  Wave  Radiation  Pattern 


Max  E-Field  Intensity  (Volts  per  meter)  max(MagElH)  =  0.02972 
Surface  Wave  Radiation  Pattern  (Horizontal  Plane) 


-1  -0.5  0  0.5  1 

wet^ ,  k2H| 


Max  E-Fiek)  Intensity  (Vohs  per  meter)  max(MagE2H)  =  1.14409*  10  * 
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Total  Power  Radiated  (Watts) 


power  =  31.69054 


Radiation  Resistance  (versus 
Maximum  Antenna  Current) 
(Ohms) 


Radres  =  63.38107 


Directivity  (or  Maximum  Power  Gain  assuming  10 0%  Antenna  Efficiency) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 

DirectivityP  ■  6.8304  DirectivityO  *  5.54301 

Effective  Isotropic  Radiated  Power  (EIRP)  (Watts) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 

DirectivityP- power  =  216.45893  DirectivityO- power  *  175.66087 


Maximum  Effective  Area  (Along  Radial  of  Directivity) 
(square  meters) 


Phi  =  pi/2  Plane 


Phi  *  0  Plane 


(X5)  2-  DirectivityP  ^  ^  ( X5)  2-  DirectivityO 

“  .  s  54.35457  —  s 

4-*  4-* 


*  44.10985 
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Maximum  Effective  Length  (Along  Radial  of  Directivity) 

(meters) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 


Radres-  (  Ac)  •  DirectivityP 

v  ^ —  -  =  6.04592  2- 


480-x2 


Radres-  2-  DirectivhyO 


480- x* 


Numerical  Distances 

Vertical  Polarization  Horizontal  Polarization 


|  PeQ  |  =  26.28165  |  Pm0 1  =  3.21687- 104 


Elevation  Angte  of  Directivity  (Maximum  Gain) 
above  the  Horizon  (Degrees) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 


AngleP  = 


/  18.542  \ 
\  18.542/ 


AngleO  =* 


89.85737  \ 
89.85737/ 


5.44643 
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HORIZONTAL  DIPOLE 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated  with 
horizontal  Urn-wire  dipole  antennas  (diameter  «  wavelength).  The  antenna  is 
mounted  above  and  parallel  to  the  x-axis  in  a  rectangular  coordinate  system.  The 
feed  is  at  the  center  of  the  antenna  at  a  set  height  directly  above  the  origin. 

Required  inputs  are  the  antenna  length,  feed  height  above  the  surface,  transmitted 
frequency,  distance  from  the  antenna,  the  conductivity  and  dielectric  constant  of  the 
surface  below  the  antenna.  The  planar  earth  model  is  assumed  in  predicting 
radiation  patterns.  Predicted  operating  frequencies  assume  that  the  antenna  is  a 
quarter-wavelength,  half- wavelength,  three-quarter-wavelength,  or  full  wave¬ 
length  dipole.  A  sinusoidal  current  input  with  a  maximum  of  unity  is  assumed.  All 
radiation  patterns  are  normalized  with  respect  to  the  maximum  electric  field 
intensity  transmitted  by  the  antenna  in  the  plane  of  interest  The  electric  field 
magnitudes  to  which  the  patterns  are  normalized  are  displayed  below  their 
respective  plots.  The  radiation  patterns  are  plotted  for  the  phi=0  and  phi=pi/2 
vertical  planes  parallel  and  perpendicular  to  the  x-axis  respectively.  A  horizontal 
radiation  pattern  is  plotted  at  an  elevation  selected  by  the  index  from  the  elevation 
angle  index  table.  Polarization  is  a  combination  of  vertical  and  horizontal 
depending  upon  spatial  orientation  relative  to  the  dipole. 

Input  the  dipole  1  =  3  h  =  - 

length  in  meters  2 

Input  the  index  of  the  elevation  angle  for 

which  to  calculate  the  horizontal  radiation  d  *  535 

pattern  (from  the  angle  index  table) 

Wavelengths  and  Frequencies 

X,  41  Xj  .  2  1  X,  »  ±!  x4  •  I 


Possible  Operating  Frequencies  (Hertz) 

^  Dipole:  ft  =  2.5- 107  ^  Dipole:  f2  =  5- 107 

^Dipole:  f3=7.5-107  A  Dipole:  f 4  =  1  •  108 
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Input  the  operating  f  ..  50. 1Q6 
Frequency  (Hertz)  5 


Input  the  Height  of 
Antenna  Feed  (meters) 


j  5  Input  the  ground 
Dielectric  Constant 


a  :*  72 


Input  the  Distance 
from  Antenna  (meters) 


R  -  3000 


Input  the  ground 
Conductivity: 


o  =  4 


Index  of  Refraction 


Rdj  =  R  -  Hq-cos(0j) 
Rtj  =  R  +  Hq-cos^jJ 


,  18000- o 


Complex  Numerical  Distance 
for  Vertical  Polarization 


Pe, 


-  j  /3  Rr, 
2-sin(0j)* 


pnc  /  A  \  i.  ^ 

-  sin^0|^  j 

Ufa  I  i/|  1  +  - 

»2  / 

Complex  Numerical  Distance 
for  Horizontal  Polarization 


Vertical  Reflection  Coefficient 


Horizontal  Reflection  Coefficient 
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Radiation  Patterns  in  Phi=Pi/2  Plane  (Perpendicular  to  Dipole) 
Space  Wave  Radiation  Pattern  (Phi=Pi/2) 


Max  E-Field  Intensity  (Volts  per  meter)  max (  MagEl  P )  =  0.03926 


Surface  Wave  Radiation  Pattern  (Phi=pi/2) 


-1  -0.3  0  0.3  1 

.  k2Pj 


Max  E-Field  Intensity  (Volts  per  meter)  max(  MagE2P)  =  8.83121  •  10 


;■* ?? =Vj H‘As?:V;Y^'  •:&%?££*  r: £?$!£  ■ -Vi •? ^ 


Radiation  Patterns  in  Horizontal  Plane 
Space  Wave  Radiation  Pattern 


-1  -0.5  0  0.5  1 


HTtt^ .  HtE| 

Max  E-Field  Intensity  (Volts  per  meter)  max(  MagElH)  =  0.02647 
Surface  Wave  Radiation  Pattern  (Horizontal  Plane) 


-1  -0.5  0  0.5  1 

wef^,  h2Hj 

Max  E-Field  Intensity  (Volts  per  meter)  max(MagE2H)  =  4.07737*10  ® 
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Combined  Space  and  Surface  Wave 


""V 


Max  E-Field  Intensity  (Volts  per 


Total  Power  Radiated  (Watts) 


power  -  41.98012 


Radiation  Resistance  (versus 
Maximum  Antenna  Current) 
(Ohms) 


Radres  =  83.96025 


Directivity  (or  Maximum  Power  Gain  assuming  100%  Antenna  Efficiency) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 

DirectivityP  =  5.50819  DirectivityO  =  5.50321 

Effective  Isotropic  Radiated  Power  (EIRP)  (Watts) 

Phi  =  pi/2  Plans  Phi  =  0  Plane 

DirectivityP- power  =  231.23457  DirectivityO-  poiver  =  231.02563 


Maximum  Effective  Area  (Along  Radial  of  Directivity) 
(square  meters) 


Phi  =  pi/2  Plane 


Phi  =  0  Plane 


(Aj)  2-  DirectivityP 
4-x 


15.77981 


(Aj)  DirectivityO 
4-x 


15.76555 
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Maximum  Effective  Length  (Along  Redial  of  Directivity) 

(meters) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 


Radrcs-  (  X«)  2  DirectivityP  Radres-  (  X.)  2-  DirectivHyO 

2.  - - 1_  »  3.74932  2-  - ^ -  =  3.74762 

«  480  x2  ,  480  x2 


Numerical  Distances 

Vertical  Polarization  Horizontal  Polarization 

|PeQ|  *  1.08943  |  PmQ  |  =  2.26469- 10S 

Elevation  Angle  of  Directivity  (Maximum  Gain) 
above  tbe  Horizon  (Degrees) 

Phi  *  pi/2  Plane  Phi  *  0  Plane 


AngteO  = 


/  89.85737  \ 
\  89.85737/ 


213 


APPENDIX  Pi 

VERTICAL  LOG- PERIODIC  DIPOLE  ARRAY  COMPUTER  OPT PUT 

This  appendix  contains  computer  hardcopies  from  the 
Mathcad  vertical  log-periodic  dipole  array  application  which 
show  the  input  values  and  predicted  radiation  characteristics 
for  two  sample  calculations.  The  configuration  is  a  given  by 
the  inputs  on  the  first  two  pages  of  each  printout.  The  first 
antenna  is  mounted  above  soil  (ef=10  and  o^lO'3)  and  the  second 
above  seawater  (er=72  and  a=«4)  .  Reference  6  [p.  114]  provides 
the  radiation  patterns  and  gain  predictions  from  several 
sources  for  the  configuration  in  the  first  example. 

The  radiation  patterns  and  maximum  directive  gain 
computed  by  the  first  Mathcad  example  are  almost  identical  to 
the  those  given  in  reference  6.  The  directivity  prediction  is 
slightly  higher  for  Mathcad,  and  the  elevation  of  the  maximum 
directive  gain  is  also  slightly  higher,  but  the  overall 
similarity  between  the  predicted  radiation  characteristics  is 
still  quite  good.  The  seawater  example  also  yields  results 
consistent  with  expectations.  With  respect  to  the  soil 
example,  the  seawater  example's  directivity  is  higher  due  to 
a  stronger  reflected  wave,  and  the  surface  wave  is  stronger  at 
grazing  angles  (0-90°)  due  to  higher  conductivity. 
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VERTICAL  LOG  PERIODIC  DIPOLE  ARRAY 

This  application  calculates  far  field  radiation  patterns  and  parameters  associated  with 
vertical  log  periodic  dipole  arrays.  The  antenna  is  oriented  such  that  the  projection 
of  its  center  axis  lies  on  the  positive  y-axis  of  a  rectangular  coordinate  system.  The 
antenna  axis  can  be  oriented  with  respect  to  the  vertical  at  any  angle  between  zero 
and  ninety  degrees.  The  feed  is  at  the  center  of  the  shortest  element  directly  above 
the  origin.  The  dipole  elements  are  bisected  by  the  antenna  axis  and  are  parallel  to 
the  z-axis.  Required  inputs  are  the  number  of  dipole  elements,  shortest  and  second 
shortest  element  lengths,  separation  between  the  shortest  and  second  shortest 
elements,  radius  of  the  shortest  element,  height  of  the  feed  above  the  surface, 
angle  of  the  antenna  axis  with  the  vertical,  characteristic  admittance  of  the  line 
feeding  the  antenna,  termination  impedance  of  the  antenna,  transmitted  frequency, 
distance  from  the  antenna,  the  conductivity  and  dielectric  constant  of  the  surface 
below  the  antenna.  The  planar  earth  model  is  assumed  in  predicting  radiation 
patterns.  A  sinusoidal  current  input  with  a  maximum  of  unity  is  assumed.  All 
radiation  patterns  are  normalized  with  respect  to  the  maximum  electric  field 
intensity  transmitted  by  the  antenna  in  the  plane  of  interest.  The  electric  field 
magnitudes  to  which  the  patterns  are  normalized  are  displayed  below  their 
respective  plots.  Radiation  patterns  are  plotted  for  the  phi=0  and  phi=pi/2  vertical 
planes  perpendicular  and  parallel  to  the  y-axis  respectively.  A  horizontal  radiation 
pattern  is  plotted  at  an  elevrtion  selected  by  the  index  from  the  elevation  angle 
index  table.  Polarization  is  vertical  for  all  vertical  log-periodic  antennas. 


Input  the  number 
of  elements 


Input  the  length 
of  the  shortest 
element  (meters) 

Input  the  distance 
from  shortest  to 
second  shortest 
element  (meters)  1 

Input  the  operating 
Frequency  (Hertz) 

Input  the  ground 
Conductivity 


N  •  12  Input  the  sKorlesI  rad  ,  .009,8 

element  s  Radius  0 

(meters) 

Input  the  length 

Iq  =  3.673  of  the  second  1  4.373 

shortest  element 
(meters) 

Input  the  height  of 

dg  =  1.0766  the  shortest  element  Hq  :=  6.87 

above  (he  surface 
(meters) 


f5  =  18  10* 


o  =  MO' 


Input  the  Characteristic 
Admittance  of  the  Line  ADM 
Feeding  the  Antenna 
(Ohms) 


Input  the  Distance  r  ?  > 
from  the  Antenna 
(meters) 

Input  the  ground 
Dielectric  Constant 

Input  the  Termination 
Impedance  connected  TIMP 
to  longest  element 
(Ohms) 


7  V  .  1  s .. .  <  -  \//: 


,4. 


Input  the  angle  between 
the  vertical  (z  axis)  and  y  -  75 
the  antenna  axis  (degrees) 

NOTE:  The  Vertical  Log-Periodic  Dipole  mutual  impedance  calculations  are  not 
valid  if  there  is  an  antenna  element  whose  length  is  an  exact  integer  multiple  of  the 
wavelength.  If  (his  occurs,  there  will  be  a  singularity  error  in  the  mutual  impedance 
calculations.  If  this  problem  arises,  K  will  be  necessary  to  vary  the  operating  frequency 
such  that  no  element  is  exactly  an  integer  multiple  of  the  wavelength. 

If  the  log-periodic  calculations  of  element  length  and  spacing  do  not  represent  the 
desired  antenna  configuration,  you  can  enter  the  values  directly  by  selecting  the 
variables  1  and  d,  using  the  define  key  (shift,  colon),  and  entering  the  lengths  and 
spacings  by  separating  each  successive  entry  by  a  comma. 

Calculated  distance 

Calculated  length  between  successive  Calculated  radii 

of  elements  from  elements- shortest  of  elements  from 

shortest  to  longest  to  longest  shortest  to  longest 

d 

v _ 

1.0766 
1.28178 
1.52606 
1.8169 
2.16316 
2.57541 
3.06624 
3.6506 
4.34633 
5.17465 
6.16084 


I 

1 


3.673 

4.373 

5.20641 

6.19864 

7.37998 

8.78645 

10.46098 

12.45463 

14.82823 

17.65419 

21.01873 

25.02447 

rad 


I 


0.00918 

0.01093 

0.01301 

0.01549 

0.01844 

0.021% 

0.02615 

0.03113 

0.03706 

0.04412 

0.05253 

0.06254 

Input  the  elevation  index 

for  which  to  calculate  the  w  r  285 

horizontal  radiation  pattern 

(from  angle  index  table) 


Radiation  Patterns  in  Phi=0  Plane 
(Perpendicular  to  Phi=pi/2  Plane) 


we£^.  xIO. 


Max  E-Field  Intensity  (Volts  per  meter)  max(MagElO)  =  0.03221 
Surface  Wave  Radiation  Pattern  (Phi=0) 


-1  -0.5  0  0.5  1 

Kfet^ ,  k20j 

Max  E-Field  Intensity  (Volts  per  meter)  max(MagE20)  =  7.64667-10 
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Radiation  Patterns  in  Horizontal  Plane 
(Parallel  to  Ground) 

Space  Wave  Radiation  Pattern  (Horizontal) 


Max  E-Field  Intensity  (Volts  per  meter)  max(MagElH)  =  0.05438 
Surface  Wave  Radiation  Pattern  (Horizontal) 


Max  E-Field  Intensity  (Volts  per  meter)  max(MagE2H)  =  1.09861*  10  * 
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/ 


N, 


Y- 


powerfl 


30 

*-RJ 


* 

2 


^  JO 


j'P-Y  wc(a2  t  a3)  (coj(O  si»(a2  +  a3)  +  srr»( ^ ) 
1  -e  7 
r 


power  '=  power# 


Radres  -  2 -power 


sqE3Pj 


^*^3F|!;2  DireclivilyP  .  lY **  ™<i3gg> 

2(120-*)  }  povycr 


sqE30j  :» 


( |  MagE30|  | ) ' 


2-(  120-*) 


Total  Power  Radiated  (Watts) 


Radiation  Resistance  (versus 
Maximum  Input  Current) 
(Ohms) 


DirectivityO  4-*  R  max(stlE30> 

power 


power  =  34.29572 


Radres  =  68.59144 


Directivity  (or  Maximum  Power  Gain  assuming  100%  Antenna  Efficiency) 

Phi  *  pi/2  Plane  Phi  =  0  Plane 

Directi  vityP  =  13.63346  DirectivityO  =  4.50794 

Effective  Isotropic  Radiated  Power  (EIRP)  (Watts) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 

DirectivityP- power  =  467.56917  DirectivityO- power  =  154.6031 
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vV-Y 


Maximum  Effective  Area  (Along  Radial  of  Directivity) 
(square  meters) 


Phi  =  pi/2  Plane 

(A*)  2  DirectivityP 

-  =  301.36553 

4-* 


Phi  =  0  Plane 

(  Aj)  2-  DirectivityO 

A-k 


99.64739 


Maximum  Effective  Length  (Along  Radial  of  Directivity) 

(meters) 

Phi  =  pi/2  Plane  Phi  *  0  Plane 


Radres-  (  A •  DirectivityP 
-  -  14.8097  2- 


48C-x‘ 


Radres-  (Aj)  2-  DirectivityO 


480-x4 


Numerical  Distance  for  Vertical  Polarization 
|ft0|  -50.7 

Elevation  Angle  of  Directivity  (Maximum  Gain) 
above  the  Horizon  (Degrees) 


Phi  =  pi/2  Plane 


Phi  =  0  Plane 


AngleP  «  14.28571 


AngleO  = 


1 3.71 429  \ 
13.71429/ 


8.51594 
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VERTICAL  LOG  PERIODIC  DIPOLE  ARRAY 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated  with 
vertical  log  periodic  dipole  arrays.  The  antenna  is  oriented  such  that  the  projection 
of  its  center  axis  lies  on  the  positive  y-axis  of  a  rectangular  coordinate  system.  The 
antenna  axis  can  be  oriented  with  respect  to  the  vertical  at  any  angle  between  zero 
and  ninety  degrees.  The  feed  is  at  the  center  of  the  shortest  element  directly  above 
the  origin.  The  dipole  elements  are  bisected  by  the  antenna  axis  and  are  parallel  to 
the  z-axis.  Required  inputs  are  the  number  of  dipole  elements,  shortest  and  second 
shortest  element  lengths,  separation  between  the  shortest  and  second  shortest 
elements,  radius  of  the  shortest  element,  height  of  the  feed  above  the  surface, 
angle  of  the  antenna  axis  with  the  vertical,  characteristic  admittance  of  the  line 
feeding  the  antenna,  termination  impedance  of  the  antenna,  transmitted  frequency, 
distance  from  the  antenna,  the  conductivity  and  dielectric  constant  of  the  surface 
below  the  antenna.  The  planar  earth  model  is  assumed  in  predicting  radiation 
patterns.  A  sinusoidal  current  input  with  a  maximum  of  unity  is  assumed.  All 
radiation  patterns  are  normalized  with  respect  to  the  maxinium  electric  field 
intensity  transmitted  by  the  antenna  in  the  plane  of  interest.  The  electric  field 
magnitudes  to  which  the  patterns  are  normalized  are  displayed  below  then- 
respective  plots.  Radiation  patterns  are  plotted  for  the  phi=0  and  phi=pi/2  vertical 
planes  perpendicular  and  parallel  to  the  y-axis  respectively.  A  horizontal  radiation 
pattern  is  plotted  at  an  elevation  selected  by  the  index  from  the  elevation  angle 
index  table.  Polarization  is  vertical  for  all  vertical  log-periodic  antennas. 


Input  the  number 
of  elements 


Input  the  length 
of  the  shortest 
element  (meters) 

Input  the  distance 
from  shortest  to 
second  shortest 
element  (meters) 

Input  the  operating 
Frequency  (Hertz) 

Input  the  ground 
Conductivity 


N  ■=  12 


In  :«  3.673 


dfl  =  1.0766 


f.  :»  18- 10' 


o  =  4 


Input  the  shortest 
element's  Radius 
(meters) 

Input  the  length 
of  the  second 
shortest  element 
(meters) 

Input  the  height  of 
the  shortest  element 
above  the  surface 
(meters) 

Input  the  Distance 
from  the  Antenna 
(meters) 

Input  the  ground 
Dielectric  Constant 


radQ  :=  .00918 


<1  - 


4.373 


Hft  --  6.87 


3000 


«r  :«  72 


Input  the  Characteristic 
Admittance  of  the  Line  ADM 
Feeding  the  Antenna 
(Ohms) 


1 

450 


Inpul  the  Termination 
Impedance  connected  T7MP 
to  longest  element 
(Ohms) 


=  450 
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Input  the  angle  between  Input  the  elevation  index 

the  vertical  (z  axis)  and  ^  73  for  which  to  calculate  the  w  r  285 

the  antenna  axis  (degrees)  horizontal  radiation  pattern 

(from  angle  index  table) 

NOTE:  The  Vertical  Log-Periodic  Dipole  mutual  impedance  calculations  are  not 
valid  if  there  is  an  antenna  element  whose  length  is  an  exact  integer  multiple  of  the 
wavelength.  If  this  occurs,  there  will  be  a  singularity  error  in  the  mutual  impedance 
calculations.  If  this  problem  arises,  it  will  be  necessary  to  vary  the  operating  frequency 
such  that  no  element  is  exactly  an  integer  multiple  of  the  wavelength. 


If  the  log-periodic  calculations  of  element  length  and  spacing  do  not  represent  the 
desired  antenna  configuration,  you  can  enter  the  values  directly  by  selecting  the 
variables  1  and  d,  using  the  define  key  (shift,  colon),  and  entering  the  lengths  and 
spacings  by  separating  each  successive  entry  by  a  comma. 

Calculated  distance 

Calculated  length  between  successive  Calculated  radii 

of  elements  from  elements-shortest  of  elements  from 

shortest  to  longest  to  longest  shortest  to  longest 


I  d  rad 
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Radiation  Patterns  in  Horizontal  Plane 
(Parallel  to  Ground) 

Space  Wave  Radiation  Pattern  (Horizontal) 


Phi=0 


Max  E-Field  Intensity  (Vohs  per  meter)  max(MagEIH)  =  0.06821 
Surface  Wave  Radiation  Pattern  (Horizontal) 


xref^ ,  h2H. 

Max  E-Field  Intensity  (Volts  per  meter)  max(  MagE2H)  =  4.93203*  10 

231 


power#  = 


r  •*  r  •w 

I<«  I  <• 

2  2 


\  Jo 
2 


jp-Y  Jec(ot2  +  a3)  (cos(£)  sh(a2  +  a3)  +-  sin(^) 

X-  Ve  Y 


power  power# 


Radres  =  2- power 


sqE3Pj 


(j  MagE3Pj|): 
2*(  l20-%) 


DirectivityP  ,  4  ’R  m,x< 

power 


sqE30j 


( |  MagE30j  | ) ' 
2-(120-x) 


DiredivityO  = 


4-*-R  -max(sqE30) 
power 


Total  Power  Radiated  (Watts) 


power  =  99.48703 


Radiation  Resistance  (versus 
Maximum  Input  Current) 
(Ohms) 


Radres  =  198.97405 


Directivity  (or  Maximum  Power  Gain  assuming  100%  Antenna  Efficiency) 


Phi  =  pi/2  Plane 
DirectivityP  =  23.15454 


Phi  =  0  Plane 
DbectivityO  =  9.2377 


Effective  Isotropic  Radiated  Power  (EIRP)  (Watts) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 

DirectivityP- power  =  2.30358- 10^  DiredivityO- power  =  919.030% 


•N  ’  Vr  r  /A'.  /  ■ 


'  ■  -  /<f f- 


Maximum  Effective  Area  (Along  Radial  of  Directivity) 
(square  meters) 


Phi  =  pi/2  Plane 


Phi  =  0  Plane 


(  Xe)  •  DirectivhyP 

-  =  511.82763 

4* 


(Xj)  ^  Directi  vityO 


4  •% 


204.19793 


Maximum  Effective  Length  (Along  Radial  of  Directivity) 

(meters) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 


Radres-  (  X.)  ^  DirectivityP 
- LJz - - -  =  32.87188  2- 


480  *4 


Radres-  ( Xj)  2*  DirectivityO 


480  x2 


=  20.76291 


Numerical  Distance  for  Vertical  Polarization 
|P*0|  =  0.14135 

Elevation  Angle  of  Directivity  (Maximum  Gain) 
above  the  Horizon  (Degrees) 


Phi  =  pi/2  Plane 


Phi  =  0  Plane 


AngieP  =  0 
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APPENDIX  E: 

HORIZONTAL  LOG- PERIODIC  DIPOLE  ARRAY  COMPUTER  OUTPUT 
This  appendix  contains  computer  hardcopies  from  the 
Mathcad  horizontal  log-periodic  dipole  array  application  which 
show  the  input  values  and  predicted  radiation  characteristics 
for  txtfo  sample  calculations.  The  configuration  is  a  given  by 
the  inputs  on  the  first  two  pages  of  each  printout.  The  first 
antenna  is  mounted  above  soil  (er=4  and  a=10'3)  and  the  second 
above  seawater  (er-72  and  cr=4)  .  Reference  6  [p.  108]  provides 
the  radiation  patterns  and  gain  predictions  from  several 
sources  for  the  configuration  in  the  first  example. 

The  radiation  patterns  and  maximum  directive  gain 
computed  by  the  first  Mathcad  example  are  almost  identical  to 
the  those  given  in  reference  6.  The  directivity  prediction  is 
slightly  higher  for  Mathcad,  but  the  overall  similarity  among 
the  predicted  radiation  patterns  is  remarkably  good.  The 
seawater  example  has  results  very  much  like  the  soil  example. 
For  horizontal  polarization,  the  higher  conductivity  surface 
below  the  antenna  does  not  result  in  a  greatly  enhanced 
surface  wave  and  increased  directivity  as  it  does  for  vertical 
polarization. 
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HORIZONTAL  LOG  PERIODIC  DIPOLE  ARRAY 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated  with 
horizontal  log  periodic  dipole  arrays.  The  antenna  is  oriented  such  that  the 
projection  of  its  center  axis  lies  on  the  positive  y-axis  of  a  rectangular  coordinate 
system.  The  antenna  axis  can  be  oriented  with  respect  to  the  vertical  at  any  angle 
between  zero  and  ninety  degrees.  The  feed  is  at  the  center  of  the  shortest  t  tement 
directly  above  the  origin.  The  dipole  elements  are  bisected  by  the  antenna  -  xis  and 
are  parallel  to  the  x-axis.  Required  inputs  are  the  number  of  dipole  elements, 
shortest  and  second  shortest  element  lengths,  separation  between  the  shortest  and 
second  shortest  elements,  radius  of  the  shortest  element,  height  of  the  feed  above 
the  surface,  angle  of  the  antenna  axis  with  the  vertical,  characteristic  admittance  of 
the  line  feeding  the  antenna,  termination  impedance  of  the  antenna,  transmitted 
frequency,  distance  from  the  antenna,  the  conductivity  and  dielectric  constant  of  the 
surface  below  the  antenna.  The  planar  earth  model  is  assumed  in  predicting 
radiation  patterns.  A  sinusoidal  current  input  with  a  maximum  of  unity  is  assumed . 
All  radiation  patterns  are  normalized  with  respect  to  the  maximum  electric  field 
intensity  transmitted  by  the  antenna  in  the  plane  of  interest.  The  electric  field 
magnitudes  to  which  the  patterns  are  normalized  are  displayed  below  their 
respective  plots.  Radiation  patterns  are  plotted  for  the  phi=pi/2  and  phi=0  vert,-,  at 
planes  parallel  and  perpendicular  to  the  y-axis  respectively.  A  horizontal  radial  on 
pattern  is  plotted  at  an  elevation  selected  by  the  index  from  the  elevation  angle 
index  table.  Polarization  is  a  varying  combination  of  horizontal  and  vertical 
depending  on  spatial  orientation  with  respect  to  the  antenna. 


Input  the  number 
of  elements 


Input  the  shortest 

element's  Radius  r8<^o  =  -00144 
(meters) 


Input  the  length 

of  the  shortest  *o  1  ^.219 
element  (meters) 


Input  the  distance 

from  shortest  to  d  1.659 

second  shortest 

element  (meters) 


Input  the  operating  r  .  <n6 

Frequency  (Hertz)  5  " 


Input  the  ground  g  ..  j.jq-3 
Conductivity 


Input  the  Characteristic 
Admittance  of  the  Line 
Feeding  the  Antenna 
(Ohms) 


ADM  = 


I 

300 


Input  the  length 

of  the  second  *1  "  8.297 

shortest  element 

(meters) 


Input  th^  height  of 

the  shortest  element  =  9.4 

above  the  surface 

(meters) 

Input  the  Distance 

from  the  Antenna  R  1  3000 

(meters)  1 

Input  the  ground  .4 

Dielectric  Constant 

T 


Input  the  Termination 
Impedance  connected 
to  longest  element  TIMP 
(Ohms) 


300 
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Input  the  angle  between  Elevation  index  for  which 

the  vertical  (z  axis)  and  ♦  =  90  to  calculate  horizontal  pattern  w  r  285 

the  antenna  axis  (degrees)  (from  angle  index  table) 

NOTE:  The  Horizontal  Log-Periodic  Dipole  mutual  impedance  calculations  are  not 
valid  if  there  is  an  antenna  element  whose  length  is  an  exact  integer  multiple  of  the 
wavelength.  If  this  occurs,  there  will  be  a  singularity  error  in  the  mutual  impedance 
calculations.  If  this  problem  arises,  it  will  be  necessary  to  vary  the  operating  frequency 
such  that  no  element  is  exactly  an  integer  multiple  of  the  wavelength. 


If  the  log-periodic  calculations  of  element  length  and  spacing  do  not  represent  the 
desired  antenna  configuration,  you  can  enter  the  values  directly  by  selecting  the 
variables  I  and  d,  using  the  define  key  (shift,  colon),  and  entering  the  successive 
lengths  and  spacings  by  separating  each  successive  entry  by  a  comma. 

Calculated  distance 

Calculated  length  between  successive  Calculated  radii 

of  elements  from  elements-shortest  of  elements  from 

shortest  to  longest  j0  longest  shortest  to  longest 


7.219 

8.297 

9.53598 

10.95997 

12.5966 

14.47762 

16.63954 

19.12429 

21.98009 

25.26233 

29.03471 

33.37041 


1.659 

1.90674 

2.19146 

2.51871 

2.89483 

3.32711 

3.82394 

4.394% 

5.05125 

5.80554 

6.67247 
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Radiation  Patterns  in  Phi=Pi/2  Plane 
(Perpendicular  to  X-Axis  and  Dipole  Elements) 
Space  Wave  Radiation  Pattern  (Phi=Pi/2) 


wet^.HtPj 

Max  E-Fidd  Intensity  (Volts  per  meter)  max(MagElP)  =  0.08512 
Surface  Wave  Radiation  Pattern  (Phi=pi/2) 


xref^ ,  m2Pj 

Max  E-Fie Id  Intensity  (Volts  per  meter)  max(MagE2P)  =  4.54242- 10  ^ 


Radiation  Patterns  in  Phi-=0  Plane 
(Perpendicular  to  Phi=pi/2  Plane) 

Space  Wave  Radiation  Pattern  (Phi=0) 


Max  E-Field  Intensity  (Volts  per  meter)  max(MagElO)  =  0.03495 


Surface  Wave  Radiation  Pattern  (Phi=0) 
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Radiation  Patterns  in  Horizontal  Plane 
(Parallel  to  Ground) 

Space  Wave  Radiation  Pattern  (Horizontal) 


Max  E-Field  Intensity  (Volts  per  meter)  max(MagE2H)  =  3.44012*  10 


1.5707963  r  1  5707963 


power#  = 


cos(O  cos(?)-Iy  .  j  p  Y  csc('f 
y  1  -  (sin(?).cos({))2 


.'-  1.5707963  .'0 


powers  = 


■:  Jo 
2 


sin(OIy 

1  -  (sin(0-cos(0)‘ 


J  P  Y  csc(  "F )  - ( cos(  r )  coj(  4 

t  1 


power  =  powers  +■  powers 


Radres  =  2- power 


E3P  ,  (|M^E3P,|)  ^  :  4-g-R*- m.Jt( sqE3P) 

'  2- (120-*)  power 


2  2 


sqEJOj  =  (l-M^3°i|)Z  DirectivilyO  ,  <*»2"«(».E30) 
1  2*(120-x)  power 


Total  Power  Radiated  (Watts) 


power  =  83.98604 


Radiation  Resistance  (versus 
Maximum  Input  Current) 
(Ohms) 


Radres  =  167.97209 


Directivity  (or  Maximum  Power  Gain  assuming  100%  Antenna  Efficiency) 
Phi  =  pi/2  Plane  Phi  =  0  Plane 


DirectivityP  =  12.93749 


Phi  =  0  Plane 
Directi  vityO  =  2.18175 


/ 


Effective  Isotropic  Radiated  Power  (EIRP)  (Watts) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 

DirectivityP*  power  -  1.08657*10^  DirectivityO*  power  -  183.23657 

Maximum  Effective  Area  (Along  Radial  of  Directivity) 

(square  meters) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 


(  •  DirectivityP 

4*x 


=  643.45795 


(Aj)  -DirectivityO 
4*x 


=  108.51135 


Maximum  Effective  Length  (Along  Radial  of  Directivity) 
(meters) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 


Radres*  (  Aj)  2  DirectivityP 


=  33.86438  2* 


Radres*  ( Aj)  2  DirectivityO 


480* x2 

Numerical  Distances 
Vertical  Polarization 


480  V 


PeQ|  =  69.28584 


Horizontal  Polarization 
3 


Pm0|  =  1.26447* 


10“ 


90659 


Elevation  Angle  of  Directivity  (Maximum  Gain) 
above  the  Horizon  (Degrees) 

Phi  -  pi/2  Plane  Phi  =  0  Plane 


AngleP  =  31.42857 


AngleO  = 


89.71429 

89.71429 
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HORIZONTAL  LOG  PERIODIC  DIPOLE  ARRAY 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated  with 
horizontal  log  periodic  dipole  arrays.  The  antenna  is  oriented  such  that  the 
projection  of  its  center  axis  lies  on  the  positive  y-axis  of  a  rectangular  coordinate 
system.  The  antenna  axis  can  be  oriented  with  respect  to  the  vertical  at  any  angle 
b.  en  zero  and  ninety  degrees.  The  feed  is  at  the  center  of  the  shortest  element 
directly  above  the  origin.  The  dipole  elements  are  bisected  by  the  antenna  axis  and 
are  parallel  to  the  x-axis.  Required  inputs  are  the  number  of  dipole  elements, 
shortest  and  second  shortest  element  lengths,  separation  between  the  shortest  and 
second  shortest  elements,  radius  of  the  shortest  element,  height  of  the  feed  above 
the  surface,  angle  of  the  antenna  axis  with  the  vertical,  characteristic  admittance  of 
the  line  feeding  the  antenna,  termination  impedance  of  the  antenna,  transmitted 
frequency,  distance  from  the  antenna,  the  conductivity  and  dielectric  constant  of  the 
surface  below  the  antenna.  The  planar  earth  model  is  assumed  in  predicting 
radiation  patterns.  A  sinusoidal  current  input  with  a  maximum  of  unity  is  assumed . 
All  radiation  patterns  are  normalized  with  respect  to  the  maximum  electric  field 
intensity  transmi'ted  by  the  antenna  in  the  plane  of  interest.  The  electric  field 
magnitudes  to  which  the  patterns  are  normalized  are  displayed  below  their 
respective  plots.  Radiation  patterns  are  plotted  for  the  phi— pi/2  and  phi=0  vertical 
planes  parallel  and  perpendicular  to  the  y-axis  respectively.  A  horizontal  radiation 
pattern  is  plotted  at  an  elevation  selected  by  the  index  from  the  elevation  angle 
index  table.  Polarization  is  a  varying  combination  of  horizontal  and  vertical 
depending  on  spatial  orientation  with  respect  to  the  antenna. 


Input  the  number  ^  ^ 

of  elements 

Input  the  length 

of  the  shortest  f()  =  ?219 

element  (meters) 


Input  the  distance 

from  shortest  to  -  1.659 

second  shortest 

element  (meters) 


Input  the  shortest 

element’s  Radius  rac*g  *  00144 
(meters) 

Input  the  length 

of  the  second  ^  =  ®  ^ 

shortest  element 

(meters) 

Input  the  height  of 

(be  shortest  element  Hq  =  9.4 

above  the  surface 

(meters) 


Input  the  operating  r  .  12-in* 
Frequency  (Hertz)  5  1 


Input  the  Distance 
from  the  Antenna 
(meters) 


R  •-  3000 


Input  the  ground 
Conductivity 


Input  (he  ground 
Dielectric  Constant 


Input  the  Characteristic  Input  the  Termination 

Admittance  of  the  Line  l  Impedance  connected 

Feeding  the  Antenna  ADM  "  ^qq  to  longest  element  TIMP  -  300 
(Ohms)  (Ohms) 


Input  the  angle  between  Elevation  index  for  which 

the  vertical  (z  axis)  and  ♦  =  90  to  calculate  horizontal  pattern  w  ‘  285 

the  antenna  axis  (degrees)  (from  angle  index  table) 

NOTE:  The  Horizontal  Log-Periodic  Dipole  mutual  impedance  calculations  are  not 
valid  if  there  is  an  antenna  element  whose  length  is  an  exact  integer  multiple  of  the 
wavelength.  If  this  occurs,  there  will  be  a  singularity  error  in  the  mutual  impedance 
calculations.  If  this  problem  arises,  it  will  be  necessary  to  vary  the  operating  frequency 
such  that  no  element  is  exactly  an  integer  multiple  of  the  wavelength. 


If  the  log-periodic  calculations  of  element  length  and  spacing  do  not  represent  the 
desired  antenna  configuration,  you  can  enter  the  values  directly  by  selecting  the 
variables  1  and  d,  using  the  define  key  (shift,  colon),  and  entering  the  successive 
lengths  and  spacings  by  separating  each  successive  entry  by  a  comma. 

Calculated  distance 

Calculated  length  between  successive  Calculated  radii 

of  elements  from  elements-shortest  of  elements  from 

shortest  to  longest  (Q  longest  shortest  to  longest 
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Radiation  Patterns  in  Phi=0  Plane 
(Perpendicular  to  Phi-pi/2  Plane) 

Space  Wave  Radiation  Pattern  (Phi=0) 


wet^ ,  hIOj 

Max  E-Field  Intensity  (Volts  per  meter)  max(MagElO)  =  0.05843 
Surface  Wave  Radiation  Pattern  (Phi=0) 


,  m20| 

Max  E-Field  Intensity  (Volts  per  meter)  max(  MagE20 )  =  4.46943*  10 
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Max  E-Field  Inteashy  (Volts  per  meter)  mex(  MagE2H)  =  5.572%*  JO 


1.5707963  r  1 .5707963 


powers  = 


1.5707963  -'0 


SC0$({)'C0$(n-J  jpY  CJC(1 

- 1 — e  Y 

1  1  -  (sin(f)  cos(?))2 


power0 


sin(0,I7 

(sin(f)-cos(O)2 


J  P-Y  cjctTJ^cosf;}  coj(1 

m  t 


power  -  power/?  *-  power# 


Radres  =  2  power 


sqE3P|  =  (|M^3F||)2  DireclivilyP  .  4' ■»:«*•■«( «q«P> 
2- (120-*)  power 


*,£30,  =  DirectivftyO  ”IE30> 

'  2(120*)  J  power 


Total  Power  Radiated  (Watts) 


power  =  261.928 


Radiation  Resistance  (versus 
Maximum  Input  Current) 
(Ohms) 


Radres  =  523.856 


Directivity  (or  Maximum  Power  Gain  assuming  100%  Antenna  Efficiency) 


Phi  =  pi/2  Plane 
DirectivityP  =  10.78579 


Phi  =  0  Plane 


DirectivityO  =  1.95488 
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Effective  Isotropic  Radiated  Power  (F.1RP)  (Watts) 

Phi  3  pi/2  Plane  Phi  =  0  Plane 

Directivity P-  power  =  2.8251-10^  DirectivityO' power  =  512.03765 

Maximum  Effective  Area  (Along  Radial  of  Directivity) 

(square  meters) 


Phi  =  ptf2  Plane 


(Xj)*  Directivity P 

4* 


*  536.4412 


Phi  =  0  Plane 


(X.')^1  DirectivityO 

SJZl - -  =  97.22772 

At 


Maximum  Effective  Length  (Along  Radial  of  Directivity) 
(meters) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 


Radres  (Xj)  •  DirectivityP 


c 


54.60488  2 


Radres-  (X j)  •  DirectivityO 


480V 


480- 1 

Numerical  Distances 

Vertical  Polarization  Horizontal  Polarization 


P*r 


0.06283 


Pm0|  =  2.2621  MO6 


*  23.24692 


Elevation  Angle  of  Directivity  (Maximum  Gain) 
above  the  Horizon  (Degrees) 


Phi  -  pi/2  Plene 


Phi  *  0  Plane 


AngleP  =  34.85714 


AngleO 


56.57143' 

56.57143/ 
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APPENDIX  Ft 

HORIZONTAL  YAGI-UDA  ARRAY  COMPUTER  OUTPUT 


V 


This  appendix  contains  computer  hardcopies  from  the 
Mathcad  horizontal  Yagi-Uda  array  application  which  show  the 
input  values  and  predicted  radiation  characteristics  for  . two 
sample  calculations.  The  configuration  is  a  given  by  the 
inputs  on  the  first  three  pages  of  each  printout.  The  first 
antenna  is  mounted  above  soil  (€,=4  and  cr^lO'3)  and  the  second 
above  seawater  (er«72  and  a»4)  .  Reference  6  [p.  107]  provides 
the  radiation  patterns  and  gain  predictions  from  several 
sources  for  the  configuration  in  the  first  example. 

The  radiation  patterns  and  maximum  directive  gain 
computed  by  the  first  Mathcad  example  are  almost  identical  to 
the  those  given  in  reference  6.  The  directivity  prediction  is 
slightly  higher  for  Mathcad,  but  the  overall  similarity  among 
the  predicted  radiation  patterns  is  very  good.  The  seawater 
example  has  results  very  much  like  the  soil  example.  For 
horizontal  polarization,  the  higher  conductivity  surface  below 
the  antenna  does  not  result  in  a  greatly  enhanced  surface  wave 
and  increased  directivity  as  it  does  for  vertical 
polarization. 
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HORIZONTAL  YAGI-UDA  ARRAY 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated  with 
horizontal  Yagi-Uda  arrays.  The  antenna  is  oriented  such  that  the  projection  of  its 
center  axis  lies  on  the  positive  y-axis  of  a  rectangular  coordinate  system.  The  feed 
is  at  the  center  of  a  single  driven  element  between  a  set  number  of  longer  parasitic 
reflector  elements  and  shorter  parasitic  director  elements.  Antenna  elements  are 
bisected  by  the  antenna  axis  and  are  parallel  to  the  x-axis.  The  center  of  the  first 
reflector  is  directly  above  the  origin  at  a  set  height  on  the  z-axis.  Required  inputs 
are  the  number  of  reflector  and  director  elements,  individual  element  lengths, 
individual  element  radii,  separation  between  adjacent  elements,  height  of  the  feed 
above  the  surface,  transmitted  frequency,  distance  from  the  antenna,  conductivity 
and  dielectric  constant  of  the  surface  below  the  antenna.  The  planar  earth  model  is 
assumed  in  predicting  radiation  patterns.  A  sinusoidal  voltage  with  a  one  volt 
maximum  across  (he  input  terminals  is  assumed.  All  radiation  patterns  are 
normalized  with  respect  to  the  maximum  electric  field  intensity  transmitted  by  the 
antenna  in  (he  plane  of  interest.  The  electric  field  magnitudes  to  which  the  patterns 
are  normalized  are  displayed  below  their  respective  plots.  Radiation  patterns  are 
plotted  for  the  phi=pi/2  and  phi=0  vertical  planes  parallel  and  perpendicular  to  the 
y-axis  respectively.  A  horizontal  radiation  pattern  is  plotted  at  an  elevation 
selected  by  the  index  from  the  elevation  angle  index  table.  Polarization  is  a 
varying  combination  of  horizontal  and  vertical  depending  on  spatial  orientation  with 
respect  to  the  antenna. 


NOTE:  The  Yagi-Uda  mutual  impedance  calculations  are  not  valid  if  there  is  an 
antenna  element  whose  length  is  an  exact  integer  multiple  of  the  wavelength.  If 
this  occurs,  there  will  be  a  singularity  error  in  the  mutual  impedance  calculations. 
This  should  not  be  a  problem  since  Yagi-Uda  Antennas  are  designed  to  operate  at 
frequencies  such  that  the  driven  element  is  approximately  one-half-wavelength  long, 
the  reflectors  are  slightly  longer,  and  the  directors  are  slightly  shorter.  If  the  problem 
does  arise,  it  will  be  necessary  to  vary  the  operating  frequency  such  that  no  element 
is  exactly  an  integer  multiple  of  the  wavelength. 


Input  the  number  of  j  Input  the  number  of 

reflector  elements  director  elements 

y  «  0„(NR  4-  ND)  0  «  0 .. ( NR  +■  ND  -  1 ) 


257 


This  entire  page  is  ailoted  for  these  entries  to  allow  for  as  many  elements  as  necessary. 


Input  the  lengths  and  radii  of  the  Yagi 
elements,  starting  with  the  outermost 
reflector  and  proceeding  inward  to  the 
driven  element.  After  entering  the  driven 
element  length  and  radius,  enter  those  of 
the  director  elements,  starting  with  the 
director  next  to  the  driven  element  and 
proceeding  outward  to  the  final  element. 

I  ■=  rad  ■= 

r  r 


15.6 

.001 

15  i 

.001 

14  ! 

i 

.001 

Input  the  distances  between 
successiv  elements  in  the  same 
order  as  the  length  and  radius 
inputs,  lliere  should  be  one  less 
entry  in  the  separation  distance 
array  than  there  is  in  the  length 
and  radius  array. 


7.5 

6.0 
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Input  (he  operating  f5  :=  10*  10^ 
Frequency  (Hertz) 

Input  the  ground  a  -  1  *  10“ 
Conductivity 


Input  the  height  of  Hq  7.5 
the  array  above  the 
surface  (meters) 

Input  the  ground  er  =  4 
Dielectric  Constant 


Input  the  Distance  ^  ..  jqqq 

from  the  Antenna 

(meters) 

Input  the  index  of  the  elevation 

angle  for  which  to  calculate  the  w  r  285 

horizontal  radiation  pattern 

(from  the  angle  index  table) 


Index  of  Refraction 


n  = 


er  -  j- 


18000-  o 

(V>o-s) 


Complex  Numerical  Distance  for  Vertical  Polarization 


*v= 


-j  g-(R  .  H„.coS(9,)) 
2-sin(0j)  2 


/ 


cos 


W) 


n2  -  sin 


in(*i)' 


Complex  Numerical  Distance  for  Horizontal  Polarization 


-j-g-(R  »  H0-cos(gt)) 
2-(sin(6|)  2) 


Vertical  Reflection  Coefficient 


Horizontal  Reflection  Coefficient 


cos(0|)  -  1 

h 

ln2-  sto(«,)2) 

cos(0,)  +  ( 

li”2  -  ‘b(9,)2) 
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Radiation  Patterns  in  Phi=Pi/2  Plane 
(Perpendicular  to  X-Axis  and  Dipole  Elements) 
Space  Wave  Radiation  Pattern  (Phi=Pi/2) 


w ,  HlP| 


Max  E-Field  Intensity  (Vohs  per  meter)  max (  MagElP)  58  5.13224- 1 
Surface  Wave  Radiation  Pattern  (Phi=pi/2) 


Max  E-Field  Intensity  (Volts  per  meter)  max(  MagE2P)  =  4.08957*  10  ^ 
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Combined  Space  and  Surface  Wave  Radiation  Pattern  (Phi-pi/2) 


Max  E-Field  Intensity  (Vohs  per  meter)  max  (  MagE3P )  *  5.13132*  10 
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Radiation  Patterns  in  Phi=0  Plane 
(Perpendicular  to  Phi=pi/2  Plane) 

Space  Wave  Radiation  Pattern  (Phi=0) 


lax  E-Field  Intensity  (Volts  per  meter)  max(MagElO)  =  1.35264-11 
Surface  Wave  Radiation  Pattern  (Phi=0) 


E-Field  Intensity  (Volts  per  meter)  max(MagE20)  =  5.32341-10 
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Combined  Space  and  Surface  Wave  Radiation  Pattern  (Phi=0) 


-1  -0.5  0  0.5  1 

Hret^ ,  h30j 

Max  E-Field  Intensity  (Volts  per  meter)  max(MagE30)  =  1. 35264 •  10 


?63 


Max  E-Field  Intensity  (Volts  per  meter)  max(  MagE2H)  =  2.73697*  10 


powerfl  - 


30 

*R: 


1.5707963 


-  1.5707963 


1.5707963 


0 


cos(*)  cos(O  I  j  p  D  sin(«: 

_ r 


1  -  (sin(f)  cos(O)2 


power#  •= 


30 

*R5 


it 

2 


.*  -'0 
2 


sin(i)I 


r 


1  -  (sin(f)-cos(O)' 


JP  D !  sfn(C)  sin(U 


power  =  power#  *  power# 


sqE3P,  = 


( |  MagE3Pj  [ ) ' 


2(  120*) 


Radres  - 


2  power 


DiredivityP  ■°»(sgE3P) 

power 


sqE30{  - 


DirectivilyO  .  I*  *2'  -»«( sqE30) 
2’(  120'*)  power 


Tolal  Power  Radiated  (Watts) 

Radiation  Resistance  (versus 
Maximum  Input  Current) 
(Ohms) 


power  =  0.00266 
Radres  =  187.94973 


Directivity  (or  Maximum  Power  Gain  assuming  100%  Antenna  Efficiency) 
Phi  =  pi/2  Plane  Phi  =  0  Plane 

Directivity  P  =  14.84639  Directi  vityO  =  1.03163 
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HORIZONTAL  YAGI-UDA  ARRAY 


This  application  calculates  far  field  radiation  patterns  and  parameters  associated  with 
horizontal  Yagi-Uda  arrays.  The  antenna  is  oriented  such  that  the  projection  of  its 
center  axis  lies  on  the  positive  y-axis  of  a  rectangular  coordinate  system.  The  feed 
is  at  the  center  of  a  single  driven  element  between  a  set  number  of  longer  parasitic 
reflector  elements  and  shorter  parasitic  director  elements.  Antenna  elements  are 
bisected  by  (he  antenna  axis  and  are  parallel  to  the  x-axis.  The  center  of  the  first 
reflector  is  directly  above  the  origin  at  a  set  height  on  the  z-axis.  Required  inputs 
are  the  number  of  reflector  and  director  elements,  individual  element  lengths, 
individual  element  radii,  separation  between  adjacent  elements,  height  of  the  feed 
above  the  surface,  transmitted  frequency,  distance  from  the  antenna,  conductivity 
and  dielectric  constant  of  the  surface  below  the  antenna.  The  planar  earth  model  is 
assumed  in  ptedicling  radiation  patterns.  A  sinusoidal  voltage  with  a  one  volt 
maximum  across  the  input  terminals  is  assumed.  All  radiation  patterns  are 
normalized  with  respect  to  the  maximum  electric  field  intensity  transmitted  by  the 
antenna  in  the  plane  of  interest.  The  electric  field  magnitudes  to  whirh  the  patterns 
are  normalized  are  displayed  below  their  respective  plots.  Radiation  patterns  are 
plotted  for  the  phi=pi/2  and  phi=0  vertical  planes  parallel  and  perpendicular  to  the 
y-axis  respectively.  A  horizontal  radiation  pattern  is  plotted  at  an  elevation 
selected  by  (he  index  from  the  elevation  angle  index  table.  Polarization  is  a 
varying  combination  of  horizontal  and  vertical  depending  on  spatial  orientation  with 
respect  to  the  antenna. 


NOTE:  The  Yagi-Uda  mutual  impedance  calculations  are  not  valid  if  there  is  an 
antenna  element  whose  length  is  an  exact  integer  multiple  of  the  wavelength.  If 
this  occurs,  there  will  be  a  singularity  error  in  the  mutual  impedance  calculations. 
This  should  not  be  a  problem  since  Yagi-Uda  Antennas  are  designed  to  operate  at 
frequencies  such  that  the  driven  element  is  approximately  cavkalf- wavelength  long, 
the  reflectors  are  slightly  longer,  and  the  directors  are  sligbuy  shorter.  If  the  problem 
does  arise,  it  will  be  necessary  to  vary  the  operating  frequency  such  that  no  element 
is  exactly  an  integer  multiple  of  the  wavelength. 


Input  the  number  of  jqp  .  j  Input  the  number  of 

reflector  elements  director  elements 

y  «  0  ..  (  NR  4-  ND )  5  *  0 ..  (  NR  +  ND  -  1 ) 
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This  entire  page  is  alloted  for  the «*  entries  to  allow  for  as  many  elements  as  necessary. 

Input  the  lengths  and  radii  of  the  Yagi  Input  the  distances  between 

elements,  starting  with  the  outermost  successiv  elements  in  the  same 

reflector  and  proceeding  inward  to  the  order  as  the  length  and  radius 

driven  element.  After  entering  the  driven  inputs.  There  should  be  one  less 
element  length  and  radius,  enter  those  of  entry  in  the  separation  distance 

(he  director  elements,  starting  with  the  array  than  there  is  in  the  length 

director  nexl  to  the  driven  element  and  and  radius  array, 

proceeding  outward  to  the  final  element. 
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Input  (he  operating  fj  *  10- 10^ 
Frequency  (Hertz) 


Input  the  ground  o  -  A 
Conductivity 


Input  the  height  of  Hq  -  7.5 
the  array  above  the 
surface  (meters) 

Input  the  ground  cr  -  72 
Dielectric  Constant 


Input  the  Distance  p  _ 
from  the  Antenna 
(meters) 

Input  the  index  of  the  elevation 

angle  for  which  to  calculate  the  w  r  285 

horizontal  radiation  pattern 

(from  (lie  angle  index  table) 


Index  of  Refraction 


.  I8OOO0 

n  -  zr  -  i  -r - - 

J  (f5' 10  ‘) 


Complex  Numerical  Distance  for  Vertical  Polarization 


-j  g  (R  .  % co,^)) 
2-sin(0j)2 


cost  9.)  » 


M9\Y 


Complex  Numerical  Distance  for  Horizontal  Polarization 

Pl"i  ’  *  ft  ■  | 


Vertical  Reflection  Coefficient  Tv, 


s  (n2-cos(g,))  -  (jn2  -  sin(g,)2) 
(n2cos(0,)J  ♦  (Jn2^  sin(0t)2) 


Horizontal  Reflection  Coefficient 


c°s(g|)  -  (Jn2  -  sin(0,)2j 
cos(e,)  *  (Jo2  -  sin(0,)2) 
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Radiation  Patterns  in  Phi=Pi/2  Plane 
(Perpendicular  to  X-Axis  and  Dipole  Elements) 
Space  Wave  Radiation  Pattern  (Phi=Pi/2) 
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Radiation  Patterns  in  Horizontal  Plane 
(Parallel  to  Ground) 

Space  Wave  Radiation  Pattern  (Horizontal) 


-t  0.5  0  0.5  1 


.  h1H. 

Max  E-Field  Intensity  (Volts  per  meter)  max(MagElH)  =  4.25018*  10  4 


Surface  Wave  Radiation  Pattern  (Horizontal) 


?  75 


power#  - 


I*  1.5707963 


-  1.5707963 


1.5707963 


JO 


cos(()cos(OIT  j-p-D 

I  -  (sin(?)  cos(  O)2 


•  *h>( 


power  =  0.00482 
Radres  =  103.76004 

Directivity  (or  Maximum  Power  Gain  assuming  100%  Antenna  Efficiency) 
Phi  =  pi/2  Plane  Phi  =  0  Plane 

DirectivityP  =  14.63143  DirectivKyO  =  1.67913 


Total  Power  Radii,  ed  (Watts) 

Radiation  Resistance  (versus 
Maximum  Input  Current) 
(Ohms) 
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>l> 


Effective  Isotropic  Radiated  Power  (EIRP)  (Watts) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 

DirectivityP- power  =  0.07051  DirectivityO- power  =  0.00809 

Maximum  Effective  Area  (Along  Radial  of  Directivity) 
(square  meters) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 


(  X5)  2-  DirectivityP 

4* 


1.0479*  103 


(Xg)  2  DirectivityO 

4* 


120.25911 


Maximum  Effective  Length  (Along  Radial  of  Directivity) 
(meters) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 


2- 


Radres-  ( X«)  2  DirectivityP 

- -  =  33.%559  2 

480-x2 

Numerical  Distances 

Vertical  Polarization  Horizontal  Polarization 

|ft0|  =  0.04363  |  PmQ|  =  2.26206*  10* 


Radres-  (Xg)  2-  DirectivityO 
480- % 2 


11.50637 


Elevation  Angle  of  Directivity  (Maximum  Cain) 
above  the  Horizon  (Degrees) 

Phi  =  pi/2  Plane  Phi  =  0  Plane 


AngleP  =  41.71429 


AngleO  = 


89.71429 

89.71429 
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